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EISENHOWER AND GRASS RIVER LOCK MODELS# 


Martin E. Nelson, | M. ASCE, and 
Harvey J. Johnson,2 A.M. ASCE 
(Proc. Paper 1582) 


SYNOPSIS 


Models of the St. Lawrence Seaway Dwight D. Eisenhower and Grass River 
cks were constructed to a scale of 1 to 24.24, model to prototype. A type of 
ydraulic system consisting of intake ports in the upper gate sill, culverts 
ad ports in both walls of the lock and lateral culvert diffusers in the lower 
sproach channel was selected prior to model studies. Vortical activity over 
ie intake ports and turbulence in the lock chamber and lower approach were 
duced to acceptable proportions. Culvert and gate pressures, emergency 
ate and tainter valve operating loads, and hawser forces were measured. 
rototype operation procedures were selected from results of the model 
sts. 


INTRODUCTION 


Model studies of the Dwight D. Eisenhower and Grass River locks, St. 
awrence Seaway, have been conducted by the U. S. Army Engineer District, 
- Paul, in facilities provided by the University of Minnesota at the St. 
anthony Falls Hydraulic Laboratory. The locks and other features of the 
merican portion of the seaway were designed for the St. Lawrence Seaway 
=velopment Corporation by the U. S. Army Engineer District, Buffalo. Dur- 
gthe course of the model tests several conferences were held at the 


te: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1582 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 

Paper presented at a meeting of the American Society of Civil Engineers, 
June, 1957, Buffalo, N. Y. 

Hydraulic Engr., St. Paul Dist., Corps of Engrs., U. S. Dept. of the Army, 
St. Paul, Minn. 

Hydraulic Engr., St. Paul Dist., Corps of Engrs., U. S. Dept. of the Army, 
St. Paul, Minn. 
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laboratory with representatives of the Corporation; Office, Chief of Engi- 
neers; U. S. Army Engineer Division, North Central; and U. S. Army Engi- 
neer District, Buffalo. Mutually helpful conferences were also held by 
Corps laboratory personnel and engineers of the St. Lawrence Seaway Au- 
thority of Canada. 


Criteria for Design of Lock Hydraulic Systems 


The obvious function of a navigation lock is to facilitate the transit of 
water-borne traffic past abrupt changes in elevation of navigable rivers or 
canals. The design engineer holds as his primary objective a structure that 
will perform this function rapidly, safely and economically, insofar as these 
characteristics are compatible and mutually attainable. The key is found in 
the design of a sound and efficient filling and emptying system for the lock. 
The first criterion by which the performance of the hydraulic system is 
gaged is the time required to elevate a vessel from the lower pool level to 
the surface of the upper pool or the reverse of this process in the case of 
downbound traffic. Increasingly greater emphasis is being placed-on the im- 
portance of expediting river traffic along navigation channels and the avoid- 

_ance of unnecessary delays in lockage. A monetary value of time gained or 
sacrificed in lockage, based on the operating costs of towing equipment and 
the value of cargo, is sometimes used to determine the amount of capital in= 
vestment that would be justified to achieve a more efficient lock operation. — 
The second criterion requires that the forces acting on a vessel during a 
lockage operation shall not produce excessive stresses in the mooring lines 
nor endanger the safety of the craft or any part of the structure. This cri- 
terion is to some extent in conflict with the first. Obviously, as the capacity 
of the hydraulic system increases the amount of energy to be dissipated in- 
creases and more elaborate devices are required to reduce to satisfactory 
proportions the attendant disturbances in the lock chamber. The evidences 
of unsatisfactory hydrodynamic performance of a lock filling and emptying 
system are manifested in excessive strains in the mooring lines due to 
strong currents or surging of the water level in the lock chamber and exces 
sive turbulence. Hawser forces are most critical with respect to large ves- 
sels while unrestrained turbulence can cause hazardous conditions for smal- 
ler craft. The third criterion to be held in proper perspective together wi he 
speed and safety is that of economy. In his efforts to achieve a high degree 
of efficiency and safety in the operation of the hydraulic system the designing 
engineer must not lose sight of the fact that construction costs generally in-— 
crease in proportion to the refinement in hydraulic design. Consequently, 
the importance of the respective criteria must be weighed by the designer 
whose judgment and experience must to a large extent enter into the selectio 
of the optimum type of hydraulic system for a specific lock project. 


Selection of Hydraulic System Type 


Lock filling and emptying systems in use today may be classified into 
three general types: (1) End filling systems admit the lockage water near 
the upstream end of the lock chamber by means of loop culverts around the — 
service gates, sluices through the sill, or by operation of the gates. A simi- 
lar arrangement at the downstream end provides the means for emptying the 
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ck. (2) The most common type of hydraulic system used in locks of nomi- 
1 height consists of culverts in the lock walls with intake ports in the upper 
ol and short ports to distribute flow into the lock chamber. In the empty- 
g operation the flow in the chamber ports is reversed and discharges into 
e lower pool through port manifolds or discharge basins. (3) Several locks 
ve been constructed in recent years with filling systems that admit the 
ckage water into the lock chamber through openings in the floor. The sup- 
y conduits with intake manifolds in the upper pool are located either in the 
ck walls or in the floor. From the culverts the flow enters the chamber 
rough a system of lateral and port manifolds or directly through vertical 
rts as in the case of longitudinal floor culverts. The emptying operation 

r this type of system is performed with any of the various types used with 
ull port hydraulic systems. 

The end filling type of hydraulic system was not considered practical for 
€ seaway locks because prior experience had indicated that satisfactory 
ckage conditions could not be attained in locks of the size in question with- 
it using excessively long periods of time to fill and empty the chambers. 
1e floor lateral system undoubtedly would have met all the design criteria 
id provided satisfactory lockage conditions for the Eisenhower and Grass 
ver locks. However, the wall port system generally can be constructed at 
substantial saving in cost and the fortunate circumstance of a deep water 
Shion made available by the necessity for removal of overburden to a con- 
derable depth in order to obtain adequate bearing for the lock walls made it 
pear possible to satisfy the design criteria with this type of system. In 
sure 1 are shown data obtained from several tests in prototype locks with 
ull port hydraulic systems which provide a guide as to permissible rates of 
ling under various depths of water cushion in the lock chamber. These 
ta, which have been used in preliminary design of several new locks, indi- 
ted that the seaway locks with wall-port hydraulic systems could be ex- 
cted to operate satisfactorily. 

Model tests on locks have repeatedly indicated that hydrodynamic forces 
erted upon craft in the chamber during filling or emptying operations are 
e to swelling of the water surface along the lock axis. Swells or surges are 
t up at the very beginning of the valve opening period when the initial pres- 
re wave is transmitted along the culvert and into the lock chamber through 
= connecting ports or laterals. The degree of surging and consequently the 
agenitude of force imposed upon a craft in the lock chamber reaches a mini- 
am value when the manifold system of ports or laterals occupies the ap- 
oximate middle third of the lock chamber walls. This feature was approxi- 
ated in the design of the seaway locks. 


General Description and Purpose 


The Eisenhower and Grass River locks located in the Long Sault Canal 
ovide a means for navigation to pass the Long Sault Dam and power plant. 

e locks are 80 feet wide and 860 feet long with a 30-foot depth over the gate 
lis. The Eisenhower lock has a normal lift of 43 feet and the Grass River 
*+k located downstream is designed for a normal lift of 49 feet. The inter- 
sdiate channel or pool between the two locks is about 3 miles long and has a 
vigable width of 442 feet. 

A vertical lift gate is located in the upper end of the Eisenhower lock as an 
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emergency closure structure. The service gates of both locks are of the 
two-leaf mitering type. The hydraulic systems for the locks are almost 
identical consisting of a culvert in each wall drawing water from the upper 
approach channel through ports in the gate sill, filling or emptying the cham 
ber through 20 ports in each wall and discharging water from the locks 
through lateral culvert diffusers in the floor of the lower approach channel. 

The model studies described in this paper were conducted for the purpose 
of developing filling and emptying systems for the Eisenhower and Grass 
River locks that will perform the lockage operations as rapidly as possible 
without violating the requirements for safety of vessels and structure. The 
criteria adopted for the design of the hydraulic systems specified a lockage ~ 
time of about 8 minutes with a maximum stress not to exceed 10 tons in the 
mooring lines of a vessel of any size and a degree of turbulence that would 
not jeopardize small craft in the chambers or in the entrances to the locks, 
The scope of the model studies was broadened during the course of the pro- 
ject design stage to obtain data to aid in the design of the vertical lift emer-— 
gency gate for the Eisenhower lock, the purpose of which is to prevent loss ~ 
of the power pool in the event of damage to the lock service gates. Other 
incidental observations were made in the model to determine whether operat: 
ing or maintenance problems in the nature of cavitation, vibration, or erosion 
might be indicated and to measure the flood discharge capacity of the hy- 
draulic systems. 

The similarity of the two locks made it possible to conduct the tests in on 
composite lock model making minor modifications to convert from physical 
properties of one lock to those of the other. The lock models were construct: 
ed to a scale 1/24.24 of the proposed prototypes. However, dimensions and 
quantities used in this paper refer to the prototype. 


Intakes 


The intakes in the sill for the wall culverts, figure 2, are two manifolds, 
each consisting of three ports, 9 feet 6 inches by 22 foet at the top convergin 
to a constricted throat area. Uniformity of flow among ports was accom- 
plished by making the constricted area progressively smaller as the flow 
distances diminished. The maximum variation from uniformity of flow in 
intake ports was 2 percent. 

Vortices over the intakes presented a problem in the first model when the 
port dimensions were only 8 feet by 12 feet. The proximity of the miter gate = 
and gate recesses in the walls was considered a contributing factor but satis- 
factory conditions were obtained by doubling the port area and enlarging the 
constricted sections. This reduced the maximum port velocity from 21 feet 
per second in the first model to 12 feet per second in the final model. All 
pressures measured in the water passages were found to be positive. 


Culvert Bends 


Flow from the intake manifolds passes through a 90-degree bend in each 
culvert. Two piezometers are located at the midpoint of each bend in oppo- 
site walls of the culvert to measure pressures and serve incidentally as flow 
metering devices for the purpose of comparing calibrations of the model and 
prototype. Pressures measured for the adopted filling operation were alway 
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the positive range at both piezometers. The formula derived for the flow 
eters was: 


Q = 1230x _/(pressure difference between opposite piezometers, c.f.s. 
in feet of water) 


Valves 


The valves which control the flow in the filling and emptying culverts are 
the reversed tainter type, that is, with the trunnions upstream of the skin 
ate. The linkages between constant speed motors and the valves produce an 
ening pattern that is not a straight line function of time; the valve opening 
te increases slightly with respect to time. The culverts at the valves are 
feet wide and 14 feet high. Stop log recesses located about 20 feet up- 
ream and 30 feet downstream from the valves are provided to facilitate 
watering for inspection and maintenance. 

Culvert pressures downstream from the valves momentarily reached a 
ightly negative pressure at the halfway point of the valve opening cycle. 
ie stop log recess in that region will be sealed to prevent the possibility of 
training air in the flow which may accumulate in the culvert and then be re- 
ased suddenly into the chamber with such force as to endanger small craft. 


Chamber Manifold 


Between the valves and the lock ports the culvert is enlarged from 12 feet 
14 feet to 15 feet by 16 feet. 

Water enters the lock through 20 venturi-shaped ports in each wall which 
e 3 feet 1 inch wide at the throat, 3 feet 7 inches wide at the wall face, and 
eet high and 15 feet long. All the ports are sloped downward 1 on 15 and 

> upstream 10 ports in each wall are turned 3.5 degrees in the upstream 
rection, the remainder being perpendicular to the lock wall. They are sub- 
srged about 34 feet below normal tail water and are about 8 feet above 
tural rock floor in the Grass River lock and about 26 feet above natural 

ck floor in the Eisenhower lock. The average jet velocity in the ports 
aches a maximum of about 20 feet per second. 

The limitation of surface turbulence was difficult in the Grass River lock 
cause of the relatively shallow water depth in the chamber when the inflow 
ring filling is the greatest. This condition is shown graphically on figure 1 
ere the solid line is an approximate envelop curve of acceptable turbulence 
‘the chamber water depth and rate of filling in a number of existing locks. 
te that the rate of filling for the Grass River lock exceeds that given by the 
velop curve. This occurs for about 3 minutes of the 8 1/2-minute filling 
1e. Experimentation in the model resulted in a design that reduced the 
bulence to acceptable proportions. 

The ports were first tested in a single port model and then a refinement of 
sign was accomplished in the composite lock model of the Grass River lock. 
e basic designs tested are shown in figure 3. The 1 on 15 port slope, flare 
the walls, number and size of ports were selected as giving the best water 
‘face conditions. The 3.5-degree horizontal angle was adopted because this 
ture indicated a reduction in hawser pulls. For stability reasons the lock 
ll was designed with a toe projecting into the chamber near the base of the 
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wall. A sloping toe was sufficient to satisfy the structural design but the I 
model tests indicated a slight reduction in turbulence during filling opera- 
tions by building the projection in steps. 

A great number of baffles, deflectors and port shapes were tested but the” 
only modification which performed as well as the design developed was an ~ 
interim test port 5 feet 2 inches square at the large end and 23 feet long, 
projecting horizontally 8 feet beyond the lock wall. The jet issuing from this 
port entered the lock on a slope of 1 on 15. The same port length to width — 
ratio and jet slope were approximated in the final port design. 


Discharge Manifold 


The emptying system is a continuation of each wall culvert of the filling © 
system which terminates in a manifold of three smaller lateral culverts 11 — 
feet 6 inches high tapered from 6 feet 4 inches wide to 2 feet in width which 
diffuse flow in the lower approach channel through a total of 96 ports 4 feet 
3 inches high and 1 foot 8 inches wide. The lateral culvert walls were 2 feet 
thick. Thus, the ports were 2 feet long but a remarkable improvement in 
water surface conditions was obtained by doubling the length of the ports by 
means of projections on the culvert walls. Velocities along the bottom of the 
lower approach channel reached a maximum of about 5 1/2 feet per second. 


Lockage Operations 


Model hulls of two ships were used in this study. The George S. Humphrey 
is considered one of the largest carriers on the Great Lakes with a total { 
length of 710 feet and a 75-foot beam, representing a gross displacement of 
37,000 tons when loaded to a maximum draft of 27 feet. The George Eaton is 
a smaller ship, length 257 feet, beam 43 feet, draft 16.5 feet, and displace- 
ment 5,000 tons. The forces required to hold these ships at various location 
in the model lock chamber were recorded for a wide range of lifts and valve 
operations. 

A recorder built especially for this type of testing made a continuous 
record of water surface levels, valve movement, surging, and hawser pulls. 
These data are shown on figures 4 to 7 for the selected operation of the two 
locks. 

To avoid exceeding a maximum longitudinal hawser pull of 10 tons by the 
largest ship a valve opening time of about 4 minutes was selected for filling 
and 2 minutes for emptying. A lifting force of about 20 tons more than fric- 
tion and weight will be required to raise the valves. Vibration of the model 
valve occurred during the first foot of opening and reached a peak frequency 
of 1/5 cycle per second producing variations in loading of between 2 tons 
12 tons, in prototype values. 

The Eisenhower lock will fill in 8.1 minutes with a valve opening time of 4 
4 minutes and will empty in 8.0 minutes with a valve opening of 2 minutes, 3 
while operating at a maximum lift of 43 feet. 

The Grass River lock will fill in 8.5 minutes when the valves are opened 
in 4 minutes and will empty in 8.3 minutes when the valves are opened in 2 
minutes, while operating at the maximum lift of 49 feet. 
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Emergency Gate 


The emergency gate in the upper approach of the Eisenhower lock, shown 
in figures 2 and 8, consists of a skin plate supported on the downstream side 
by horizontal plate girders spanning the 80-foot channel. The bearings of the 
gate roll on the downstream side of a recess in each wall and the gate when 
open will be submerged in a pit in the bottom of the channel. 

The gate model simulated 50 feet of the 80-foot channel width to a scale of 
3 to 100. 

The gate as designed with a sharp crest functioned satisfactorily. The 
force required to raise it varied from 364 tons, the submerged weight while 
resting in the recess of the channel bottom, to 880 tons near the top of the 
closure operation. The three conditions of flow and total lifting forces are 
illustrated in figure 8. The nappe becomes aerated as the crest reaches 
elevation 226 when air is drawn in through the wall recesses. The nappe 
falls on the top plate girder when the crest is near elevation 239. The loss 
of buoyancy when the nappe becomes aerated and the weight of water trapped 
on the webs of the girders are also shown on figure 8. No attempt was made 
to simulate prototype friction factors in the model and consequently the 
friction loads indicated by the model are excessive relative to those antici- 
pated in the prototype. 
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DEPTH OF WATER CUSHION IN LOCK CHAMBER - FT. 
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RATE OF RISE OF WATER SURFACE, FT./ MINUTE 


FIG. |1—- COMPARISON OF FILLING RATES 
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FIG. 5- EMPTYING CHARACTERISTICS 
EISENHOWER LOCK 
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FIG. 6— FILLING CHARACTERISTICS 
GRASS RIVER LOCK 
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FIG. 7- EMPTYING CHARACTERISTICS 
GRASS RIVER LOCK 
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HERMAL DENSITY UNDERFLOW DIVERSION, KINGSTON STEAM PLANT 


Rex A. Elder,! M. ASCE and Gale B. Dougherty,2 A.M. ASCE 
(Proc. Paper 1583) 


SYNOPSIS 


During the design of the Tennessee Valley Authority’s Kingston Steam 
lant it became desirable to develop facilities which would divert cold waters 
hich pass the site as thermally stratified density underflows and cause them 
} be drawn approximately 2 miles up a side channel at the same depth as the 
riginal channel, after which they would cross an area one mile long and 20 
) 25 feet higher than the original channel. The unique solution to this prob- 
‘m consisted of a canal across the high ground area at the end of which was 
ced a skimmer wall structure which allowed only the cold waters in the 
xttom of the side channel to flow into the canal. A diversion dam was also 
aced underwater immediately downstream from the side channel. This pa- 
sr gives the methods used in determining the required design. Much of the 
udy concerned density current flows where analytical methods were not 
yailable and had to be developed. The efficacy of these methods and the 
lopted design is proven by the presentation of field operational data obtained 
ter the structures were placed in use. The financial benefits attributable 
j}the design were found to be high. 


INTRODUCTION 


During the construction of the Kingston Steam Plant, TVA engineers were 
ced with a challenging problem in developing a design which would make 
railable for condensing purposes cold waters which passed the plant site as 


te: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1583 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 
Presented at a meeting of the American Society of Civil Engineers, June, 
1956, Knoxville, Tenn. 
Head, Hydr. Lab., Div. of Water Control Planning, Hydr. Data Branch, 
Tennessee Valley Authority, Norris, Tenn. 
Hydr. Engr., Southern Services, Inc., Birmingham, Ala. 
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ermal density underflows. The general site plan for the Kingston Steam 
lant is shown in Figure 1. The plant is located on the neck of a peninsula 
rmed by the Emory and Clinch River embayments of the Watts Bar Reser- 
ir. During the summer months the only appreciable flows in the reservoir 
-cur as thermal density underflows in the Clinch River arm.(1) The source 
these underflowing waters is TVA’s Norris Reservoir some 78 miles up- 
ream. Because of stratification in this large storage reservoir, relatively 
id water is discharged through the powerhouse turbines from spring until 
te fall. Over-all steam plant considerations required that the condensing 
aters be drawn from the north or Swan Pond side of the peninsula and dis- 
1arged on the south or Clinch River side. Under normal conditions the 
mory River flows are negligable during the summer months. As a conse- 
lence the cold Clinch underflows push their way up the original Emory River 
lannel. This fact is shown by Figure 2, which presents the maximum tem- 
sratures observed at Clinch River Mile 4.0 and Emory River Mile 2 over a 
1/2-year period. Except for the surface the two sets of curves are quite 
milar, indicating that the Clinch River waters penetrate and thus control 

e Emory River temperatures in the immediate area of the Kingston Plant. 
igure 2 also shows that in the Swan Pond Embayment the surface tempera- 
res become very hot, reaching 95° F. on occasions. 

Figure 3 shows a section taken along the southern edge of Swan Pond at 
sction A-A of Figure 1. This section is typical for Swan Pond and shows 
at the water depth in this area is comparatively shallow and the water is 
us susceptible to stagnation and marked temperature increases due to solar 
ating. In its natural state, the high bottom elevation of this area in essence 
rms a dam which would block out the cold underflowing Clinch River water. 
nder this condition, only the hot surface waters would be drawn into the 
imps. 

For the best over-all plant design it was desirable to set the maximum 
lowable condenser water temperature at 75° F. The plant as originally 
mceived was to have a maximum condenser water demand of 2470 cubic 
et per second. This value was used for all underflow diversion design pur- 
ses. Later changes reduced the maximum demand to 2310 cubic feet per 
‘cond. There was no question about the normal Clinch flows being sufficient 
id of the proper temperature to meet the plant demand. The problem was 
us one of developing a design which would economically make available to 
€ pumps the cold underflowing Clinch River waters by bringing these up the 
mory River channel and across Swan Pond to the condenser water pump in- 
kes. 


Proposed Design 


One obvious solution to the problem would have been to excavate a channel 
ross Swan Pond so that the cold Clinch waters could be reached. This solu- 
In was not economical since over 1 1/2 million cubic yards of material, a 
eh percentage of which would have been rock, would have had to be exca- 
ted under water. It was proposed, therefore, to provide an intake canal 
ross the high Swan Pond area and into the original Emory River channel by 
eans of dikes such as shown in Figure 1. A wall structure, as shown on 
gures 1 and 3, would close the east end of the canal except for a series of 
enings along its bottom. These openings would allow the cold bottom water 
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flow into the canal to replace water pumped from the canal for condensing : 
‘poses. In effect the wall would skim off the hot top waters and it was 
refore termed a skimmer wall. Figure 4 shows the details of the final 
ll design. Figure 5 is a photograph of the finished structures. 
In addition, it was desired to investigate the feasibility of an underwater 
n located on the Clinch River just downstream from the Emory River 

uth. Since the releases from Norris Dam are intermittent and subject to 

reme regulation,(1) especially during the summer months, it appeared 

t such a dam might be required to divert the Clinch flow into the Emory 

innel at a high enough level to make full use of the available cold water. 


sign Considerations 


The design required three separate but interrelated studies: (1) The study 
flow conditions in the Emory River channel from its mouth to the skimmer 
1, (2) the study of the diversion of underflow waters into the Emory from 
Clinch, both with and without an underwater dam, and (3) the determina- 

n of the basic dimensions for the skimmer wall openings. The first study 
ermined the capability of the Emory River channel, and when combined 

h the second study it determined the anticipated temperatures. The 

‘ond study also established the effectiveness of an underwater dam. 


jw in Emory River Channel 


At the time of these studies, TVA had a large number of field measure- l 
nts of the flow of thermal density currents obtained in Watts Bar and Fort 

idoun Reservoirs. An analysis of these data indicated that, in reservoir 

mnnels such as those at the project site and for thermal stratifications of 
magnitude expected at the project site, a relationship exists between the 

ximum velocity at which density underflow can occur and the depth of the 

lerflow. Figure 6 shows the data and the envelope curve established for 

1 in the Kingston study. This curve indicates that, if at any given depth the 

ocity is equal to or less than the value given by the curve, thermal under- 

w can be sustained. 

A second item of data available from the TVA field measurements was 

t, at water intakes, the velocity in the intake could be double the underflow 

ocity without breaking the stratification and pulling in warmer overlying 

ers. Velocities greater than double appeared to break the stratification. 

By using Figure 6 the capability of the Emory River channel was deter- | 
1ed in a preliminary manner. The section of minimum cross-sectional | 
a in the channel was determined and a curve of depth versus velocity was i 
tted for the maximum condenser demand. This curve, cross-plotted on 

ure 6, showed that a depth of 23 or more feet would be required to sustain 

ensity underflow and the velocity at this depth would be 0.35 foot per 

ond. The river channel bottom at the point of minimum cross-sectional 

a was at elevation 704. Therefore the maximum requirement of water 

ld be drawn from below elevation 727. Reference to the Clinch River 

e 4 temperature curves of Figure 2 makes it apparent that an average 

perature of 75° F. or less would always exist throughout the depth be- 

en elevations 704 and 727. 

A more comprehensive survey of the capability of the Emory River chan- 

was made after preliminary dimensions for the skimmer wall openings 

e established. In this study the density underflow interface profiles were 


Figure 4.--SKIMMER DETAILS--HALF PLAN AND ELEVATION 


KINGSTON STEAM PLANT 


‘SCE 


1583-8 HY 2 April, 195 | 


ree ia te te 


O 


) 0.5 
VELOCK Tiare: S 


Figure 6.--CRITERION FOR THERMAL DENSITY UNDERFLOW 


SCE KINGSTON STEAM PLANT 1583-9 


smputed by a modified standard step backwater computation method. A 
umber of authorities have shown that stratified flow calculations can be 
jade by using normal hydraulic equations if the gravitational term for the 
nderflow is multiplied by the factor 4” where ? is the density of the under 
ow fluid and A? is the difference in density between the underflow and the 
verlying layer. 
For the Emory River channel backwater computations the Chezy equation 

as used for computing head loss for each reach. The equation was of the 


V= Cy VRS 


hezy’s C can be written as 
c=7/88_ 
f 


erefore C,, = Gis Sg Az =C AJ? 


In Chezy’s equation the hydraulic radius R is computed by using the cross- 
‘ctional area of flow and the wetted perimeter. For these calculations the 
tal length around the flow section, including the width at the interface of the 
msity underflow and the overlying water, was taken as the wetted perimeter. 
is the slope of the energy gradient. The coefficient f is the Darcy-Weisbach 
iction coefficient for the channel reach under consideration. 

The field data from Watts Bar and Fort Loudoun Reservoirs indicated that 
| f value of about 0.016 should be used. This gives a Chezy C of 130 and for 
e of 0.0016, which is the value expected in midsummer, a C, value of 5.05. 
bsequent operations have proven these values to be reasonable. 

These backwater calculations determined the total depth of the underflow 
the confluence of the Emory and Clinch Rivers which was required before 
= diversion at the Emory River mouth could be studied. 


Vision of Flow at Mouth of Emory River 


To determine the anticipated temperature of the underflow the amount of 
Ww originating from the Clinch River underflow and the amount originating 
9m the warmer overlying waters had to be determined. These amounts 

ire found by considering the flow at the junction of the Emory River and the 
inch River to be similar to the branched flow problems in pipes. As a first 
proximation of the division of the density underflow, the discharge up the 
nory River was assumed. Thus, the interface elevation at the confluence of 
> two rivers was obtained from the modified backwater curve computations. 
was then assumed that the depth of underflow thus determined at the 
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confluence was the depth at which the underflow continued down the Clinch — 
River. It was further assumed that the flow in the lower Clinch River was ’ 
normal flow and the calculated interface depth was the normal depth. With — 
the normal depth known, the discharge in the lower Clinch River could be de- 
termined. The density underflow in the upper Clinch River was then the sum 
of the two parts. The normal depth of the density underflow was computed : 
from the equation : 


4 


3/ 2Q2 | 
SCy* Ww? 


where W is the width of the density underflow interface. Since the width-to- 
depth ratio of the underflow channel is relatively large, the wetted perimeter 
was taken as two times the width of flow at the interface. This gave the num- 
ber 2 in the above equation. 

The effect of a submerged dam located just downstream from the mouth of 
the Emory River upon the division of flow was calculated by assuming that — 
critical depth of flow would occur over the weir. The critical depth was de- 
termined from the equation 


2 
dco = oy i 


where q is the flow per unit width. : 

Figures 7 and 8 show schematized flow and temperature conditions and the 
results obtained for two computations. Figure 7 is for conditions with only 1 
the skimmer wall constructed, while Figure 8 shows conditions with the un-_ 
derwater dam added. The dam for these calculations had its top at elevati 
722, the value selected for construction. The percentage of the maximum 
condenser requirement that any particular Clinch River underflow could sup- 
ply for the two design conditions is shown in Figure 9. These curves show 
clearly the need for the underwater dam, especially when it is recognized 
that Norris turbines normally discharge an average of about 5000 cubic feet 
per second during the summer months. 

With the division of flow computed, temperature estimates were made 
based upon the curves of Figure 2. Figures 10 (a) and 10 (b) show the antici 
pated inflow temperatures for Clinch River underflows of 1500 and 5800 cubi 
feet per second, respectively, with the skimmer alone and with both skimme! 
and underwater dam. These plottings show that at a discharge of 1500 cubic 
feet per second the submerged dam would reduce the temperatures obtainabl 
from only the skimmer structure from 1 to 2 1/2 degrees Fahrenheit. Ata 
discharge of 5800 cubic feet per second the dam similarly would effect a re- 
duction of from 1 to 3 1/2 degrees Fahrenheit. Figure 10 (c) further shows 
the benefits from the submerged dam in making available at the intake mini- 
mum water temperatures at Clinch River discharges much lower than would 
be possible without the dam. With no dam it would require 9500 cubic feet 
per second of flow in Clinch River to give the same minimum temperatures 
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which, with the dam, are produced with Clinch River flows of only 4400 cubic 
feet per second—less than half of the quantity required without the dam. T. 
higher discharges are above the normal operating range for Norris Dam, 
whereas the required 4400 cubic feet per second of flow with the dam is les 
than the usual operating discharges. Thus the dam makes it possible to 
achieve minimum intake temperatures that would not otherwise be obtainable. 
The underwater dam also provides a material benefit with respect to in- 
termittent releases from Norris. Daily fluctuations of the Norris discharg 
are not noticeable at the general location of the Kingston Steam Plant, but 
shutdowns of two days which occur normally over weekends result in the 
minimum water temperature rising for several hours in midsummer to ap- 
proximately 78° F. at the elevation of the skimmer wall intake. Longer 
shutdowns result in higher temperatures for longer periods. The underwatel 
dam, however, will provide a storage pool of cold water which will require 
approximately 2.5 days to deplete. Thus, for a weekend shutdown condition 
the temperature rise would occur gradually throughout the period of draft ane 
reach the 78° F. maximum only after 2.5 days. As a-consequence the storag 
pool will materially reduce the effect of fluctuation in discharge from the { 
Norris powerhouse. 


Skimmer Wall Design 


ae Fem 


é 
The basic dimensions for the skimmer wall openings were determined by 
first computing the area of opening required in order not to exceed the allow- 
able velocity and then determining the permissible maximum height of open- 
ing. The Watts Bar-Fort Loudoun underflow data indicated that the maximum 
average velocity that could be sustained in the Emory River channel immedi- 
ately upstream from the skimmer would be 0.35 foot per second without 
breaking the stratification. These data also indicated that the velocity 
through the intake openings could be 0.70 foot per second without disturbing 
the stratification. The necessary area of opening was determined, therefore 
to be 3530 square feet. 
Physical characteristics at the wall site indicated that the bottom of the — 
opening could be at elevation 700 and that a structure about 350 feet long 
could be used. It was considered desirable to keep the height of opening as 
small as possible to minimize the opportunity for the overlying waters to be 
drawn through the skimmer. Conversely, of course, it was desirable to kee 
the length of the structure as short as possible to minimize the cost. In co- 
ordination with the design engineers, a structure consisting of five bays eae 
48 feet wide and 15 feet high with circular sheet piling piers 17 feet in diam 
ter was adopted for final study. This is the design shown in Figure 4 and is 
the design adopted for construction. The backwater curves indicated that th 
15-foot high opening would be satisfactory but that any increase in opening — 
height would result in some increase in the expected temperatures. 
Final check of the design was obtained by model studies which were spon 
sored by TVA at the Massachusetts Institute of Technology, Hydrodynamics 
Laboratory. The results of these tests are being published concurrent with 
this paper in a paper by Harleman, Gooch, and Ippen, entitled “Submerged 
Sluice Control of Stratified Flow.” Upon the basis of the results from the ~ 
M.I.T. tests, the maximum allowable velocity through the proposed opening 


would be 0.59 foot per second in the early spring and late fall when ay value 
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f 0.0005 are experienced and 1.06 feet per second in midsummer when the 
4 value may reach 0.0016. These values generally confirm the satisfactori- 


ess of the adopted skimmer wall design. 

This design provides a flow area of 3600 square feet and the maximum 
lant demand is 2310 cubic feet per second. With the allowable velocity of 
-59 foot per second from the M.I.T. tests for spring and fall, the flow 
hrough the skimmer wall opening without drawing in water from the overly- 
ng warmer water would be 2120 cubic feet per second. Since the maximum 
lant demand is 190 cubic feet per second more than this, this deficit over 
nd above that available from the underflow would have to be supplied from 
he overlying water. In the early spring and late fall seasons this is not im- 
ortant because the general over-all water temperatures during those periods 
re comparatively low and no significant loss in plant efficiency would result 
rom the small amount of slightly warmer water being drawn through the 
kimmer into the intake canal. 

The velocity of 0.70 foot per second used in the design is far below the 
alue of the allowable velocity of 1.06 feet per second found at M.I.T. for 
1idsummer periods. Hence, during the major portion of the stratification 
eason when cold water is at a premium, the structure is capable of supply- 
ng all of the demand from the coldest bottom waters. 


Performance and Benefits 


The skimmer wall and Swan Pond dikes were completed in October 1954 
nd the underwater dam was constructed during the winter of 1955-56. Per- 
ormance data are available, therefore, on the operation of the system with 
nd without the underwater dam. The efficacy of the skimmer for widely 
arying degrees of stratification is proven by the two vertical temperature 
rofiles shown in Figure 11. These profiles were measured immediately up- 
tream from the skimmer. The data in Figure 11 (a) show operation when 
= 0.0011, while Figure 11 (b) shows a profile for a oa value of 0.0004. The 
rrows indicate the temperature measured in the pump forebay at the re- 
order location marked with a triangle on Figure 1. In June the incoming 
Olar radiation heats the water flowing in the Swan Pond canal and can raise 
1e temperature by as much as one degree Fahrenheit. If this is taken into 
ccount it is apparent from the Figure 11 (a) data that only the cold bottom 
aters were passing through the skimmer wall opening. It is equally appar- 
nt from Figure 11 (b) that at very small density differences only the coldest 
ottom waters are drawn through the wall openings. 

On June 20, 1955, velocities and flow direction were also measured at a 
ection across the Emory River channel at the temperature record site, 
hown in Figure 1 just down the river from the skimmer wall. The measure- 
ents are plotted in Figure 12. The data clearly show the cold Clinch River 
aters being drawn up the Emory River. The data also show the sensitivity 
fa density underflow interface. It will be noted that the interface occurs at 
much higher elevation on the left or south bank than on the opposite bank. 
leasurements made at other stations down river showed that the channel 
irvature caused these interface variations. When it is remembered that the 
fective gravitational factor is only 0.001 of that normally encountered in 
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¢ngineering problems it is not too difficult to appreciate these interface 
inomalies. 

The real proof of the worth of these structures is contained in Figures 13 
ind 14 which show the Emory River surface, Clinch River underflow, and 
yump forebay temperatures for the stratification seasons 1955 and 1956, 
respectively. In 1955 the underwater dam was not constructed and the pump 
orebay temperatures were, as prophesied, several degrees warmer than the 
clinch River underflow temperatures. From Figures 10 (a) and 10 (b) it was 
yredicted that during August with a steady flow of 3000 cubic feet per second, 
he temperature at the intake would be about 73.59 F. Since about 19 F. of 
eating would occur in the Swan Pond canal the pump forebay temperature 
vas thus predicted to be 74.5° F. Figure 13 shows that the actual measured 
emperatures ran from 73 to 75 degrees Fahrenheit. 

Figure 14 clearly demonstrates the worth of the underwater dam. Where- 
is, the pump forebay temperatures were 3 to 5 degrees Fahrenheit higher 
han the Clinch River underflow temperatures without the dam, with the dam 
he two temperatures were nearly the same. The predicted temperatures 
shown on Figure 10 prove to be slightly higher than the actual. In late June 
[956 a steady flow of about 2500 cubic feet per second existed and resulted 
n the temperature rising to 70.59 F. The predicted temperature for this 
sondition, allowing for Swan Pond heating, was about 72° F. During August, 
vhen the Norris flows averaged about 5500 cubic feet per second, the forebay 
emperatures ran about 65.59 F. The predicted value was 69° F. 

The temperature benefits derived from the two major structures during 
1956 are graphically shown in Figure 15. 

Prior to operation, and as a part of the economic justification of the works, 
he curves of Figure 10 were used in conjunction with an anticipated Norris 
“elease schedule and plant operating cost characteristics to estimate the 
inancial benefits that would result from the diversion structures. The bene- 
its were based wholly upon the purchase cost of the additional coal (not in- 
‘luding costs for handling and processing through the plant) that would have 
o be burned if the plant produced nine-tenths rated capability while using 
-ondensing water at the temperature available without the structures. The 
kkimmer wall and the dike system, without the dam, were estimated to save 
375,000 per year. It was estimated that the underwater dam would add at 
east $24,000 per year to these savings. In actual operation, the measure- 
nents of temperature reduction show that the savings realized from the 
kimmer wall alone in 1955 amounted to $105,000. In 1956, with both skim- 
ner wall and underwater dam effective, the savings amounted to $155,000. 

f the costs for handling and processing coal through the plant were included— 
yhich it is reasonable to do—then the savings realized would be increased by 

bout 20 percent. Furthermore, if the costs were evaluated upon the basis of 
dditional energy sales, the benefits would be still greater. 

The construction costs of the skimmer wall were $251,000, of the dike 
ystem $191,000, and of the underwater dam, $31,000. These total $473,000 
or all the works. Assuming that the year 1956 is representative of the 
nnual benefits from the system, then in the short course of only three years 
r less, the total benefits will equal the entire cost. Thus, the skimmer and 
nderwater dam are very profitable additions to the steam plant facilities. 
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SUBMERGED SLUICE CONTROL OF STRATIFIED FLOW@ 
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SYNOPSIS 


Two layers of liquid differing slightly in specific gravity are contained in a 
iannel or reservoir upstream of a submerged sluice gate. The liquid down- 
ream of the gate is homogeneous and has a specific gravity the same as the 
wer, heavier layer in the upstream portion. The denser liquid is caused to 
yw through the gate opening under the action of a slight head differential. 
1e problem is to determine the maximum discharge of the lower layer 
rough the gate without steady state withdrawal from the upper layer. Experi- 


ental and analytical results are presented to define this condition of incipient 
‘awdown. 


INTRODUCTION 


Engineering interest in density currents has increased rapidly since the 
rly days of Lake Mead when observations were made of the intermittent 
ssage of silt-laden underflows through the reservoir. At the present time— 
er twenty years layer—analytical and experimental information on the dy- 
mics and behavior of density currents is still relatively meager. Density 
rrents or gravity underflows occur as a result of slight density differences 
lich may be due to temperature differentials, dissolved salts or suspended 
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sediment. The resulting stratified flow may contain one or more interfacemy 
or zones in which large density gradients are found. 

The control of density currents has long been considered economically de-. 
sirable; as an example, various barriers have been proposed for the control 
of salt water intrusion in rivers, canals and tidal estuaries. At the present 
time, the profession is beginning to realize the possibility of large scale con- 
trol of density current phenomena of even more diverse nature. The purpose | 
of this paper is to report on some experimental and analytical work dealing 6 
with gravity underflow control by means of a submerged sluice. 

The problem which resulted in this study arose during the design of the © 
condenser water intake for the Kingston steam power plant of the Tennessee 
Valley Authority. During the summer months, a thermal stratification occurs’ 
in the reservoir from which the condenser water is drawn. Since the thermal. 
efficiency of a steam plant increases as the condenser water temperature de- 
creases, there was economic justification for a control structure which would. 
result in the supply of water having the lowest possible mean temperature. A’ 
submerged sluice or skimmer wall was proposed as a means of accomplishing 
this objective. A study of the effects of density currents in the T.V.A. system 
and a discussion of the Kingston site in relation to the thermal stratification 
has been given by Fry, Churchill and Elder.(1) A companion paper being 
published concurrently by Elder and Dougherty entitled “Thermal Density 
Underflow Diversion Works for Kingston Steam Plant” discusses the Kingston: 
skimmer wall design in detail. 

Inasmuch as no attempt was made in this study to reproduce the geometric | 
or geographical conditions at the Kingston intake, the reader is referred to { 
the above papers for further details. 

The experimental and analytical studies described in this paper may be 4 
considered as a fundamental investigation of the submerged sluice gate con- 
trol of a simple two-layer system having a definable interface or density dis— 
continuity. In particular, it is desired to establish the conditions for maximt 
discharge from the lower-heavier layer without withdrawal from the lighter 
layer above. Fig. 1 is a schematic diagram showing the notation used in this 
investigation. 


Sluice Gate 
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Experimental Equipment and Procedure 


The experimental equipment consisted of a parallel glass wall flume 4.5 
inches wide, 18 inches deep and 8 feet long connected to an entrance tank or 
reservoir 18 inches square and 36 inches deep. Between the flume and reser- 
voir plastic cylindrical quadrants in both the horizontal and vertical planes 
form a smooth transition section. A schematic diagram of the flume is shown 
in Fig. 2. The adjustable sluice gate can be placed in any longitudinal position 
along the flume. Throughout the test program, fresh water was used as the 
upper layer liquid and sodium chloride solutions of varying specific gravity 
were used to simulate the colder, lower layer. The addition of vegetable dyes 
to the two liquids made tke interface between the two layers distinct. 

In preparing the equipment for an experimental run the flume and reser- 
voir were initially filled with fresh water to a depth of about six inches above 
the flume floor. The salt solution was then transferred by means of a small 
filling pump from the mixing tank to the bottom of the reservoir. By introduc- 
ing the salt solution at a very slow rate, it could be made to lift the fresh 
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water without appreciable mixing. The salt solution entered the reservoir 
through a small diffuser and baffle which served to reduce the entrance ve- 
locities and to maintain a well defined interface between the two layers. This — 
step of the process (which required from three to four hours) was continued 
until the horizontal salt-fresh water interface was a distance h, above the 
bottom of the flume. At this point, the original free surface of the fresh water 
had been raised to a distance y above the flume bottom as shown in Fig. 1. 
The sluice gate remained closed during the filling stage and salt water was 
poured into the flume downstream of the gate until it was at a level only slight~ 
ly below the free surface of the fresh water on the upstream side of the gate. 
The small difference in surface levels was necessary to avoid a reverse flow 
through the gate upon opening. After the required amount of salt water was in 
the flume, the filling pump was turned off and the larger circulating pump was — 
started. Appropriate changes in the valves of the system were made to isolate 
the mixing tank and to open the circuit connecting the suction side of the main ~ 
pump to the downstream end of the flume and the discharge side to the bottom ~ 
of the reservoir. The sluice gate was slowly raised a distance b which was 
always less than the elevation of the horizontal interface h,. The control 
valve was opened gradually and the denser liquid began to circulate smoothly 
in a clockwise direction. As long as no water from the upper layer was drawn 
through the gate the density of the circulating lower liquid remained constant. ~ 
The total discharge through the gate was obtained by means of a small Venturi 
meter in the pump discharge line. The discharge was increased slowly by 
small incremental valve adjustments. After each change several minutes 
were allowed to elapse while the interface adjusted to the new steady state 
condition. This procedure was continued until the point of incipient drawdown 
was reached. Incipient drawdown is defined as the condition at which a slight — 
increase in flow results in a depression of the interface below the gate lip 
with a corresponding steady state withdrawal of liquid from the upper layer. 
For a given geometry, density difference and interface elevation, incipient 
drawdown therefore represents the maximum rate at which liquid can be dis- 
charged from the lower layer alone. aT | 
At the time of incipient drawdown the discharge was recorded and samples j 
and temperatures of both layers were taken so that the density difference | 
could be obtained for this condition. Specific gravities were determined on a 
Westphal balance reading to four decimal places. In addition, the instantane- 4 
ous reservoir interface elevation, h,, and the profile of the interface at draw-_ 
down was obtained by point gage measurement at stations along the flume. ‘ 
During the course of a test, it was noted that the reservoir interface elevation 
increased slightly and the density of the lower fluid decreased by a small 
amount; hence, the necessity for noting these readings at the time of incipient © 
drawdown. This density change was due to the increased volume of the dense: 
fluid caused by transient discharges from the upper layer as the drawdown 
condition was approached. It was found after some experience with the flume 
that the gradual increase in the reservoir interface elevation proved to be a 
useful means of increasing the amount of data which could be obtained from a 
single filling of the equipment. In other words, for a constant gate opening, 
the condition of incipient drawdown could be obtained for several reservoir 
interface elevations during arun. Fig. 3 shows a typical plot of correspondi 
values of discharge and specific gravity differences as a function of the time 
required for a single test. At the conclusion of this run, the difference in 


specific gravities of the two layers was 78 per cent of the value at the begin- 
ning of the test. 


j : 
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Experimental Results 


A summary of the experimental runs and data for incipient drawdown con- 
ditions is given in Table I. The fourteen runs in series A were made with the 
sluice gate near the entrance to the flume (3" downstream from the point 
where the parallel walls begin) and the remaining fifteen runs in series B 
were made with the sluice gate 82 inches downstream from the entrance. Gate 
openings ranged from 1.0 to 2.5 inches and density differences between the 
two layers varied from a maximum of six-tenths of one per cent (.0060) to a 
minimum of one-quarter of one per cent (.0025). These values correspond to 
density variations caused by temperature differences between 20 and 40 de- 
grees (F) for fresh water. In addition, Table I gives the instantaneous value 
of the reservoir interface elevation and the discharge at the time of incipient 
drawdown. A dimensionless grouping of the pertinent variables (which ap- 
pears in the theoretical studies) is given in the last column. 

The following is a brief description of the observed flow conditions in the 
flume for a typical run prior to the time at which incipient drawdown took 
place. 

At the start of a run, the density interface was horizontal and coincident 
with the interface elevation in the reservoir end of the flume. After the es- 
tablishment of a small steady state flow through the gate a longitudinal line 
defining a zone of high velocity gradients immediately became evident. This 
line, which is called the separation streamline, is shown schematically in Fig. 
4, It was noted that the line of high density gradients (interface) remained 
essentially horizontal. For convenience, the wedge shaped region delineated 
by these two lines is called the separation wedge. The outline of the wedge 
was clearly observed by introducing dye below the density interface upstream 
of the gate. Any coloring injected below the separation streamline was rapid- 
ly carried away by the flow while that above remained and served to make the 
wedge visible. Due to the shear stresses at the separation streamline a low 
velocity circulation pattern was developed in the wedge region. As the flow 
rate was successively increased the curvature of the separation streamline in 
the vicinity of the gate increased and the length of the separation wedge gradu- 
ally decreased due to entrainment with the moving lower layer. The disap- 
pearance of the wedge was also hastened by transients accompanying the suc- 
cessive changes in discharge. At the time of incipient drawdown, the density 
interface coincided with the separation streamline. Fig. 5 shows a series of 
interface profiles at the drawdown condition for runs 63, 65-68. In this series, 
the density difference and gate opening were constant while the reservoir 
interface elevation increased. At flow rates beyond the point of incipient _ 
drawdown, the density interface passes through the gate below the lip and is 
maintained for some distance downstream. The two-phase flows were not 
investigated quantitatively due to the difficulty of maintaining steady state 
conditions with the existing equipment. 


Theoretical Considerations 


A completely general analytical approach to the submerged sluice gate dis- 
charge of a stratified fluid is complicated by the presence of the separation 
wedge previously noted. Hence, the theoretical considerations will be con- 
cerned with flow conditions at incipient drawdown when the separation wedge 
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TABLE I 
SUMMARY OF RUNS AT INCIPIENT DRAWDOWN 
Run No. Gate Interface hy hy Discharge at Fi = 
Opening Elev. in 5a = Incipient b 
b(in.) Reservoir Y, Drawdown Y 
h, (ina) Q, (in? /sec) g7b 
A, Sluice Gate at Flume Entrance (St 0 
ami 1.0 1.60 1.60 .0050 5.6 0.90 
12 1,0 2.25 2625 0049 9.0. > 1.51 
ai 1.5 1.90 1.27 60062 6.8 0.53 
22 1.5 2.35 1.57% 0060 10.3 0.82 
23 NEE) 2.60 1.73 0059 12.6 1,01 
24, 1.5 3.30 2.20 0050 17.4 1.52 
25 1.5 3.95 2.64 0046 23.2 Zell 
31 2.0 240 1.20 8.0047 Tee 0.42 
32 2.0 2.90 1.45 0047 11.4 0.68 
33 2.0 3.20 1.60 .0046 14.2 0.84 
34 2.0 3.50 1.75 0043 18.7 1.15 
Al PRLS) 2.40 1.36 0025 11.4 0.65 
42 205 3.80 L529 00R5 15.4 0.88 
43 en 4290 1.96 .0025 22,0 1,26 


B. Sluice Gate at Sta. 82 


Du 1.5 2205 1.37  .0040 3.0 0.29 
ae 1.5 2.20 1.47 0040 4.0 0.39 
5 1.5 2.40 1.60 ,.0040 6.0 0.59 
5A 1.5 2.50 1.73 0040 Tel 0.69 
DD 1.5 2.90 1.93 0040 9.0 0.88 
56 1.5 3.15 2.10 .0040 11.0 1,07 
57 1.5 3245 2.30 .0040 12.9 1.26 
61 2.0 2.60 1.30 0056 4.0 0.21 
62 2.0 2.70 1.35 0056 6.0 0.33 
63 2.0 2.80 1.40 .0053 8.0 O44 
64, 2.0 2.95 1.48 0053 9.0 0.49 - 
65 2.0 3,10 1.55 0053 11,0 0.55 
66 2.0 3035 1.68  .0053 12.0 0.66 
67 2.0 3BV65 1.82  .0053 14.0 0.77 
68 2.0 3.85 1.92 .0053 17,0 0.93 
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has disappeared and the interface has the form shown in Fig. 1. In the analo- 
gous single-fluid submerged-sluice problem, the normal analytical approach 
is to write the total head equation between a section of parallel flow upstream 
of the gate and the section through the vena-contracta below the gate. There- 
fore, the surface elevation difference across the sluice is introduced as a 
primary variable in the discharge equation. This procedure is not advan- 
tageous in the two layer system under consideration, since (as will be shown 
later) the presence of the upper liquid greatly reduces the gravitational driv- 
ing force required to produce flow and results in surface elevation differences 
across the gate which are too small to measure with any degree of accuracy. 
Thus, it is desirable to avoid introducing this term which, aside from the 
measurement difficulties, is not particularly useful in practice. The down- 
stream section for evaluation of the total head is therefore chosen to coincide 
with the plane of the sluice gate. It is recognized that in choosing this section 
the curvilinear nature of the flow will have to be considered in the total head 
equation. 


i 


Referring to the definition diagram, Fig. 1, the total head equation is written © 


between the reservoir and the vertical section at the gate using the flume 
bottom as a datum. It is assumed that the stratification originates in an area 
of such extent that the flow velocities in the reservoir region may be neglect- 
ed by comparison with those in the vicinity of the gate. If AY represents the 
difference in specific weights of the two fluids, then Y and Y - AY are the spe- 
cific weights of the lower and upper layers respectively. The total depth of 
the two layers upstream of the gate is y and h, equals the elevation of the 
density interface in the reservoir. The mean velocity Vp at the gate is given 
by the discharge per unit width dividided by the gate opening b. Using the 
above notation, the total head equation becomes: 


2 
eS as h, = ealyce) ae (b- Ahm) + Sal + h (1) 


The left side of the equation represents the piezometric head and also the total 
head in the reservoir in accordance with the assumption of negligible reser- 
voir velocities. The first two quantities on the right side represent the mean 
piezometric head in the plane of the gate. Due to the curvature and inclination 
of the streamlines in this region the piezometric head varies between the gate 
lip and the flume bottom. The length Ahp is therefore introduced to express 
the amount by which the mean piezometric head differs from the value b which 
would be obtained for hydrostatic pressure distribution. The kinetic energy 
head at the gate is represented by the third quantity on the right side of Eq. (1) 
where o is the kinetic energy correction factor required by the non-uniform 
‘velocity at the gate section and the subscript i refers to the critical condition 
of incipient drawdown. Finally, hg is the total resistance loss in the approach 
channel between the reservoir and the sluice gate. The total loss is a function 
of the boundary shear on the bottom and side walls of the channel and the 
frictional resistance along the interface. 

Expanding the terms containing (7 - AY), Eq. (1) reduces to: 


2 
Abn, = use Ahm + oo + hy (1a) 


which can be put in the following dimensionless form after dividing each term 


an 
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hy j Atm + vy + h 
== = : 
b a 25 oo (2) 


It should be noted that Eq. (2) represents the interrelationship of the basic — 
variables h,, b, A?/Y and Vp only at the critical point of incipient drawdown. ~ 
Prior to the drawdown condition, the total head equation should contain a term 
accounting for the additional pressure at the gate due to the separation wedge. 
Beyond the incipient point the equation should again be changed to account for 
the two-phase flow through the sluice gate. Eq. (2) is therefore useful in 
correlating the experimental data which was obtained at precisely the flow 
condition for which the equation is applicable. The term containing Vp is 
readily seen to be a generalized Froude parameter in which the acceleration 
of gravity is reduced by AY/Y due to the bouyant effect of the upper layer of 
fluid. Defining the reduced gravitational constant as g' = g(A?/Y), then 


; 


V;): 
iv — (3) 

g'b 
and Eq. (2) becomes: 3 ! 
hy Ah oC Fj h¢ ; 
—=j-Sm+ St +a f 
b at, 2 ab, (4) 


b 
Eq. (4) can be shown to reduce to the elementary concept of free surface flow 
on a broad crested weir in which the critical depth (b) is assumed to exist 
over the entire length. For the case in which the upper fluid is air AY ~7, 
neglecting curvature Ah,, = 0, for uniform velocity distribution @ = 1 and neg- 
lecting friction hg = 0; then hy/b = 3/2 since Fj = 1. 

In order to simplify the correlation of the experimental data with the above 
analysis, emphasis will be placed primarily upon the situation in which the ~ 
sluice gate is located near the reservoir entrance. This is equivalent to the — 
assumption that the head loss term in Eq. (4) is negligible. Approximately 
one half of the experimental runs (section A of Table I) were conducted with 
the gate three inches from the flume entrance and may be assumed to satisfy — 
this condition. 

Various aspects of the problem of pressure distribution in curvilinear fre¢ 
surface flow have been treated by Rouse, (2) Jaeger(3) and more recently by 
Khafagi and Hammad. (4) The latter have developed expressions for pressure 
distributions appropriate to the type of curvilinear flow obtained near the 
critical sections of free surface Venturi flumes. The influence surface incli- 
nation (streamline convergence) and surface curvature were analyzed inde- 
pendently with the combined mean departure from hydrostatic conditions bei 
given by an equation of the form: 


Vp 
Ahm = 5e(f(R) + £,(2)) (5 


where R is the radius of curvature and 6 the angle of inclination of the inter- 
face (or surface) as shown in Fig. 6. The Khafagi-Hammad functions for R 
and @ are given by: 


i 


= i 
a (Tae 
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Sy ea = YL. eee 
f,(@) = | rE where @G= 6 (In tan( > + )) (7) 
bstituting Eq. (5) into the total head Eq. (4) and noting that hr = 0, the follow- 
tis obtained if R and @ are evaluated at the drawdown condition: 


2 2 

hy Fj oF: 

a i= S (4 (R) + ,(@)) + oe (8) 
ice it is known that for given boundary conditions the surface configuration 
free surface flow is a function of the Froude parameter it may be assumed 
it the geometric functions given by Eqs. (6) and (7) are therefore dependent 
on Fj. For convenience, a new parameter 7 is defined such that, 


ai F(t (R) + £,(@)) 
OS 2M 2 (9) 


len R =co and 6 = 0, the pressure distribution must be hydrostatic and it 

n readily be shown that for this condition 7 reaches a maximum value of 
ity. According to the foregoing reasoning, the dependence of 7 on F; can be 
restigated if R and @ can be obtained experimentally. In two of the experi- 
tal runs at the drawdown condition photographs of the interface at the gate 
re obtained which were sufficiently clear for evaluation of R and @. The 
9erimental data and the plotted values of 7 and Fj are shown in Fig. 7. 

me indication of generality and an extension of the relation of higher values 
F; can be obtained by making use of some experimental results of Khafagi 
i Hammad for free surface flows having combined inclination and curvature 
a Venturi flume. These values are also shown in Fig. 7 and compare favor- 
ly with the stratified flow results. 

Combining Eqs. (8) and (9) the basic total head equation is given by: 


2 
h AF; 
rake VK + = (10) 


ere the dependence of 7 on Fj is given by Fig. 7. 

Eq. (10) can be plotted using hy/b and F; as independent variables and a 
nily of curves depending on the kinetic energy correction @ will be obtained 
shown in Fig. 8. The corresponding experimental data from part A of 

ble I are also shown for comparison with the analytical results. In’general, 
experimental data follow the trend indicated by Eq. (10). For the low 

ues of Fj the flows are laminar and values of o between 1.75 and 2.0 are 
unusual for a rectangular channel having a width to depth ratio of approxi- 
tely 1 to 2. For increasing flow rates a comparison of the experimental 

a with the curves of Eq. (10) indicates a progressive decrease ino toa 
nimum value of approximately 1.10. The decrease in the kinetic energy 
‘rection factor is attributed to increasing curvature and acceleration rates 
the onset of turbulent flow. No attempt was made to measure o due to the 
ee-dimensional character of the flow and the low velocities involved. 

The experimental data for the sluice gate 82 inches from the flume 

rance (part B of Table I) are also plotted in Fig. 8 for comparison with the 
ve. In accordance with Eq. (4) the theoretical curves for this condition 
tld be displaced upward due to the friction loss in the flume. For identical 
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values of F, the vertical displacement between points obtained at the flume — 


entrance and at station 82 should equal the dimensionless term be . No pre: 
AY ; 


¥ : 
cise methods for determining the energy loss in non-uniform subsurface flow 
have as yet been developed. Only the resistance problem for laminar uniform 
underflows has been treated in any detail. 5 


Illustrative Example 


As an illustration of the use of Fig. 8 in determining the drawdown con- — 
dition, the computations for the Kingston skimmer wall will be presented. 
TVA had determined that in early spring and late fall AY/Y values as low as 
0.0005 exist while in mid-summer the density differences may produce value 
of 0.0016. The proposed TVA design was for a skimmer wall opening b = 15 
feet. Their measurements and calculations indicated the elevation of the inte 
face upstream of the skimmer wall (hy) to be about 24 to 25 feet (zone of max 
mum rate of change of density with depth). For the relatively high Reynolds 
numbers of the prototype flow the kinetic energy factor might reasonably be 
assumed at ao = 1.10. From Fig. 8, for hy/b = 1.63, the corresponding value | 
of F, = 1.20. Hence, in the spring and fall with gt = 0.0005 (32.2) = 0.016 ft./' 
sec.’, the critical drawdown velocity (Vp); = 0.59 ft./sec. is obtained directly: 
from Eq. (3). For the mid-summer condition when g' = 0.0016 (32.2) = 0.052 | 
ft./sec.” the critical drawdown condition will occur at a velocity of 1.06 ft./ _ 
sec. The maximum discharge per foot of width from the colder underlying 
layers would be therefore qj = (Vp)j (b) = 8.8 cfs/ft. and 15.9 cfs/ft., re- 
spectively. If a discharge in excess of the critical values given above is ob- 
tained, steady state withdrawal of the warmer overlying fluid should begin. 

Summary and Conclusions | 

The experimental and analytical work described herein indicates that the! 
control of stratified flows involving very small density differences is both 
possible and practical. The primary aim in the presentation of the experi-— 
mental data for the small scale laboratory work was to achieve a correlation 
which would permit application of the information to the design of engineering 
structures. The family of curves in Fig. 8 delineate a zone of incipient draw 
down which depends upon the kinetic energy correction factor. The region t 
the right of the curves represents a combination of the parameters hy,/b and 
Fj which result in simultaneous draft of fluid from both the upper and lower 
layers. In the region to the left of the curves, flow from the lower layer onl 
will occur at less than the maximum rate. 

In engineering applications it can reasonably be assumed that high values 
of @ will never be obtained. The uncertainty in the choice of the proper ki 
energy factor, however, is by no means as important as that involved in the 
assumption of a discrete density discontinuity for the determination of hy. 1} 
practice the assumed interface might be logically chosen to coincide with the 
elevation at which the rate of change of density with depth is a maximum. ~ 
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MECHANICAL ANALOGS AID GRAPHICAL FLOOD ROUTING* 


Max A. Kohler,** M. ASCE 
(Proc. Paper 1585) 


The dampening effect of valley storage on a wave moving through a river 
ach is a familiar phenomenon, and hydrologic literature abounds with 
scriptions of routing methods—analytical, graphical, electric analogs, and 
chanical devices.(1,2,3,4,5) A similar dampening is apparent if one com- 
res the rear-wheel track of a moving vehicle with that made by a front 

eel, although closer examination will reveal that the two effects are not 

ke in all respects for a vehicle with fixed wheel base. Nonetheless, study 
the two phenomena reveals some rather interesting results, pointing the 

y for development of mechanical analogs and simplified graphical solutions 
flood routing problems. 


Linear Storage-Discharge Function 


Let us first consider the case of a linear storage-discharge relation, i.e. 


S =KQ a) 


ere S is the reach storage, Q is the reach outflow, and K is a constant hav- 
‘the units of time. Figure 1 shows schematically the requirements for an 
log which will satisfy Equation 1. In this diagram, the time scale for out- 
w is shifted K time units to the left of that for the inflow, and the tongue 
necting the inflow and outflow pens is adjustable so as to maintain the time 
placement between the two at K units. This is the underlying principle of 
patented Harkness flood router.(4) It will be seen that the slope of the out- 
w hydrograph dQ/dt is equal to the difference between inflow and outflow 
ided by the storage factor K. Since I - Q is the time rate of change in 

rage 


8: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1585 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 

Paper presented at the Annual Convention of ASCE, New York City, Octo- 
ber, 1957. ; 

Chief Research Hydrologist, U. S. Weather Bureau, Washington, D. C. 
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dQ (2) 
dt 
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Qa. 
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. 
K 
which is equivalent to Equation 1. 
A simple, graphical solution for Equation 1 is also apparent from Figure ? 

If the outflow up to any initial time is plotted K units to the left of the inflow 
hydrograph to be routed, a straightedge connecting simultaneous values of 
and Q has the slope dQ/at, and the outflow can be extended for a reasonable 
time with this slope. (3) Moving to successive later times, the inflow hydro- | 
graph is routed segment by segment. The routing period (i.e., the length j 
with which the outflow curve is extended in one step) can be varied at random 
depending upon the rate of change of dQ/dt. Moreover, the outflow and inflow 
time scales can be made to coincide by pla placing the straightedge so as toc 
nect the outflow with simultaneous inflow lagged K time units (see Figure 

In this case, it is preferable to have a straightedge calibrated in time units 
to assist in lagging the inflow. 


Muskingum Storage Function 


The Muskingum method of routing developed by McCarthy(1) is based 
the assumption that storage is proportional to a weighted value of inflow 
outflow, i.e. 


S=Kilxl tex ror 


where x is visualized as a weighting factor which takes into account wedge 
storage. Figure 3 illustrates schematically the requirements for a mec 
cai analog satisfying the Muskingum storage function. In this case the to: 
is jointed at some point, the location of which is determined by the relative 
weights assigned to I and Q in the storage function. In addition to tracing tl 
inflow hydrograph, the inflow portion of the tongue must always be pivoted 

as to have the slope dI/dt. Then 


which is equivalent to the Muskingum storage function. 

Graphical Muskingum routing can also be accomplished with the aid of 
jointed straightedge of the general character shown in Figure 4. Actually, 
experience has shown that the Muskingum equation provides no more flexi 
bility than does the linear storage-outflow relation employed in connection 
with a lagged inflow. Since graphical solution of the simpler equation can 
readily incorporate fixed lag L, there is little to recommend use of the 
Muskingum function. 


Variable Lag and Storage Factor 


There are numerous reaches for which neither the Muskingum nor the 
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ag-and-route (constant K and L) techniques provide satisfactory results. In 
act, it is quite generally - found that the results of these approaches can be 
mproved upon if K and L are permitted to vary with flow in the reach. It 

an be shown that inflow and outlfow volumes do not check when K is assumed 
0 be a function of inflow and, therefore, K must be related to outflow. Lag 
an be made a function of either inflow or outflow. That is, the inflow hydro- 
raph can be lagged as a function of the inflow discharge and the result routed 
vith a variable K, or the inflow can be routed with K as a function of outflow 
nd the resulting hydrograph lagged. While reasonably comparable results 
an be obtained in either manner, practical considerations, both in develop- 
nent and application, favor the assumption that lag is a function of inflow. 

It is not the purpose of this paper to describe in detail the actual develop- 
nent of routing procedures, but a brief outline for the variable K and L tech- 
ique will aid in following the subsequent discussion. If storage is a function 
f outflow only, then it follows that the peaks and troughs of the routed outflow 
i “ = OandI=Q). Therefore, 
, good approximation to the lag function can be derived by analysis of a series 
f inflow and outflow hydrographs (adjusted for local inflow) as shown in 
figure 5. The time difference between the outflow peak and the correspond- 
ng value of inflow is the lag to be used for that inflow. The mean line 
hrough a plotting of such points constitutes a first approximation of the re- 
uired lag function. 

Replotting the inflow hydrographs lagged in accordance with the relation 
erived as above, the values of K can be obtained as shown in Figure 6. Plac- 
ng a straightedge tangent to the outflow curve (at any point except where 
2 —> 0), Kis read as the time difference between the lagged inflow (at 
ime of selected outflow) and the straightedge (Equation 2). Having determined 
; at each of a number of points for selected rises, the relation K = f(Q) can be 
erived by plotting K vs Q. 

Relations for K and L developed in the manner described are usually high- 
y satisfactory. Nevertheless, comparison of routed and observed hydro- 
raphs may indicate improved results through minor adjustments. The analy- 
is for determination of K and L may show there is little to be gained by per- 
litting these factors to vary with flow, and in this case the simpler solution 
3 to be preferred. Limited experience has shown that K = f(Q) and L = f(1) 
re quite similar for a particular reach—if L or K varies appreciably with 
low, then both vary in a similar fashion. 

Graphical routing with variable K and L, although more complex than when 
sing constant factors, is more rapid than many less flexible methods. It is 
reatly facilitated by the aid of a transparent template of the type shown in 
‘igure 7. The left-hand edge of the cutout section represents L = f{(I), and 
ie right-hand edge represents K = {(Q), both as measured from the vertical 
ase line. The template is moved horizontally along a straightedge to deter- 
line a segment of the outflow hydrograph centered at the base line. Use of 
1e template is illustrated in Figure 8. Intersection of the inflow hydrograph 
ith the left-hand edge of the template represents lagged inflow at time indi- 
ated by the base line. Eye extrapolation of the outflow curve to the base line 
etermines the value of K on the right-hand edge of the template. The point 
A” thus obtained is connected by straightedge to the previously-determined 
egment of the outflow hydrograph and a reasonable extension to the curve is 


all on the inflow hydrograph (i.e., when so oO 


Sa fa 
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constructed. This process is repeated until the routing operationis com- __ 
pleted. 
The mechanical analog for accomplishing this form of routing must 
naturally be more complex than that for the simple, linear function. The 
time difference between the inflow tracer and the outflow pen must be — 
in accordance with L = f(I), while the time difference between the outflow pen 
and the anchor point for the outflow tongue must follow K = f(Q). Modifica- 4 
tion for variable lag is reasonably simple, but variations in K must be ge 
to the outflow pen, and this presents quite a different problem. The required | 
mechanism can be readily designed, but vertical movement of the outflow pen | 
would be forced to “work” the mechanism, which would tend to induce side — 
thrust on the wheels and thus possibly distort the outflow hydrograph. .. 
Although a more difficult task, it was believed an operator could soon ac= 
quire skill to follow the inflow hydrograph and simultaneously make mechani-: 
cal adjustment for variations in K. Figure 9 is a photograph of an analog ~ 
which has been designed and constructed accordingly. The left-hand knob is 
used for maintaining the crane assembly at the intersection of the inflow and 
lag curves while the carriage is moved to the right along the time axis. The 
right-hand know is used to vary K in accordance with the routed outflow. 


? 


Preliminary operational tests indicate this pilot model to be quite satisfac- — 

tory, but it does have certain limitations which will require minor modifica= 

tions for improved efficiency and convenience. ry 
*) 


SUMMARY AND CONCLUSIONS 


Graphical techniques for routing directly on the plotted hydrograph chartil 
have been discussed for (1) the case of simple, linear storage-discharge 
function, (2) the case of the Muskingum equation, (3) the case of linear stor- — 
age function in conjunction with constant lag (translation), and (4) the case 
both lag and the storage factor varying as complex functions of flow in the 
reach. The fourth category represents what is believed to be one of the mos 
flexible methods of routing yet devised. Amazingly reliable results have 
obtained for reaches which are incompatible with storage functions of the 
Muskingum and other recognized methods of routing. The time required for 
procedure development compares favorably with any other known technique. 
Moreover, the graphical solutions illustrated are operationally rapid and re 
liable. The basic features of mechanical analogs comparable to the graphii 
solutions described have also been discussed. An analog which will facilitat 
variable K and lag has now been constructed and is undergoing tests. 
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TIDAL MOVEMENT IN THE CAPE COD CANAL, MASSACHUSETTS2 


B. W. Wilcox! 
(Proc. Paper 1586) 


SYNOPSIS 


Variations in the form of the tide wave at selected points, analyses of ob- 
srvations to obtain harmonic constants, and a practical method for predict- 
ig the tidal currents in the canal are discussed. Diagrams are included to 
iow the shapes of the tide curves and also the agreement between predic- 
ons and observations of the tidal current at the railroad bridge. 


INTRODUCTION 


For many years the U. S. Corps of Engineers has been engaged in the im- 
“ovement of the Cape Cod Canal. The land cut has been widened and deep- 
led. Dikes have been built and new channels dredged in Buzzards Bay. 
ig. 1 shows the canal and the greater part of Buzzards Bay. The improve- 
ents have changed the tide and the tidal current in the canal necessitating 
~w observational data to establish the pattern for predictions which are of 
msiderable importance to the navigator. 

This discussion of tidal movement divides naturally into two parts: the 
rtical rise and fall of the water which is ordinarily called the tide; and the 
rizontal flow or tidal current. 

Automatic tide gages furnished by the Coast and Geodetic Survey have 
en operated by the canal engineers at a number of selected places during 
e years 1955-1957. During this same period tidal current observations at 
th ends of the land cut were made by the Coast and Geodetic Survey. 


te: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1586 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 
Presented at meeting of the ASCE, New York, N. Y., October, 1957. 
Chief, Predictions Branch, Tides and Currents Div., Coast & Geodetic 
Survey, U. S. Dept. of Commerce, Washington, D. C. 
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Tides 


The difference in range and phase of the tide in Buzzards Bay and Cape 
Cod Bay produces a periodic change in the slope of the water in the canal 
which in turn causes the periodic movement of water in the canal. Slope 
lines were drawn for 12 lunar hours between stations 625,379,200 and 35. 
These slope lines determined the high and low water profiles of figure 2. 
Observations at station 452 and 320 fit well into the profiles thus furnishing © 
a check on their validity. The high and low water profiles in figure 2 show 
the slight decrease in mean range from 3.9 feet at Abiels Ledge, Buzzards — 
Bay to 3.5 feet at the entrance to the land cut (RR. Bridge) and then the sharp : 
rather uniform increase to 8.7 feet at the mooring basin (Sandwich) near the 
Cape Cod Bay entrance to the canal. 7 

The great difference between the tides at the two ends of the canal is 
brought out clearly in figure 3. The tide at Sandwich (station 35) is ofthe 
regular semidaily oceanic type with the duration of rise equal to the duration © 
of fall. The distortion appearing in the Buzzards Bay wave is due to what is 
commonly called “shallow water effects.” The duration of rise is not the 
same as the duration of fall. Note the change in the low water part of the ~ 
tide wave as it progresses from station 625 to station 379. At Abiels Ledge 
(station 625) the rise of the tide takes longer than the fall whereas at the : 
railroad bridge (station 379) the duration of rise is less than the duration of 
fall. Before the change in the channel and the building of the dikes the shape 
of the curve for station 379 was similar to that for station 625. 

The tide changes continually as the tide producing forces vary with the 
motions of the moon, sun and earth. Inasmuch as the equations of motion of 
moon, sun and earth are known it is possible to compute the tide producing 
forces for any desired time and hence the resulting tide. Present theory is 
not sufficient for mathematical prediction of the tide without observational 
data. Observations are necessary because each body of water reacts to the — 
tide producing forces in its own way. i 
The well known harmonic equation(1) of the periodic tide is ; 

» 


y = Zo +L fH cos (at + Vg + - K) 


in which 


y = height of tide at time t 

Zo = height of mean sea level above datum of predictions 

H = mean amplitude of any constituent A 

f = factor for reducing mean amplitude H to year of prediction 
a = speed in degrees per solar hour of constituent A 

t = time in solar hours reckoned from some initial epoch 


Vo+u = theoretical phase (equilibrium argument) of constituent A 
when t =0 


K = difference between theoretical and actual phase of constituent A 
as determined from observations 


The number of constituents involved will vary with the complexity of the 
tide. 
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The mean amplitude (H) and the epoch (x) for each constituent are the 
armonic constants to be obtained from observed hourly heights of the tide. 
he problem is simplified by the fact that the speeds of the constituents 
hich may be involved are known in advance. There are a number of ways 
) make the analysis. The one used by the Coast Survey involves the assign- 
rent of each observed hourly height to a constituent hour. There are 24 con- 
tituent hours in a constituent day. Twenty four sums (one for each constitu- 
nt hour) are obtained from the entire series. Modern electronic computing 
1achines are now used in the summation of long series. The 24 hourly sums 
xr constituent A contain the augmented amplitude of A but the constituents 
f other periods are largely cancelled. Short period meteorological effects 
hich are contained in the observations are filtered out in the process of 
nalysis. From the 24 hourly sums the amplitude H and the epoch x of A 
re readily obtained. Some refinements are required to eliminate inter- 
srence of constituents which the process of summing can not complete. 

The harmonic constants thus obtained may be used in making predictions. 
his is an arduous task unless special machines are used. Many of the 
rincipal maritime nations have had analog computers built especially for 
dal predictions. By the process briefly outlined above the Coast and 
eodetic Survey has analyzed the tidal observations made in 1955 at stations 
79 and 35. The principal tidal harmonic constants are listed in table 1. 


Table 1 
Tidal Harmonic Constants(1) 


M2 S2 N2 Ki O1 Py M4 M6 


fa. 35 H ft. 4.08 0.62 0.91 0.44 0.35 0.14 0.12 0.14 
K deg. 326 002 286 135 111 135 106 243 


ta. 379 H ft. 1.62 0.28 0.46 0.27 0.18 0.09 0.29 0.07 
K deg. 268 276 244 113 128 113 139 236 


he tidal harmonic constants for these two stations adjusted by current ob- 
srvations are used in predicting the tidal current at station 379 - the method 
ill be discussed later under currents. 


Currents 


Strong tidal currents flow back and forth periodically through the land cut, 
ie westbound current being the stronger. A velocity of 5.3 knots has been 
served in midchannel at the railroad bridge. In entering Cape Cod Canal 
nd cut from Buzzards Bay it is natural to call the current from Buzzards 
ay into the land cut the flood current; but the canal may also be entered 
om Cape Cod Bay and it is therefore just as natural to regard the west- 
ard current as the flood current. In the canal, clearly, upstream and down- 
ream cannot be used in designating the flood and ebb current, and so the 
me relations between current and tide must be used. The observations 
1ow that the strength of current setting eastward at the RR. bridge occurs 
yout three hours before local high water; that is, on a rising tide. Hence 
e flood current at the RR. bridge is the one setting eastward from 
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Buzzards Bay into the land cut. In fact the eastward current is the flood cur-. 
rent throughout the canal for the tide and current data at the Cape Cod Bay © 
entrance also show that the eastward current is on the rising tide. 

The following table gives the average time intervals referred to the 
Greenwich transit of the moon for various phases of the tidal current and 
the average maximum velocities for the entrances to the land cut as deter- 
mined from observations in 1955 and 1956. 


Table 2 
Slack Maximum Slack Maximum Velocities 
Flood Ebb Flood Ebb 

hours hours hours hours knots knots | 
RR. Bridge 8.25 11.06 reeves 4.68 4.0 4.5 
Sta. 379 
Mooring Basin 8.03 10.86 2.03 4.45 2.4 yD | 
Sta. 35 


Observations at station 379 were made with Price meter at depths 10 feet 
and 20 feet below the surface. The currents at the two depths were so nearly 
equal that the values were meaned for table 2. At Station 35 the currents 
were measured at 15 foot depth by a Price meter. The time of the current 
is practically the same throughout the canal which verifies the theory that 
the flow is almost entirely governed by the changing surface head in the 
canal. The velocity at the RR. bridge is much greater than at the Mooring 
Basin. This is due to the tide being about two hours later at the Mooring 
Basin than at the bridge. 

Advance information about the currents in the canal is extremely important 
to the navigator. Daily predictions of the slacks and strengths of current con- 
tained in the annual tidal current tables published by the Coast Survey are 
based on observations of both tides and currents. The differences in tidal 
range and the times of high and low water at the two ends of the land cut 
cause the water surface at one end to alternately rise above and fall below 
the level at the other end, thus creating a periodic reversing current inthe 
canal. Theoretically, disregarding friction or inertia, the velocity of the 
current would vary as the square root of the surface head, being zero when — 
the surface is at the same level at both ends and reaching a maximum when | 
the head is greatest. Actually there will generally be a lag of some minutes © 
in the response of the current movement to the surface head which must be © 
determined by observations. The problem of predictions therefore, reduces — 
to one of predicting the surface head between the two ends. This may be done 
by means of harmonic constants derived from the tidal constants for the 
entrances. 

For any one constituent let T represent time as expressed in degrees 
reckoned from the phase zero of its Greenwich equilibrium argument. Also 
let Yg and Yb represent the height of the constituent tide for any time T < 
referred to the mean level at stations 379 and 35 respectively and let Yg = 
surface head = Yy-Yp. Tidal harmonic constants for stations 379 and 35 ar 
obtained from analyses of observations as outlined under tides. H, is the 
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mean amplitude and G, is the Greenwich epoch of a constituent for station 
379; Hp, and Gp are the corresponding harmonic constants for the same 


constituent at station 35. Formulas for the heights and surface head may 
now be written 


Yq = Hg cos (T-Gg) for station 379 
Yb = Hp cos (T-Gp) for station 35 
Yq = Ha cos (T-Gg) -Hp cos (T-Gp) 


Harmonic constants (Hg and Gg) for the surface head may be obtained by tak- 
ing the vector difference between H,,G, and Hp,Gp. Hence the surface head 
for any one constituent may be written Yy = Hg cos (T-Gg). 

Separate computations are made for each of the tidal constituents. The 
values obtained for Hg and Gg are the corresponding amplitudes and 
Greenwich epochs in a harmonic expression for the continually changing dif- 
ference in elevation of the water surface at the two ends of the land cut. 

For predictions of the current the amplitude Hg must be modified 
through the formula (Velocity)2 = constant (C) x Hg. Furthermore the epochs 
Gq need to be changed to allow for the lag in response of the current to the 
changing difference in water level at the two entrances. 

Following the above outlined method the tidal harmonic constants from 
table 1 were combined vectorially and adjusted for the observed current at 
the RR. bridge. The constants are listed in table 3. 


Table 3 


Tidal Current Harmonic Constants 
M2 S2 N2 Ky O; Py Mg Mg 


RR. Bridge Hft. 8.73 1.65 1.62 0.52 0.47 0.17 0.51 0.18 
Sta. 379 K deg. 185 223 150 353 281 352 189 117 


With the modified harmonic constants direct readings from the tide 
predicting machine would give velocity squared. However a special square 
root dial on the machine makes it possible to read the velocities directly. 

Figure 4 shows a comparison of predicted tidal current with observed. 
The period was picked at random. The differences are small with one ex- 
ception. Part of the deviation of predictions from observations is due to non i 
periodic variations in the observations and part is due to the omission of 
factors other than surface head involved in the current. But it is evident 
from Figure 4 that the other factors have a very small effect and that the 
predictions through surface head are dependable. 


REFERENCE 


1. Paul Schureman, Manual of Harmonic Analysis and Prediction of Tides, 
Revised (1940) Edition, U. S. Coast and Geodetic Survey Special Publica- 


tion No. 98. 


HY 2 


1586-6 


wimoL 190 10010 


‘= ag SFL 
} : r- =~ AY 
— Cie Cra 
> = he 
te = St ys 
Says: Bae 
ws. : 
ey oO > 
ae 3 


‘SSVW “IVNVD GOD Adv9d 


CAPE COD CANAL 1586-7 


ASCE 
° 
BE yoimpues 
° 
ro) 
N 
a 
<= 
— 
<< 
oO a 3 Pie 
a st aspiug H 9uinog 
oO 
© 
WW o 
fal 6 espug “YY 
<= 
oO 
® pues} 30H 


x) oO 
= = 
° 2 
oa a 
5) 
lace ~ 
oO oO 
= = 
s 

& é 
= =H 


28pe7 sjaiqy 


Station 625 


April, se 


HY 2 


1586-8 


cl 


0 


T 8 9 
SSAUNND ACIL NVA 


Ja}eEM MO] — NM] 
49}eM YSIY- MH 


eveccccccese CE UOI}E}S 
ome 6/£ U0!}21S 


lee ime CO NUOIETS 


0 sunoH 


1586-9 


CAPE COD CANAL 


ASCE 


cl 


"SOI}IDOJAA PaAJOSGO ase SO ‘peay adejuns Ysnoiy} paydipesd sem aAund pijos 
qqg wnwixey| 
O 


OT 9 v é 


8 
6ZE NOILVLS SLNAYYENO WIL 


N 
es ee ey 


s}ouy 
Ay00/8A 


O sunoy 


Paper 1587 April, 1958 HY 2 


SSS mnmmmmememnmmneneeeeeeeeeeeeeee ee ST 
Journal of the 


HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 
a 


DISTRIBUTION OF SEDIMENT IN LARGE RESERVOIRS2 


Whitney M. Borland,* A.M. ASCE 
and Carl R. Miller,** A.M. ASCE 
(Proc. Paper 1587) 


SYNOPSIS 


In planning and developing multi-purpose reservoirs it is necessary to 
make an allocation of storage space for sediment accumulation. This sedi- 
ment accumulation, in turn, influences the design of the dam because of the 
effect on active storage capacity requirements, outlet sill elevations, recrea- 
tional facilities, and backwater conditions. It is, therefore, necessary to 
predict the probable distribution of the sediment in the reservoir with time. 

Two methods of predicting the probable sediment distribution are given 
and an illustrated example shown. One method is basically mathematical 
(Area-Increment Method), whereas the other is a mathematical procedure 
developed on the basis of actual experiences (Empirical Area-Reduction 
Method). The area in which additional study is needed to evolve procedures 
that will give consistently reliable predictions of sediment distribution under 
varying conditions is pointed out. 


INTRODUCTION 


Unfortunately, when a dam is constructed to store water for useful pur- 
poses, sediment being transported by the water is also stored. Unlike water, 
the sediment cannot be effectively released through the reservoir or put to 
useful purposes. It is inevitable, in most cases, that the useful storage ca- 
pacity behind a dam will be reduced in time by sediment accumulation until 
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part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 

a. Presented at meeting of ASCE, New York, N. Y., October, 1957. 

* Head, Sedimentation Section, Hydrology Branch, Bureau of Reclamation, 
U. S. Dept. of the Interior, Denver, Colo. 

** Asst. Head, Sedimentation Section, Hydrology Branch, Bureau of Reclama- 


tion, U. S. Dept. of the Interior, Denver, Colo. 
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replacement storage is necessary. In many countries reservoirs construct- 
ed years ago have become completely filled with sediment deposits rendering” 
them practically useless. As most large reservoirs in this country are rela- 
tively new, only a small number have yet been substantially reduced in ca- 
pacity. However, future generations will be faced with replacing storage 
capacity lost to sediment in our major reservoirs. Replacement of many 
reservoirs may not be easy since the most economical reservoir sites are 
currently being used. Removal of sediment from a reservoir by mechanical © 
means is not economically feasible at present but increases in the value of 
an acre-foot of water may change this. Various methods have been suggested 
and tried for removing, reducing, or by-passing sediments. Although none 
has yet supplied us with a solution, particularly in large reservoirs, perhaps 
some day when the need becomes acute a method will evolve. Of course, the — 
economic value of a reservoir may become less with new methods of power 
production and water supply. 


The Problem 


In present day analysis we seldom look beyond one or two generations in 
planning the water supply to be provided by a structure and often the struc- 
ture size is dictated by available water supply and economical considerations. 
Indeed the pay out period on all projects built by agencies such as the Bureau 
of Reclamation is usually less than 100 years. However, in planning the 
storage allocations it is the Bureau’s policy to provide for 100 years of sedi- 
mentation unless it can be demonstrated by economic analysis that a lesser 
allotment is justified. 

There was a time when it was believed that sediment always deposited in 
the bottom elevations of a reservoir rather than depositing throughout the 
full range of reservoir depth. It is now fully realized that sediment deposits 
will spread throughout the reservoir reducing the incremental capacity at all © 
elevations (Figure 1). Factors that influence the mode of deposition include: 


1. Reservoir shape 

. Sediment characteristics 

. Reservoir operation 

. Sediment-reservoir volume ratio 
. Inflow-capacity relationship 


oP & DO 


There are other factors, of course, such as salt cedar infestation, amount 
and location of tributary inflow within the reservoir, etc., but these, when 
existing, require special considerations. 

Of immediate concern in planning and designing a dam and reservoir, is 
the effect of sediment accumulation on design considerations. Major design 
considerations include: 


1. Elevation to which sediment will accumulate at the dam in a given 
period of time. This information will affect the design elevation of the river 
outlet and power penstock sills. Sometimes an adjustable stoplog type trash- 
rack arrangement is used which provides flexibility in the outlet elevation to 
allow complete use of the water storage at alltimes. Stoplogs are added as 
the sediment accumulates. 


2. Reduction of active storage capacity. If this loss can be predicted the ; 
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ructure height can be set to provide the required water supply storage. 
iten, the effect of sediment accumulation on the water supply is evaluated 
rough the operation studies. 


3. Deposits of sediment in the upper reaches of the reservoir.(1) These 
posits can result in a significant increase in the water surface elevations 
stream. This will, in turn, endanger upstream installations and develop- 
ents, particularly during periods of high water inflow. The dangerous con- 
tion results because the reservoir water surface is usually high when the 
flowing stream is at flood stage. This condition must be considered in the 
signs. 


4, Effect on recreational developments around the reservoir. The sedi- 
ent deposits could influence the location of such facilities, especially when 
e reservoir is one that will have considerable drawdown exposing sediment 
‘posits. 


5. Reservoir planned solely for sediment accumulation. In this case the 
‘obable sediment accumulation is of primary importance. A considerable 
tion of the total sediment stored may be above crest elevation of the 
‘servoir. 

It is therefore, evident that estimating only the amount of sediment inflow 

a reservoir is not sufficient. We must try to predict how and where this 
diment will deposit. 


Methods for Predicting Sediment Distribution 


Having discussed the need for predicting the manner in which sediment 
ill be distributed, two methods that can be used in making this prediction 
ill now be presented and illustrated. The first method which is strictly 
athematical is known as the Area-Increment Method, and was developed by 
_ A. Cristofano.* The second method, which is a mathematical procedure 
veloped on the basis of actual occurrences in large reservoirs, is identified 
rein as the Empirical Area-Reduction Method.(2) Both procedures involve 
e adjustment of the original surface areas to reflect the decrease in area 
ith sedimentation. (Figure 2). 


Area Increment Method 


The area increment procedure is based on the assumption that the sedi- 
ent deposition in a reservoir can be approximated by reducing the reservoir 
ea at each reservoir elevation by a fixed amount. This involves a series of 
proximations. By using the average end area or prismoidal formulas the 
Servoir capacities are computed on the basis of reduced surface areas un- 
the total reservoir capacity below maximum normal water surface is 
sentially the same as the predetermined capacity obtained by subtracting 
6 sediment accumulation with time from the original capacity. Using 
amogordo Reservoir on the Pecos River in New Mexico(3) as an example, 
= step procedure for the Area-Increment Method follows. The procedure 


Former Hydraulic Engineer, Project Office, Bureau of Reclamation, 
Albuquerque, New Mexico. 
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illustrated is an adaptation developed by W. T. Moody.* The basic infor 
tion for Alamogordo Reservoir as obtained in the 1944 resurvey is: 


Original capacity at El. 4275—156,750 acre-feet 

Capacity after sedimentation—132,170 acre-feet 

Sediment accumulation—24,580 acre-feet (1936-1944) 
Original Reservoir depth at dam—125 feet (El. 4150 to 4275) 
Elevation of sediment at dam in 1944—4190 (approximately) 


mcr mate ill ayy LD, samntidied 


Step Procedure—Area-Increment Method (See Table 1) 


The area increment procedure can be expressed by the following basic — 
equation: 


Vg = Ao (H - ho) + Vo : 


. . : : ‘ 
where: Aj = area correction factor in acres which is the original reser- 
voir area at the new zero elevation at the dam, 


Vo = sediment volume below new zero elevation (acre-feet), ; 


Vg, = sediment volume to be distributed in the reservoir (acre-_ 


feet), " 


H = reservoir depth at the dam—streambed to maximum norma 
water surface (in feet), 


ho = depth in feet to which reservoir is completely filled with a} 
sediment—new zero elevation 


This equation states mathematically that the total volume of sediment Vg, Hl 
consists of the portion which is uniformly distributed vertically over the 
height H - ho, plus the portion below the new zero elevation of the reserv 

Using this equation the step procedure for the Area-Increment Method 


Step 1. From the equation, ho can be obtained by trial and error as iol = 
lows: 


Vg = 24,580 acre-feet 
H = 125 feet 


Vo, hg, and Ap are obtained by assuming h, and reading Vo and Ag fro: 
the original area-capacity curves (Figure 6) at that indicated elevation. 


First trial—assume ho = 25 feet 
then Ao = 150 acres Vo = 1,600 acre-feet 
so 24,580 = 150 (125-25) + 1,600 + 16,600 


Second trial—assume h, = 34 feet 
then re = 235 acres Vo = 3,300 acre-feet 
so 24,580 = [ 235 (125-34) + 3 ,300] = 24,680 


The new zero elevation at dam = 4150 + 34 = 4184 
Area correction factor = 235 acres 


Step 2. Compute accumulative sediment volume by applying the area c 


* Engineer, Technical Engineering Analysis Section, Division of Design, 
Bureau of Reclamation, Denver, Colorado. 
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rrection factor at each depth increment (Col. 4) and computing volumes by 
erage end-area formula (Col. 5). Computed volume should check the pre- 
termined sediment volume within 1 percent. 


Step 3. Reduce original area at each increment by the area correction 
tor of 235 acres to give revised area table, (Col. 6 = Col. 2 - Col. 4). 


Step 4. Determine the revised capacity table by reducing the original 
pacity at each increment by the sediment accumulation, (Col. 7 = Col. 3 - 
m5). 


The results obtained are compared with the actual resurvey curves in 
igure 6). 

The basic equation for this method is very useful in making a quick esti- 
ite of the depth of sediment at a dam over a given time period. By refer- 
ig to the original area-capacity curves and knowing the amount of sediment 
it will accumulate in x years, we can, with a few approximations, balance 

> equation and thereby derive the height to which sediment will accumulate 
the dam. Once this is established the area reduction can readily be applied 
each interval of depth up to maximum normal water surface and the proba- 
} area-Capacity curves after sedimentation can be determined. 

Tt should be noted that there is no smooth reduction in area from the 
iginal to the revised condition between maximum normal water surface and 
} next increment above in the flood control pool. This can be corrected, if 
sired, by allowing some sediment to deposit above maximum normal water 
rface and by making the area correction factor progressively smaller 
irting a few intervals below the maximum normal water surface elevation 
i extending a few intervals above. The surface area for the before and 

er condition should be the same at the uppermost elevation. Some addi- 
nal trial and error will be necessary here to get the computed sediment 
lume to equal the predetermined volume. 

The Area-Increment Method is only an approximation of the probable loca- 
n of the sediment deposits. It is known that sediment does not generally 
jJuce the surface areas by the same amount at every elevation. However, 
mparisons have been made between the results obtained by this method and 
- actual resurvey results and good agreement has been obtained in most 
ses. It should be noted that the applicability of the method decreases with 
increase in the sediment/capacity ratio. Also the accuracy decreases 

h any deviation of the reservoir in question from a more or less standard 
e. Usually the arbitrary stipulation is that if the 100-year sediment 
sumulation exceeds 15 percent of the original capacity a more exacting 
thod of distributing the sediment should be applied. In any case, the area 
rement method is a good first check procedure for determining the prob- 
e sediment distribution. 


Empirical Area-Reduction Method 
The Empirical Area-Reduction Method for determining the probable sedi- 
nt distribution is accomplished through two main steps: 


1. Classify the reservoir using 4 basic standard type curves which were 
eloped from actual resurvey data. 


2. Make a trial and error type computation using the average end area or 
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prismoidal formulas until the capacity computed equals the predetermined ; 
capacity. 4 


The resurvey data for 30 reservoirs have been used to develop 4 standard 
type curves of percent sediment deposit versus percent reservoir depth 
(Table 2 and Figure 3). The reservoirs included in establishing the 4 stand-- 
ard types vary widely in location, type operation, and sediment characteris-= | 
tics. Capacities range from 40,000 to 30,000,000 acre-feet (with a few ex- ; 
ceptions). The depth is based on the difference between streambed elevation : 
at the dam (zero capacity) and elevation at maximum normal water surface, 
The general classification resulting from the analysis of the data from the 
30 reservoirs is as follows: 


4 


Reservoir Standard 
“M” type classification 

1.0-1.5 Gorge IV 
1.5-2.5 Hill il 
2.5-3.5 Flood plain-foothill I 
3.5-4.5 Lake i 

(“M” is the reciprocal of the slope of the line obtained by 
plotting reservoir depth as ordinate against reservoir capacity Pi 


as abscissa on log-log paper (See Figure 4). 

It should be noted that a lake type does not necessarily have 
to be out on the plains nor does a gorge type have to be in the 
mountains. The gorge type (IV) is the poorest defined of the four 
classifications.) 


In applying the standard type curves to a proposed reservoir, it is first — 
necessary to determine the type that fits the reservoir involved. First | 
selection of type is obtained by plotting depth versus initial reservoir ca- b 
pacity on log-log paper (Figure 4). This usually gives a straight line rela- — 
tionship or in some cases 2 straight lines.(4) The slope of the line (or lines) 
gives an indication of the type. Generally a curve showing rapid increase in 
capacity with depth indicates a large basin type reservoir whereas a small if 
increase in capacity with depth indicates a gorge type reservoir. This in- 4 
formation classifies the reservoir according to shape. The standard type 
curve can then be selected and adjusted as necessary to reflect other condi- 
tions, such as anticipated reservoir operation, type sediment inflow, or 
capacity/inflow ratio. The type curve of percent depth versus percent sedi : 
ment deposit finally adopted may be a combination of 2 types or could fall i 
between 2 of the 4 basic types, but generally 1 of the 4 basic types will fit the 
conditions involved. It can be seen that considerable judgment is sometimes 
involved in this selection. 

Having selected a type sediment deposit curve the rest of the procedure i 
mathematical. The 4 standard type sediment versus depth curves have bee! 
converted to area design curves for use in the computations (Figure 5). The 
area design curves have relative sediment storage as the ordinate and rela- 
tive depth as the abscissa. The area under the curve in each case equals 
unity. Other curves can be produced as long as the area under the curve is 
unity. The conversion from the standard type sediment deposit curves to th 
area design curves has been made by Moody by applying the equation: 


Ap =C p™(1-p)" 
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where: Ap represents a dimensionless relative area at relative distance 
Pp above streambed 
C, m, and n, are dimensionless constants which 
are determined by the type of reservoir 


merical values of m, and n, were computed by trial and error using a 

ast Square procedure so as to make the analytical curve the best possible 
_to the field survey data. With m, and n, determined C is fixed by the con- 
deration that the total area under the Ap curve must be unity. The charac- 
ristic constants C, m, and n doterwiney for the 4 types of reservoirs are 


follows: 
Type Cc m i Sediment Storage 
near 
I 3.4170 1.5 0.2 Top 
0 2.3240 0.5 0.4 Upper middle 
Il 15.882 el 2.3 Lower middle 
IV 4.2324 0.1 2.5 Bottom 


ere is no logical progression of these constants as the silt factor moves 
om top to bottom. 

The procedure from this point can be handled in tabular form. Again 
amogordo Reservoir data are used in the illustration. The steps are as 
lows. 


Step Procedure—Empirical Area Reduction Method (See Table 3) 


Step 1. Determine the relative depth for each increment at the dam (Col. 

. This is merely the ratio of the incremental depth to the total depth (maxi- 
im normal water surface to streambed—El. 4275 to 4150). This is deter- 
ined from values in Col. 1. 


Step 2. Determine values of relative sediment area (Ap) from the standard 
xe curve selected (Figures 3 and 5). For Alamogordo a type II has been 
termined. Enter the curve with relative depth for each increment and read 
rresponding Ap values. (Col. 5). 


Step 3. Select a first approximation of the probable sediment elevation at 
> dam after sedimentation. Areas at and below this elevation will equal 
eas in Col. 2. Sediment areas for each depth increment above the estimat- 
new zero elevation are obtained by dividing the original area at zero ele- 
tion (Col. 2) by the corresponding Ap value (Col. 5) and multiplying this 

tio (K) by the Ap values at each succeeding increment. For example—re- 
‘ring to Table Z: 


Assume zero elevation of 4190 
Surface area = A, = 300 acres 
Ap at El. 4190 = 10125 
300 4 
K= 1125 = 266. 
The new area at each succeeding elevation is the Ap at that elevation 
aes 266.7. 
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Step 4. With the sediment areas established (Col. 6) the incremental sedi- 
ment volumes can be computed by the average end-area formula and entered | 
in Column 7. Columns 6 and 7 represent the first approximation. If the sedi. 
ment volume summation in Column 7 exceeds or is less than that desired a 
second approximation is made. Second and third approximations are shown 
in Columns 8 through 11. Note that the new zero elevation may not occur 
an even contour interval. In the example shown the new zero elevationis — 
4183.5. When making the evaluation at odd increments the original area 
capacity curves (Figure 6) and Figure 5 are referred to. 


Step 5. Accumulate the sediment volumes shown in Col. 11 and enter in ; 
Col. 12. 


q 
Step 6. Determine the revised areas for after sedimentation (Col. 13 = — 
Col. 2 - Col. 10). 


Step 7. Determine the revised capacities (Col. 14 = Col. 3 - Col. 12). 


Usually 2 or 3 trials are all that are necessary to obtain the desired 
balance. If the summed volume is within 1 percent of the predetermined 
volume, sufficient accuracy has been attained. 

There are several advantages to the Empirical Area-Reduction Method. — 
It permits the classification of a proposed reservoir into a standard type 
that immediately gives an indication of the manner in which the sediment will 
distribute. A portion of the sediment can readily be deposited above maxi-— 
mum normal water surface (flood control space) in the computations, by ’ 
selecting a top elevation above maximum normal water surface in determin- 
ing relative depths. (In a large reservoir the sediment deposits above maxi- 
mum normal water surface will seldom exceed 10 percent of the total ina 
50- to 100-year period, except for those dams which are for sediment reten- 


tion only.) It gives a very good approximation of the height of sediment at “ 
dam at various time intervals. 


of 
1 


CONCLUSION 


A comparison of the Area-Increment, Empirical Area-Reduction Method 
and actual reservoir resurvey area and capacity curves demonstrates the 
adequacy of the methods, (Figure 6). However, in the two methods demon- 
strated it is evident that a procedure in which all the major factors affecting 
sediment disposition is still needed. Also, several major factors affecting — 
reservoir sediment distribution, such as density currents, delta development 
and density of deposits with time, but not discussed herein, need additional 
study and evaluation to develop procedures for predicting the effect of these: 
items on conditions in and upstream from the reservoir. 

Current work of the ASCE Task Force on Reservoir Sediment Distributio 
includes an evaluation of the effect of the various factors on sediment deposi: 
tion and distribution, and the nature of delta development. Investigations int 
these problems are also being carried out by agencies such as the Corps of 
Engineers, U. S. Geological Survey, and Bureau of Reclamation. Factors in- 
fluencing the sediment deposition in numerous reservoirs are being collectet 
consolidated, and analyzed, and mathematical relationships and criteria will 
be developed therefrom to give procedures that will yield consistently reliab 
results in predicting sediment disposition in proposed and existing reservoil 
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Table 1 


SEDIMENT DEPOSITION COMPUTATIONS 


Project: Carlsbad 

Reservoir: Alamogordo 

Stream-bed elevation at Dam: 4150 Annual sediment inflow: 3,600 ac.ft. Computed by: CRM Date 
Max. water surface elevation: 4275 Period of sedimentation: 6.8 years 6-15-54 

Spillway crest elevation: 2 Total sediment for period: 24 ac.ft. Checked by: JML Date 6 


AREA INCREMENT METHOD 


Flevation _orietnal Pee PT kien La 
(feet) (Acre feet) (Acre) (Acre feet) (Acre) 
1 2 3 4 5 6 7 { 
275 y 156,750 5 , 15 i3e, 
4270 4100 133,500 235 23,510 3865 110, 
4260 3200 97,000 235 21,160 2965 755 
4250 2450 68,750 235 18,810 2215 4g, 
42ko 1750 47,750 235 16,460 1515 31, { 
4230 1250 32,750 235 14,110 1015 18, 
4220 900 22,000 235 11,760 665 10, 
4210 650 14,250 235 9,410 415 4, 
4200 450 8,750 235 7,060 215 1,690 
4190 300 5,000 235 4,710 65 290 
#418) 235 3,300 235 3,300 ° 
4180 200 2,500 200 2,500 ) 0 
4170 100 1,000 “100 1,000 to} oO 
4160 50 250 50 250 ° 
4150 0 ° fo} ° ° 


*Determined new zero elevation 


AREA INCREMENT PROCEDURE COLUMN DESCRIPTIONS 
Basic data 
Vs = 24,580 ac.ft. (sediment volume) 1, 2, 3. Basic data. 
H = 125 ft. 

Vo, ho, Ao = obtained by approximating ho 4, Area correction factor 
and reading Vo and Ao from determined above. 
original area - capacity 
curves (Fig. 12) 5. Sediment volume as 

determined by average 
Basic equation end area formula, as 


based on area correcti 


Vs = Ao(H-ho) + Vo factor in Column 4. 


assume ho = 25 feet 6. Revised areas obtained b 
Ao = 150 acres subtracting Colum 
Vo = 1600 acre feet from Column 2. _ 

then 24,580 = (150)(125-25) + 1600 7. Revised capacity obtai 


24,580 16,600 by subtracting Colum 


from Column 3. 


assume ho = 34 feet 
Ao = 235 acres 
Vo = 3,300 acre feet 


then 24,580 = 235(125-34) + 3,300 
24,580 24680 


The area correction factor is therefore 
235 acres and the new zero elevation is 4184.0. 
The resulting area and capacity data are show 
in Columns 6 and 7. 
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Table 3 


ALAMOGORDO RESERVOIR 
Sediment Deposition 
Empirical Area Reduction Method 


3 3 ee ie : :First:Tri 3 cond: ia : ird = ie : 3 : 
: z Alea : 3 3 oss 3 76 : oe om H ss 
Se od 3 d 3 By cova ERG : e: : 3 : : S258. & : 7 
aeordc no * Ce ed ® H te Pee o » we . 2 as oc ov" z 2 vg & 
2 ent aoo . a0 4 s tS es re a: © ws c @; ce sw sc 9: cg s: aH cl 3 = 
eS REE: Bae Gut ae ate whe, OSe, Wee ees ae ‘ale ary | 
s- BSS 58s ge: eB ipehgeg: aed: geg seas. gee 2eee ok ee 
a : : aS : ~ 30 :a7S:a8 rece $e gt = - = 
(2) 2) GB) Se) 5) 6) ee) a) Oe oye OTe ee (13) 7 (24) 
34,550 : 2100 31.00 3:0 20 5 nu ; ooh 2 324, 24,050 =: 
F : : : > M5: £ 92501 3 3. B34. 38 3 7 
4270 :4,100° :133,500 : .96 : .c2 :1€6 : Fa + Tavs 224,180 23,970 :109,300 
: _ 3 a : 32,070 : 71,615 > 31,675" 3 3 3 
4260 :3,200 : 97,000 : .88 : .93 :eh8 : : igh : $a = 222,505 :3,000 : 74,50 
q fi fz : A 32,690 : 32,105" = 32,190 ¢: ¢ 3 | 
4250 :2,450 : 68,750 : .60 :1.090 :290 : ey ae a a 320,315 :2,c10 : 46, 
: : : : 2 33,030 ; [2,310 2 = cg k65 | x 3 : 
kekho 31,750 ; UR, TSO 3 Te 32.185 3316 = : 2k7 . : 256 3 317,850 :1,490 : 29,900 
i : A E . 23,220 ; t2,520,. 2 32,615 3 & 3 
4230 :1,250 : 32,750 : -64 :1.235 :328 : 225i bis 2/267) > 315,235 : 980 = ifpamm 
: : ; G : 23,310 : 82,590 ¢ 32,6490 : ; i 
4220 ; 900 : 42,000 ; .5€ :1.255 :334 : ft26h) | fe Pek 222,545 3 630) 5 ay 
é : : : : 33,320 ;: cp As = ae 22,695 : : : 
4210; «4650 : 14,250 ; .4B :1.2h0 :330 : : 258 : ShOBO 1 ts : 9,850 : 380 : 4&,kOO 
: & : : t 23,240: 32,535 3 22,630 : : : 
4200; 450 : 8,750 : .-4O 31.195 :318 : : eho: 7258: : 7,220 : 190 =: 1,5 
4 : : 7 t 23,090 : 22,415 22,505 : 3 : 
4190 : 300 : 5,000 : .32 31.125 :300 : eiegie: : 2h3 25 = es 
: ‘ o F : Seg OO os Seal pace se5eI5 > : $ 
4180 : 200 : 24,500 ; .ch 31.02 :.00 ;: : 200 ; 2: 20OO : 2,500 : o ¢ 
A fi : : = 31,500 : 21,500 : 31,500 : : ¥ 
4170 : 100 ; 1,000 : .l6 : .87 :100 : : 100 : : 100 :; : 1,000 : oO g 
: ts i : a 3. laws 3 wae) & es eee 3 tu 
4160 ; 50 : 250: 08 +: 63 +350 : se SO TS t 500s :. 250 3: QO 
fe “ : g H Seas e250) oes 5 PE) 7s : = 
4150; O'S: Ome O" pers 108 Ose Ona CP es 2) BOF : QO: (an 
: : : : : £9,305: Br: et at : : 
Supplement : 3 3 3 3 - 3 : : : 3 : 
4183; : : .264 31.06 ; : : : : : : : : 
4183.5; : : .268 31.065 : : : ie Z : 5 3 : 


Sediment to be deposited = 24,600 acre-feet 
Standard Type II selected 


Column Descriptions 


1, 2, 3. From original area-capacity data. x 
4, Relative depth--Ratio of contour interval depths to total depth (spillway crest to stream bed) from 
elevations in Column l. : 
5. Ap--Values obtained from Figure 5 after proper standard-type curve (Figure 3) has been selected. 
Enter curve with relative depth values to obtain Ap values. 
€,8,10. Approximations of zero elevation at the dam after sedimentation. Areas at and below approximated 
zero elevation et each increment will be equal to those in Column 2. New areas for each contour 
elevation above zero elevation are obtained by dividing the original area at zero elevation 
(Column 2) by the corresponding Ap value (Column 5) end miltiplying this ratio (K) by the Ap: 
values at each contour elevation above. 


Example: Assumed zero El. = 4190 
Surface area = 300 acres 
Ap at El. 4190 = 1.125 


= 00 = ~ 
K as 266.7 


Note that zero elevation may not occur at an even contour interval (see supplement data) 
7,9, 11. Delta sediment volume determined from average-end area formula: 
n/2 (Aj+ Az) =V 
Summation of this column should equal (within 1 percent) the total sediment being deposited. 
Accumulated sediment volumes from Column 11. 


Revised areas obtained by subtracting values in Column 10 from those in Column 2. 
Revised capacity obtuined by subtracting values in Column 12 from those in Column 3. 


PRE 
wh 


1587-13 


DISTRIBUTION OF SEDIMENT 


SCE 


|) Ns} 


"Sjajino 1p Huloinjs Aq paulojulow si abouojJs ajqojinAo Huluiowes 


-puo ouaz Sawooeg Adualdisja dos} abofs siut ty 
“JUaWIpas ae Pajif JIOASaSay 


zy) 
0) 
"UMOPMDIJP JIOAJOSON ajqOJaPISUOD PUD “UMOPMDIP JIOAJaSaJ afOuapowW AjjowJOU puD 
‘sazis ajdijuod Aojo pun 4y!S ul AjpuoUIWOp ‘SOZIS aj|{u0d PUDS Ul} PUD 4)IS Ul Aj,UDUIWOp 
“fuawiipas ae pally Ajjolyund wlovuasay “LUBWIPAS YLIM palfly AjjDlyuod wioAuasay 


“YOLOM JO |/NJ AjjoW4OU UlOAUaSes puDd 
‘SaZzIS ajdlyu0d aSupo0d ul AjjuoulWwOp ‘paisodap juawipas on 
*} Ua! pas ULIM palit Ajjolysod JIOAJasay SUOI4IPUOD |DIYIUI JApUN JIOAJESeY 


Py 
_--P2quDad4s 
EA \e 


D fO9LJ9 JALDMYODG --* 
$}9]41NO aHdDIOJS podeq aHDJOLS IAILODU| dHDJO4S UOIL DAU8SU04 99DdS JOULU0D poo|4 
2555 p J q D 


NOILVNV IdX4 
SNOILIGNOD NOILVLNAWIGSS SNOIYVA YSGNN SYIOAYSSSY INIMOHS 


AIMLTIEL WALI CAI! NII WIALIALWiIOAWIN 


catia - seliadaamaiamee 


| L 7 
2 PRN 
WDdd4s |DUIBIUQ-~___ 
ey ‘ 
pousaw uolponpay 
- DY “a, | 
. Poyjaw juawauguy - Dady-~- Yy 
ir AA 
LLEEz ay 
On 4 
aos Ww NON 
V4, Zio fa) 
: YYZ ons VY NI LNSWIGSS JO NOILNEIYLSIG 
x Zz, SNILIIGSYd JO SGOHLIW LNSW3SYONI-V3NV 
S JHL GNV NOILONGSY-V3SYV IWOIYIdWS 
3 JHL 30 NOSIYVdWOD DJILVWWYYSVIG 


1587-15 


DISTRIBUTION OF SEDIMENT 


ASCE 


Ol 


02 


O£ 


06 


Q3L1ISOd3aqd LNS3WIGSS LN390Y3d 


00! 06 08 OL 09 0S Ov 0¢ 02 Ol 


JO SAAUND AdAL 


cara 


RMGmerMres 
BRAM Sess Se’ 


SV ACR eS 
W/ a cavawaile sce en! 
Pole bo katahrusiatae 
[poet bidokay betel cl 


FEM 


Hey 


0 


Ol 


oO 
3 i 


oO 
t+ 


Hid3d YIOANSS3SY LN39YAd 


(1$3Y9 AVMI1IdS OL G3IEGWVayLs) 


06 


HY 2 


1587-16 


TL adhy 


(111H) 


Tr adhy~ 


-82'0 Of Ob'0 = (U) ado| ee 
SY Of Gt * i, te 
220 04 82'0 = (u) 7 LA 


(FTIHLOOS-NIVId GOO74) 


dIHSNOILV 134 


ALIOVdVO SNSYAA H1d5d 
Ag 


SYIOAYSS3SY JO NOILVOISISSV 19 


(9) ALlovdwo 


Se Oh Sg = \w, 


S’¢ Of 1 "ow ' 
Ob'0 04 29'°0= (u) adojs-” 


Hld4jd 


(U) 


SCE DISTRIBUTION OF SEDIMENT 1587-17 


om IW Ap= 4.2324 pl |-p)25 mise: 
eee 


|X) 


AUTURC Meee 
ak oe 
V4 Me 


Geet 
Ws RELATIVE DEPTH (p) \_| a 
(MEASURED FROM Ro EN 
Oz2. 0.4 0.6 0.8 1.0 
SEDIMENT DISTRIBUTION 


AREA DESIGN CURVES 
(BASED ON RESERVOIR STORAGE CURVES ) 


FtGss 


HY 2 


1587-18 


4334 - NOILVA313 


Oeeb 


Ofep 


Over 


OS2r 


YSAIY SO03d 
YIOAYSSIY OGYOSOWYV IV 


S3AUND ALIOVdVO - V3YV G3SIAZY 


NOSIYVdWO9 


(pp6l) ASAYNSAY 
GOHLIW LNSW3YONI V3YuV 
JAYND ALIOVIVD TWNIDIYO 


GOHLIW NOWONGSY Vay WOI4idW3 ——»— 


NOILVNV1dX3 


a a ae AAUND WAYV IWNIDINO 
JAYND ALIOVdVD TWNIDIYO 


47 


Ober 


OS2b 


09¢b 


Oley 


0 00¢ 


008 


002! 009! 0002 O0r2 0082 ooze oo9¢ 000b 
S3udV - VauV 


O0bb 


008+ 


O82r 


41334 - NOILVA313 


aper 1588 April, 1958 HY 2 
ee 
Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 
2 


WATER DISTRIBUTION DESIGN AND THE 
McILROY NETWORK ANALYZER* 
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and J. V. Radziul*** 
(Proc. Paper 1588) 


SYNOPSIS 


Illustrates time-distribution of demand rates, design rates and equalizing 
Orage requirements for certain Philadelphia districts. Efficacy of 
rfectly-balanced system characteristics and utilization of the McIlroy 
lalyzer for design in Philadelphia are discussed. Merits of different types 
computers are compared. 


INTRODUCTION 


Inasmuch as the conditions and requirements of a given water system are 
nerally somewhat unique, samples of data for some of the high service 
stricts in the City of Philadelphia will be presented as a background for 
> descriptions of design practices and experiences which follow. After 
scussing design experiences with the McIlroy Network Analyzer, compari- 
ns will be made with digital and analog computers, from the standpoint of 
e by both utilities and consultants. 


tes Used in Distribution Design and Their Justification 


In most of the discussion of design criteria in this paper, data from a 
cently completed study of the City of Philadelphia’s Roxborough High 


te: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1588 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 
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‘Research Engr., Philadelphia Water Dept., Philadelphia, Pa. ; 
‘Water Distribution Design Engr., Philadelphia Water Dept., Philadelphia, 
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Service District will be used as an example. The Roxborough High Service 
District is predominantly residential, with a moderate mercantile develop- | 
ment and practically no industry. 

The average Roxborough consumption rates for each day in 1955 and 1956 
have been arranged in order of magnitude and appear in Figure 1. The t 
“normal,” or demand equalled or exceeded half of a year, is very close to : 
the yearly mean, or average, for both years. The maximum-day demand of © 
record occurred in 1955; the 1955 data are generally more diverse. t 

Hourly demand variations over a 24-hour period for both a typical day 
with an average demand equal to the 1955 average annual demand, andthe ~ 
maximum day of 1955, are shown in Figure 2. Here, the normal is definitely 
greater than the mean or average. The curves are similar in shape. Note 
that demands included between 0.8 and 1.2 of the average demand in both 
cases are representative of about 12-hours (50%) of the 24-hour period, and 
rates above 1.2 of the average demand for each case are representative of 
only about 6-hours (25%). , 

Figure 3 is to some degree a composite of the previous two figures, for 
1955. In addition to the daily average demand for the 365 days, the maximum | 
and minimum hour demands have been graphed. To be especially noted is 
the fact that variations in daily average demands within 90% to 110% of the 
average annual demand cover 83%, or about 300 days, of the year. 

Because of the moderate variations about the average annual demand for 
the major share of the year, design has been based primarily upon the maxi- 
mum and minimum hour demands of the maximum and average day. Since 
Philadelphia is supplied entirely from river sources by pumpage (annual 1 
average consumption, 370 MGD), efficient pumping at average day rates is i 
particularly important in all districts. The maximum hour of the maximum _ 
day is normally the critical demand in terms of service pressures at points © 
of high elevation. if 

The variations, daily and hourly in Figures 1 to 3, give a good indication ~ 
of intermediate pump needs for adequate coverage of demand fluctuations. _ 


Uniform Pumping Storage Requirements 


In an effort to provide more uniform service pressures and to improve 
supply reliability, equalizing storage is planned for all high service distric 
One of the major objectives in providing equalizing storage is to permit 
pumping rates as near as possible to constant or uniform 24-hour rates. 

Figure 4 shows the storage requirements for five districts recently 
studied. Average day and maximum day storage requirements are similar 
in trend. In Figure 4(b) the greatest equalizing storage requirement shown 
is 14.8% of the maximum day demand. The zero storage pumping ratios 
represent maximum hour demand rates. In contrast, Schmid(1) reports a 
25% uniform pumping equalizing storage need for Royal Oak, Michigan, for 
a maximum day demand of 26.4 MGD (where the maximum hour demand is 
1.99 times the maximum day demand). Percentagewise, much less storage 
is needed for uniform pumping in the Philadelphia districts shown. 


Utilization of the McIlroy Analyzer in Distribution Design 


There are approximately 2900 miles of pipe in the system, ranging in siz 
from 4'' to 93''. However, the McIlroy Analyzer of the Philadelphia Water 
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epartment, installed in January, 1956, has been used mostly in connection 
ith arterial mains. Field measurements of Hazen-Williams C-values for 
‘terial mains are made by the Pitometer Company. To date, the rate of ob- 
ining field measurements has set the pace in rehabilitation design. 
Distribution rehabilitation design includes piping, equalizing storage, 
imping stations, routing and districting. Many design variables and sys- 
m complexities become involved. For example, in the Roxborough High 
srvice District previously mentioned, the existing district is presently fed 
y one pumping station; revisions include two standpipes recently constructed 
id a third future standpipe at another location, in addition to a second, en- 
rely new, pumping station and a modification of the district boundary. 
Many of these revisions are necessary for abandonment of an existing filter 
ant). 
The following conditions were investigated using the analyzer for present 
id future demand rates in the Roxborough High Service arterial system: 


1) Minimum and Maximum Hour of Maximum Day - latter usually deter- 
mines reinforcement needs (new mains or cleaning and lining). 


2) Minimum and Maximum Hour of Average Day. 
3) Arterial pressures for fire flows. 
4) Pumping station outages. 

or present demand rates, also investigated were: 


5) Effects of premature change in district boundary (prior to new, second 
pumping station), and 


6) Schedule for satisfactory service under which mains can be laid or 
cleaned and cement lined, connections can be made and valves installed 
or rehabilitated. 


Uniform pumping was one of the major design objectives, from the stand- 
int of improved service, efficiency and amenability to automatic control. 
orage volumes and rates were computed for perfect 24-hour uniform pump- 
g. Using constant pumping rates and the computed storage water levels, 
imping station discharge heads were obtained with the Analyzer (after piping 
rstem had been designed). The resulting discharge heads for one of the two 
xxborough High Service Pumping Stations is shown in Figure 5, plotted in 
rms of pump T.D.H. Hypothetical pumps are shown for the purpose of illus- 
ation. On the maximum hour of the average day (lower point) a greater 
lantity than the average rate will be pumped, with a lesser quantity than 
mputed coming out of storage. On the minimum hour of the average day 
pper point) a lesser quantity than the average rate will be pumped, with a 
sser quantity than computed going to storage. Adjusting the storage curve 
this instance would not noticeably change the storage volume used, and 
uld modify only slightly the water levels computed for uniform pumping. 

If the curve for pump No. 1 was raised, then less than the storage volume 
mputed for uniform pumping would be used; for a lowered curve, more 
orage volume would be used. For the maximum day under practically all 
nditions shown, the pumping rate would be less than the average of 9.5 MGD, 
d the difference would have to be made up by the second pumping station or 


'a storage depletion. 
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¢ 
In Figure 6 is plotted part of the maximum day storage output-versus- 
water-level data for three Roxborough High Service standpipes, computed 
for uniform pumping and assuming a specific filling schedule. The Maximum: 
Hour rate out of storage with uniform pumping is 8.0 MGD. In the table below 
the graph, are the T.D.H. for pumps of each of the two stations, as designed, 
for uniform rates on the Maximum Hour of the Maximum Day (the second and} 
third cases are outages of each of the two stations). In the “initial balance” 
columns are shown T.D.H. and rates for the sources indicated, which were — 
obtained by extrapolating system test points to the pump curves; the storage ~ 
rates shown would, by default, be required to balance the rates. A revised — 
storage rate requires a new storage curve. Since uniform pumping can be ; 
approached but seldom perfectly attained, the only way in which a perfectly 4 
balanced system analysis could be achieved would be by making a series of 3 
runs starting at an assumed full storage time, following the demand curve ~ 
hour by hour and balancing heads with pump characteristics. This is much © 
too involved for the minor modifications which would be obtained in this 
instance. The balancing of pumped systems with and-without equalizing stor- 
age will be discussed later in terms of computer operation and analysis. 


en - 


Other Computers for Analysing Distribution Systems 4 


With the advent of the McIlroy Analyzer and the successful programming ~ 
of digital computers for network analyses, any discussion of approximation 
methods should be confined to the simplest of distribution systems. by 

Features of the McIlroy Analyzer have been presented previously.(2,3,4,56 
Suffice it to say that the principal feature of the McIlroy Analyzer is a system 
of resistances through which voltage drop (analogous to head loss) is propor= 
tional to the 1.85-power of the current (current flow is analogous to flow rat 
The McIlroy Analyzer thus provides a “model” or direct, working analog of © 
the system under study. In operating the Analyzer, “drawoffs” are set at 
prescribed full load values and system inflow rates (from pumps or storage 
are brought up from zero to system capacity. Often, in design analyses, 
needed system revisions are evident at inflow rates well below the design 
rates; obviously needed piping improvements are then made prior to any 
formal or complete test runs. Trial changes can be made in a matter of 
minutes. For example, changing the distribution of inflow between two pu 
ing stations can be accomplished by adjusting two or three dials. The most 
valuable attribute of the McIlroy Analyzer in this respect is the elimination 
of unrealistic or superfluous or discordant arrangements with quick, simpl 
direct dispatch. 

Electronic digital computers have been employed in gas and water distri 
bution network analyses. The punched card program for a low-pressure ¢: 
system using a digital computer, with a wired program, has been presented 
by Sickafoose.(7) (However, a 2.00-power rather than 1.85-power is involv 
An outline of the programming for a water network analysis using a stored 
program type of digital computer has recently been presented by Hoag and 
Weinberg.(8) In both of these references, a form of the Hardy Cross relax: 
tion technique was included in the analysis program, the computer perform- 
ing the trial-and-error iterations. The digital computer thus adjusts flow 
rates initially assumed by the designer for each pipeline into a balanced 
distribution to almost any degree of desired accuracy in loop headloss 
balance. Other tests of the same system can often be made by replacing onl 
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few of the punched cards from the initial test data instructions. The ma- 
uine itself is tied up only a very few minutes in any one run. However, 


bulation and manipulation of data instructions is tedious and time -consuming. 


urther, the results are normally presented as tabulated head loss and flow 
ates (with direction) for each pipeline, and these must be summed 
igebraically before the analysis can be evaluated. An error in programming 
r data instructions is not evident until the run is made and the analysis is 
Jmpletely summated and checked, and can be very difficult to locate. 
punched data card is needed for each loop (identifies mains and flow 
rections in mains comprising the loop), and for each main (“C”, diameter, 
mgth, initially assumed flow rate) for each loop (mains often common to two 
¢ three loops). About 50 to 60 punched cards are needed for the basic pro- 
fam instructions, almost regardless of the size of the network being studied; 
e basic program instructions can be used over and over with different net- 
orks. The possibility of errors in data instructions is evident, particularly 
ith medium to large networks. Further, the “model” visualization feature 
‘the McIlroy Analyzer is completely absent from the digital computer. 

Analog computers are being used in water-hammer analyses, (9,10) but 
svices for the simulation of the non-linear characteristic of pipelines in an 
talog circuit are cumbersome, and network analyses for either a steady or 
transient state must of necessity be restricted to simple systems. As with 
ie McIlroy Analyzer, meter adjustment and wiring of components must often 
> charged against “machine time,” since in both cases the machines are 
sually not then available for other work. (The McIlroy Analyzer is a 
ecialized type of steady flow, or passive, analog computer). 

Network analyses can be run, by a series of trial runs with intermediate 
ymputations and adjustments, on electric service power-system calculating 
yards with linear resistors. (Refer to first four references listed in 
eference (3) ). Only a very few of these boards have been used for water 
stwork analysis. 


nalyses with Balanced Pumping Station Characteristics - 
ithout Equalizing Storage 


Shown in Figure 7-A are two hypothetical pumping station pump charac- 
ristics. For direct pumpage the demand would be supplied by the combined 
itput of the two stations. For any given pumpage rate the discharge hydrau- 
¢ gradient, h, of a given pump is fixed by its characteristic curve. For an 
itput QR there is a singular hp for a given Roxborough pump. Since Qw 
ust equal the demand less Qp, by default there is a specific hy for the 
ump delivering Qw. Hence, for a given demand rate, either QR or Qw 

a function of hy-hp, or hp-hy. As an illustration, for an 18.9 MGD 
mand in Roxborough High Service (1955 Maximum Date Rate), the com- 
ned pump characteristics for the two stations have been plotted in 
igure 7-B. (A maximum difference of 20-feet is the practical limit in 
is district. “Pumps b and d” represent the only realistic combination 
this illustration, if only from the standpoint of efficiency). ; 

Information of this type (Figure 7-B) could be fed into a digital or analog 
mputer, and provided that the system including pump characteristics has 
e capacity of being balanced, a perfect balance would be achieved. With the 
rrent model of the McIlroy Analyzer, balance can be achieved only by suc- 
ssive meter adjustments, converging on pump characteristic balance 
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curves similar to Figure 7-B. From observations of direct-feed low-pressury 
gas network analyses using a McIlroy Analyzer, achievement of a balanced ~ 
system is not particularly difficult. Experience with direct-pumped system ‘ 
analyses at the Philadelphia Water Department is limited inasmuch as all ~ 
high service areas have or will have equalizing storage, and the only district 
which will be served by direct pumpage in the future has not yet been . 
analysed. : 
Where the design of new pumps is at issue, balance curves can be derived 
from network analyses and pumps selected to match or approximate the 
balance curves. , 
Analysis of a multiple gravity-supply system without equalizing storage 
would be similar, but more complicated if the water levels varied over a 
wide range. 


Analyses with Balanced Pump and Equalizing 
Storage Characteristics 


; 

In Figure 5 it was noted that a program of uniform pumping with equalizing | 
storage would not match perfectly a set of typical pump characteristic curves.: 
If the pumps and the equalizing storage are existing system components (not. 

a subject of design) and the system is inherently capable of balance (a very © 
restricted case) a step-by-step integration of the demand curve and the 
presentation of an automatically derived storage curve could be achieved 

with both a digital and an analog computer. The normal analog computer 
would be restricted to a simple system because of the complex assemblies 
necessary for each pipeline. The digital computer would be limited, in ex- 
tent of complexity of system to be analysed, to the patience and exasperation 
limits of the design engineer and to the ingenuity of the programmer. Use of 
a digital computer for a balanced system with pumped supply and equalizing 
storage has not been mentioned in the literature. 

Step-by-step integration of a 24-hour demand curve using the McIlroy 
Analyzer would be tedious and would require considerable meter adjustment 
to get each hour of demand in balance. 

In general, when the design of either the characteristics for a pumping 
station or for equalizing storage are under investigation, the necessity or 
desirability of simultaneously achieving a perfectly balanced system is 
highly questionable. This is particularly true when new pumping stations 
and new equalizing storage facilities are both to be designed. : 

In attempting a perfect system balance, using any type of machine, some 
very pertinent questions must be raised: 


1) Is the field data on pipe lengths, diameters, locations and C-values 
practically perfect? 


2) Are drawoffs as evaluated by population indices and field flow-data 
typical and near-perfect? 


3) Will or can the pumps, tanks and other improvements defined by the 
designer be exactly reproduced in procurement and construction? 


4) Will the system always be operated on identically the same schedule 
presumed by the designer ? 


Inasmuch as righteous affirmative answers to all these questions will lack 
complete conviction, attainment of a perfectly balanced system analysis is 
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ften not justified. The objective here is to promote a quality or level of 
nalyses reasonably consistent with the quality and accuracy of the data. 
‘herefore, attempting to attain a perfectly balanced system with a digital 
omputer would seldom be warranted despite the capacity of the machine to 
andle the problem. Another factor worthy of consideration is the large seg- 
nent of additional design time involved in balance refinements, using any 

ype of machine. 


ferits of Various Machines for Network Analysis 


The adaptation of electric service power-system calculating boards with 
inear resistors, by means of computation and manual adjustment iteration, 

s not very satisfactory. The calculating board is similar to a digital com- 
uter to the extent that the final result is all that is obtained; it lacks most 

f the analog features of the McIlroy Analyzer. More important, the operating 
taff of only a very few of these calculating boards have run water network 
nalyses. 

The current types of analog computers available for analyses are too 
estricted in pipeline capacity to handle typical complex network problems. 
Intil a small and simple non-linear resistance component is developed for the 
nalog computer, exploitation of this excellent tool will be quite limited in net- 
ork analysis. 

There are at present only 15 McIlroy Analyzers in use, 7 of which are fully 
mployed on gas network analyses, one on air and only 3 on water. Only 4 
machines are readily available to private consulting engineers. 

However, various models of advanced types of digital computers are in 
se in most metropolitan, business and university areas. Because of its 
10re widespread availability, the digital computer is the machine which 
olds the most promising future for use by private consultants. As soon as 

standard set of instructions is published (particularly for the stored- 
rogram type of computer), together with an outline for setting up data for 
roblem mains, loops and special operating instructions, the digital com- 
uter will become an extremely valuable design tool for all consultants. 

jith a set of standard program instructions available, the time required 

f the designer to set up and evaluate a test run should not differ appreciably 
-om the time required for the same test run on a McIlroy Analyzer, pro- 
ided that special instructions (such as balanced pump characteristics) are 
ot too involved. If few runs are made in a given study, or if the designer 
as little experience with the procedure, it will no doubt take more of his 
me to do a test run on the digital computer. A Texas computer firm(11) 
ill perform digital computations via correspondence. This service is a 
tep in the right direction. Waiting time between test runs will not be re- 
iced until either similar services are available across the country or until 
standard program is published for general use. Analysis by means of the 
icIlroy Analyzer has outstanding advantages over analysis by digital com- 
iter, if the design engineer is present during the test run. There is, on the 
her hand, no particular advantage to the design engineer in being present 
iring the performance of iterations by a digital computer. 

Opposed to the private consultant in general practice to whom a network 
lalysis is only an occasional problem, is the distribution designer or planner 
‘a public utility. It is the steadfast opinion of the authors that the McIlroy 
nalyzer is the superior device for distribution design of gas or water 


: 
1588-8 HY 2 April, 1958 


networks when such work constitutes an important part of engineering : 
functions. § 

The analog features of the McIlroy Analyzer are outstanding where dis- 
tribution design is a vital or full-time occupation. The design engineer gets 
the “feel” of the whole system as the required demand is approached. Sys- 
tem revisions can be made in a matter of minutes. Decisions concerning 
revisions or alternate schemes can be made most effectively by immediate 
appraisal of the results as modified or evaluated in terms of observed, 
mentally recorded trends. Proportionate head-loss in mains is indicated 
by the comparative brightness of the resistors. Thus the resistors and flow 
and head meters provide a special auxiliary service as instantaneous, simple 
visual aids. 

With the digital computer, on the other hand, the designer gets only a 
tabulated summary which is meaningless until analysed in toto. The elapsed 
time between conversion of the problem into simple arithmetic and interpre- 
tation of the final results contributes nothing to the problem evaluation. Pos- 
sible revisions or alternate schemes are not readily suggested. Physical 
identity with the problem is lost and the solution is reduced to abstract num-— 
bers in tabular form. 

The Mcllroy Analyzer of the Philadelphia Water Department is used an _ 
average of only one or two full days per week. The remaining work week is — 
devoted to studying and evaluating data for analyses (about 50% of rehabilita- 
tion design work-load), preparing reports, advising on small main local grid ¥ 
changes, and additions consistent with overall distribution planning. The 
Water Distribution Design Unit is keeping up with the Pitometer Company 
C-value surveys and the capital budget program. The machine could be used 
more if the 2-1/5-man staff was increased, but to double machine use it would 
be necessary to double data analysis and the whole operation would become 
quite involved. 


CONCLUDING REMARKS 


To repeat, the accuracy of network analyses, regardless of type of ma- 
chine employed, is of the order of the accuracy of the data fed into the ma- 
chine. The value of the results obtained depend entirely upon their proper 
utilization and interpretation by the design engineer. 

The McIlroy Analyzer was tailor-made for Water Department network 
distribution design and is heartily recommended for use by medium and 
large water departments. The more prevalent digital computers will 
probably become increasingly adapted to network design by private con- 
sultants and will probably be extensively used as soon as general instructions — 
are made available and until such time as there are a much greater number _ 
of McIlroy Analyzers available on a rental basis. If a simple pipe-simulating 
component for analog computers is developed, water-hammer characteristics 
of networks can then be determined with little difficulty and vastly reduced 
approximation. 

As evidence of the progress made in water network analysis procedures, ~ 
one need only consider that as recently as two years ago the question at issue 
was the comparison of “hand” trial-and-error calculations versus utilization 
of the McIlroy Analyzer. However, the gas industry remains well ahead of 
the water field. Six of the monthly issues of Gas, in 1955, gave pertinent 
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discussions related to a series of symposia on the use of various computers 
in gas distribution network piping design (particularly J anuary, May and 
September). 
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SYNOPSIS 


A review of important work performed in the field of sediment densities is 
ven in this paper. The factors affecting the initial density of newly deposited 
diments and the rate of consolidation of deposited material are discussed. 
sults of extensive sampling of deposited sediments are summarized, and 
ethods of estimating sediment densities for particular sediment conditions 
id methods of reservoir operation are presented. Consolidation theories 
veloped in the field of soil mechanics are presented, with a discussion of 
eir possible application to the sediment problem. 


Purpose 


This review is presented as a contribution to the work of the Task Force 
Reservoir Sedimentation, Hydraulic Division, ASCE, which is making a 
view of work in the field of reservoir sedimentation. 

The purpose of this paper is to present a review and summary of the work 
at has been done in the determination of densities of deposited sediments. 
me of this work has been presented as incidental to sediment studies of 
ore general nature and, therefore, has not been given the attention it de- 
rves. Other work has been presented in Government publications that are 
t widely available to the non-Government engineer. Still other work is 
corded in literature of other fields of engineering which is not widely read 


i i i ing date one month 
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Presented at ASCE meeting, New York, N. Y., October, 1957. 
Hydr. Engr., Harza Eng. Co., Chicago, Il. , 
Hydr. Engr., Bureau of Reclamation, U. S. Dept of the Interior, Denver, 
Colo. 
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by engineers working in the field of sedimentation. No original contributions 
are made in this paper. 4 


* 


The Problem { 
: 


The determination of the volume of sediment expected to accumulate in a 
reservoir is a problem which occurs in the planning of all reservoirs. Where’ 
the expected sediment volume is large, it becomes necessary to set aside 
reasonable amount of space for sediment storage, so that the purpose for ~ 
which the reservoir is intended can be served without undue shortages during) 
the useful economic life of the reservoir. Even when the sediment volume 
anticipated is so small that it constitutes no problem it is still necessary 
most cases to demonstrate proof, through studies or verification by me 
ments. 

Estimates of sediment accumulation in reservoirs consist of three pha: 


(1) Estimate of sediment inflow to the reservoir 
(2) Estimate of sediment inflow that is retained in the reservoir 
(3) Conversion of retained sediment from a basis of weight to volume 


Much attention has been devoted to the first phase, the estimate of sedi- © 
ment inflow. Intensive research has been performed on development of i 
proper instrumentation for measuring sediment load in the field, and for 
analysing the sediment characteristics in the laboratory. Much effort has « 
been expended in the routine observation of a network of sediment sampli Ca 
stations throughout the country and in processing the data obtained for gen- 
eral use. Many studies have been made of the proper method of interpre 
these data to arrive at an estimate of the sediment inflow. The frontal attach 
made on this phase has been very essential, very effective, and has justified 
its cost. 

The second phase involves an estimate of the trap efficiency of the reser 
voir. This phase was covered in another paper sponsored by the Task Foret 
The third phase is the subject of this paper. 

The sediment inflow that is retained in the reservoir is determined in 
terms of weight, since all field observations are made on the basis of sedi 
ment concentration by weight. To complete the estimate of sediment accum 
lation in a reservoir, there remains the very short, but highly important, st 
of converting the weight of the retained sediment inflow into volume. This: 
quires an estimate of the space occupied by a unit weight of sediment, which 
is variously referred to as specific weight, unit weight, or density. In this 
paper, the term “density” is used to define this conversion factor. Density 
is expressed as the weight of solid material per cubic foot of mixture of 
solids and water. 

If the density of the deposited sediment is not correctly estimated, the 
timated volume of deposited sediment will be in error by the same percenta 
as the error in the estimated density. The density of deposited sediment in 
reservoirs has been observed to vary from 18 to 125 lbs. per cubic foot. 
Therefore, it can be seen that the determination of a proper density value 
is quite important. Unfortunately, while many observations have been made 
of density under various conditions, relatively little concentrated study has 
been devoted to relating density to the various casual factors. 
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Equipment for Measuring Sediment Density 


This discussion of equipment is not intended to cover all the types of 
nplers that have been developed for collection of samples of deposited 
iment. A description of a number of mechanical samplers is given by 

|. S. Government interagency publication.(1) Two pipe samplers developed 
the U. S. Bureau of Reclamation are typical of the mechanical samplers. 
The first pipe sampler was developed about 1947 and used in a number 
reservoir surveys through 1949. It consisted of a one and one-half-inch 
meter pipe about eight feet long with a check valve near the upper end and 
eaded at the bottom end for attaching sample containers, which were about 
ht inches in length. A rope was attached to the top of the sampler for 
fering it through the water into the sediment at the reservoir bottom. 

The second sampler was designed and fabricated in 1950. It consisted of 
tandard two-inch diameter pipe with a removable cutter head at the bottom 
soil penetration, and a check valve and removable lead collar weights at 

top. Either a five-foot or ten-foot pipe section may be used. A plastic 

e inserted in the pipe sampler is used to secure the sample. This sampler 
jropped over the side of the boat and allowed to fall to the reservoir bot- 

a 

The latest instrument developed by the Bureau of Reclamation is the 
lioisotope densitometer, which was developed and tested under the sponsor- 
p of the Sedimentation Subcommittee of the Federal Interagency Committee 
Water Resources.(2) It consists of a probe containing cobalt 60 as a source 
yamma rays and a self-ionization chamber for the detector. A lead cylin- 
'is provided for radiation shielding and weight to drive the probe into the 
liment. To utilize the phenomenon of Compton scattering of gamma rays 
sediment density measurement, the radiation scattered from the source 
yamma rays back to a detector is measured. When the probe is placed in 
sediment, the intensity of the radiation scattered by the surrounding 

dia to the detector is a function of the electron density of the saturated 
liment. Preliminary field tests of this instrument have shown that the 

ic method of measuring sediment density with gamma ray scattering is 

nd and practical. Further tests of the instrument are contemplated. 


Factors Influencing Sediment Density 


It can be assumed, with reasonable accuracy, that all sediment particles 
e the same specific gravity, about 2.65. The density of deposited sedi- 
nt thus depends upon the amount of void space occupied by water. This 
nfluenced by the size of material, degree of sorting and gradation, and by 
amount of consolidation that has oceurred. 

The rate of consolidation in itself varies with a number of factors. Engi- 
rs in the field of soil mechanics have developed well proven theories for 
solidation of continuously submerged materials. However, engineers in 
field of sedimentation have not yet solved the problems of compaction of 
ler submerged or exposed sediments. 


Definition of Sizes 


nerous classifications have been used to describe the various size ranges 
sediments. Except where otherwise specified, the classification of the 
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U. S. Bureau of Soils is used in this paper. These are as follows: i 
Sand - particles with sizes between 0.05 and 1.0 mm 
Silt - particles with sizes between 0.005 and 0.05 mm 


Clay - particles with sizes less than 0.005 mm : 


i 


Variation in Observed Densities ; 


Many samples of deposited sediments have been collected for purposes ¢ 
determining the density of the deposits. These samples have covered a wide | 
range of sediment characteristics and conditions of deposition. Comprehen= 

sive summaries of these data are given by Lane and Koelzer(3) and by ri 
Hembree, Colby, Swenson, and Davis.(4) The first named summary presents | 
data on over 600 samples collected at nearly 100 locations by various inve 
gators throughout the world. Hembree, Colby, Swenson, and Davis gave d 
on about 225 samples from 38 locations in the United States. The variation — 

that exists in sediment densities is demonstrated by the data in Table 1, ' 
selected from these summaries. These data indicate a general relationshig 
between particle size and density. However, it can be seen from the com- 


parative data on similar size classifications that factors other than size in 
fluence the density. 


Consolidation Theories 


In the field of soil mechanics, consolidation theories for submerged 
materials have been well developed, starting with the pioneering work of 
Terzaghi. Consolidation of submerged material takes place primarily be- — 
cause of imposed load. In the consolidation process, water must be squeezed 
from the deposits; therefore, the permeability of the deposits is a significa ni 
factor. Thus permeable sands have a higher rate of initial consolidation tha 
impermeable clays. 

Terzaghi first advanced the theories under which consolidation under 
constant load varies with time. Figure 1 shows the results of typical 
laboratory tests, carried out at Swarthmore College, which demonstrates 
Terzaghi’s theories. Many similar curves are given in standard references 
in soil mechanics literature. The percent consolidation is the amount of | 
compression at any time divided by the total primary compression. Primary 
compression is that which results from squeezing out the excess water from 
the voids as the load is being transferred to the soil particles. The time of 
primary settlement varies from a few seconds for coarse material to years 
for very fine material. Secondary compression starts when all the load is 
taken by the soil particles. The rate of secondary consolidation is extreme! 
low, with settlement being due primarily to readjustment of the particle 
structure. 

The relationship of imposed load to consolidation has been developed 
through years by many engineers working in the field of soil mechanics. 
Figure 2 demonstrates the relationship between pressure and voids ratio, © 
for clay and sand sizes. The voids ratio is the volume of voids divided by 
the volume of soil particles and is thus inversely proportional to density, 
assuming a given specific gravity. A greater degree of consolidation occur 
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with clay than with sand, but clay sizes maintain a higher percentage of vo 
than the coarser sand sizes. These standard curves developed by soils me= 
chanics engineers appear to be significant to the sediment engineers. Con 
clusions similar to those reached from these curves have been reached by 
Parker Trask, in the field of geologic sedimentation, and by Russel and 
Dickey, in investigating petroleum shales. 

Trask(5) states that compaction of sediments, with resulting decreases 
in porosity, is clearly related to the imposed load. For sediments to com- 
pact, the water must be squeezed from the deposits; therefore, the permea~ — 
bility is also a big factor. Thus, permeable sands have a high initial rate 
compaction compared with impermeable clays. All sizes compact quickly at : 
first, and then the rate of compaction declines. Trask concludes that the es-_ 
sential factors in compaction are load, grain size, permeability, and time. 

Russel and Dickey(6) discuss the compaction of fine silts and clays, sta’ 
that these sizes, when settled out of still water, have porosities greater than ! 
75 percent, the porosity increasing with decreasing particle size. The fine — 
character of the particles causes them to settle into a cellular structure. — 
As the pressure of the overlying sediments increases with continued deposi- 
tion, compaction occurs through collapsing of the cells and expulsion of the t | 
water. Soft clays are stated to have high porosities, but very low permeabili- 
ties, so that the water is expelled slowly. They show curves which indicate 
an increase of density with depth of deposits, in petroleum shales from 600 
feet to 5,000 feet of depth. 

Happ(7) states that, in alluvial sediments, there is no doubt that air drying 
is the major factor in determining density of unsubmerged sediments, at least 
of the finer sediments. He cited cases of two samples of similar fine 
material, one collected a day or so after exposure and the other after at 
least four months of exposure. Density of the first sample was 48.3 pounds — 
per cubic foot, and that of the latter sample was 94.2 pounds per cubic foot 

From the above, it is apparent that the primary factors influencing the 
rate of compaction of particles of a given size are, (1) the weight of overlyin 
sediment, (2) the degree of exposure to drying, (3) particle size, (4) permea- 
bility (probably related to particle size), and (5) time. 

Some of the foregoing theories are verified by the results of the compre- 
hensive sedimentation survey of Lake Mead, carried out in 1948-1949.(8) In 
this survey, the deposited sediments were intensively sampled. The Lake 
Mead deltas are geographically graded, the coarsest particles accumulating 
at the head of the deltas, with the finest particles being carried further out 
into the lake. As the delta has built downstream, new beds of coarse mater! 
have deposited over old beds of fine materials in some areas. The sand de 
posits, although almost continuously submerged, are firmly packed. On the 
other hand, the silt and clay deposits have an extremely high water content 
which gives them the consistency of a soft ooze. Only the samples of the 
most deeply buried layers are solid enough that they do not flow of their owt 
weight. There are also numerous gas cavities in some zones of the finer _ 
deposits. The water content of the new finer deposits is commonly more 
than 80 percent. As additional material is deposited, part of the water is — 
driven out. Water content in a vertical section has been observed to vary — 
from 16 percent in the sand at the top to 75 percent in the silt and clay of 
the bottom. The water content in all sections decreases with increasing 
depth. The density has been observed to vary from 18 to 108 pounds per 
cubic foot. 
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Initial Density 


The sedimentation engineer is interested in the variation in density 
hroughout the economic life of the reservoir. He, therefore, wants to know 
he initial density, the density at the end of a given period, say 100 years, 
nd the average density throughout the period. From his viewpoint the 
lensity at the end of any reasonable short time, say up to a year, can be 
onsidered initial density. Since primary consolidation of large sized 
naterial has been demonstrated to occur very rapidly, with only minor 
econdary consolidation occurring thereafter, the data from practically all 
amples of coarse material can be assumed to represent initial density, 
egardless of the time the material has been deposited. This will not be 
rue of the finer materials. Initial densities of sediments, as determined by 
arious investigators, are shown in Table 2. 

The values of density developed by the California State Division of Water 
tesources(9) have gained considerable acceptance for coarse sediments. 
these values were developed after a very thorough study of the voids and 
yater retaining capacity of sediments, in which more than 200 samples were 
nalyzed. It should be noted that the sizes are identified by the maximum 10 
ercent size range, whereas all other data in the table are in terms of the 
omplete size range or the median size. 

Trask(10) separated sediments into various size ranges and determined 
heir density after what he termed a “very short settling time.” The data 
f Hembree, Colby, Swenson, and Davis are from samples collected near the 
urface of submerged sediment deposits, which were believed to be formed 
fithin a few years of sampling time. 

Data for 65 samples taken of alluvial deposits in the Middle Rio Grande 
‘alley are recorded by Happ. Although he found great variation for similar 
izes, he found a general relationship in which density increased with median 
article size. He attributed most of the variation to differences in compac- 
ion, and drew envelope curves which he believed represented the initial and 
ompacted densities. His values in Table 2 are taken from an envelope curve 
f minimum densities which represent what he believed to be initial densities. 

A comparison of the laboratory values determined by Trask with the field 
amples determined by the other investigators indicates higher densities for 
ne field samples. This could be due to some consolidation occurring prior 
) the field sampling, even though the deposits were assumed to be relatively 
ew and unconsolidated. 


Variation of Density with Depth of Deposit 


The variation of density of marl sediments with depth of deposit was 
etermined in Lake Niedersonthofen, Bavaria, By Reisinger.(11) These ob- 
ervations reflect the combined effect of consolidation due to time and the 
eight of overlying sediment. The results are shown in Table 3. 


Variation of Density with Drying 


Happ states that his envelope curves are representatives of the approxi - 
late magnitude of the factor of compaction by drying in the Middle Rio 
rande. His results are shown in Table 4. The higher compaction shown 
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IN LAKE NIEDERSONTHOFEN, BAVARIA 


Depth Density 
(meters) (lbs. per cu. ft.) 

Upper Layer 21.6 
1 39.6 

Ps 41.5 

3 34.8 

+ 4201 

5 45.1 

6 52.6 

7 53.9 

8 56.0 

9 73.6 
10 74.0 
ll 80.4 
ie 85.2 
rT Dia it 
14 OZ ent 
15 84,0 
16 Sie2 
17 85.3 
18 O2ac 
19 88.1 
20 89.6 


"ABLE 4~ ESTIMATED COMPACTION OF MIDDLE RIO GRANDE SAMPLES 


(BY HAPP) 
Median 
Particle Size Initial Density Terminal Density Ratio of Compacted 
(mm) (lbs. per cu, ft. ) (lbs. per cu, ft.) to Initial Density 
ak 83 104 L330 
05 78 102 P33 
Ol 72 98 36 
005 68 ot 1.43 
003 59 5 1. 61 


20012 48 94 1,98 
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for smaller sizes bears out the consolidation curves shown on Figure 2. 


Consolidation of Sediment with Time 


Lane and Koelzer studied the sediment density data available and develo 
a formula relating density to particle size, time, and method of reservoir 
operation. Their classification of reservoir operations was as follows: 


(a) Sediment always or nearly always submerged - includes reservoirs 
normally kept full and emptied only for a short time, at infrequent intervals. 


(b) Normally a moderate reservoir drawdown - includes cases where 
reservoirs store over long periods, with small drawdown in normal years. 
It also includes run-of-river hydroelectric plants, where the basin is used 
for pondage. 


(c) Normally considerable drawdown - includes reservoirs with relativel 
large drawdown in normal years. 


(d) Reservoir normally empty - includes retarding basins and reservoirs 
that are empty in normal years. 


The Lane and Koelzer formula is shown in Table 5. Values of Wy were 
determined by examination of over 600 samples, giving consideration to 
degree of drying considered compatible with the method of reservoir opera- 
tion. Values of K for the first three methods of operation were chosen to 
give the same density over a 1,000 year period as would be reached under 
condition (d) in one year, where the reservoir is normally empty and the 
sediment would have ample opportunity to dry out. The logarithmic correc- 
tion for compaction was used to give a rapid consolidation in the first few 
years, with rapidly decreasing rates later. This follows the general trend 
found by soils mechanics engineers, although it seems probable that equa- 
tions could now be developed which more closely approximate the primary 
and secondary settlement rates shown on Figure 1. For material covering 
a range of particle sizes, it was recommended that the equation be used with 
appropriate weights applied in proportion to the percent by weight in each 
size classification. An illustration of this computation is shown on Table 5. 

Miller(12) of the U. S. Bureau of Reclamation has made studies which he 
suggests indicate that the initial densities in the above formula are too high. 
He believes that the values determined by Trask, as shown in Table 2, are 
more applicable, particular for the finer sizes. Trask’s values are for the 
always submerged condition. It should be noted that Trask’s values are for 
density after “a very short settling time,” in contrast to density values of 
Lane and Koelzer which are represented as being after one year. 

The Bureau of Reclamation uses Miller’s modification, considering that 
the “always submerged” condition is more representative of initial densities 
and that Trask’s values for clay sizes are more applicable. 


Methods of Computing Average Sediment Density 


The data presented lead to computation of sediment density by applying 
specific densities to different size grades. Miller presented one method and 
Hembree, Colby, Swenson, and Davis presented two methods of applying thes 
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TABLE 5 
LANE AND KOELZER FORMULA 


W=W)+ K lug1oT 


where W =density after T years in lbs. per cu. ft. 
W) = density after 1 year in lbs. per cu. ft. 
K =a constant for each sediment class and operation 
condition to reflect consolidation 
T =number of years of consolidation 


Values of W and K 


SS ee ee 
We 


d Silt 
Reservoir Operation w 


a. Sediment always submerged 
or nearly submerged 93 


b. Normally a moderate 
reservoir drawdown 93 


c. Normally considerable 
reservoir drawdown 93 


d. Reservoir normally empty 95 82 0.0 | 78 


Example of computation of density of deposit "always submerged"! 


1603-13 


Deposit consists of sand = 20 percent 
silt = 40 percent 
clay = 40 percent 


Density of deposit after 100 years consolidation ~ 0.20 (93) 


+ 0.40 (05+5.7 log)g 100) + 0.40 (39 + 16 logjg 100) = 74.0 lbs. /cu-ft. 
ESS ESET 
! MILLER EQUATION 


| Average density of all deposits after T years of compaction = 


T 
Wave = Ww) + 0,434 K aay (loge T) ol 


: 
j 
| 
4 
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data to suspended sediment data to arrive at an overall value representing 
the average for the suspended sediment load. 

The Lane-Koelzer equation gives the density of the first years deposition 
after T years of consolidation. With sediment continuing to enter the reser= 
voir, the average density of all sediment deposited would be the average of 
all the densities from one to T years of consolidation. Miller integrated th 
equation from one to T years and developed the equation shown at the botto. 
of Table 5. Waye is the approximate density for the sediment deposits after 
T years of compaction. 

In the first method by Hembree, Colby, Swenson, and Davis, median pe ) 
cle sizes for each sample were plotted against the suspended sediment load. 
A curve was drawn through all points on the curve so that about 50 percent 
by weight of the particles in each size range would be finer than the median 
size indicated on the curve. Median particle sizes were taken for predeter 
mined size ranges and combined with their values of density given in Table 
to obtain an average density for the entire sample. 

The second method of Hembree, Colby, Swenson, and Davis was to divide 
the sediment into size fractions and compute the density from the percentage 
in each size fraction and the density of each size. To determine the averag 
percentage of sand, silt, and clay in the suspended sidement, correlation 
curves between percentages of each class and suspended sediment discharge ~ 
were prepared. Frequency tables of suspended sediment discharge were ap= 
plied to the correlated values to determine the tonnage in each classification 
and the average percentage of each classification thus computed. 


a a a a 


CONCLUSIONS 


It is concluded that adequate data are available for determination of densi- 
ties of coarse sediments, in which consolidation is not a problem. However, — 
it is the opinion of the writers that the problem of estimating sediment densi- 
ties for finer sediments has not yet been satisfactorily solved. The Lane- 
Koelzer method, with the modification by Miller, is a rational approach, but 
does not appear to utilize all the techniques of the soil mechanics engineer ~ 
that are applicable. It appears that the “always submerged” condition could 
be readily solved by application of the consolidation theories of solid me- 
chanics. While this does not represent the only type of sediment deposition, _ 
it is certainly one of the more important. The effect of drying on exposed or 
partially submerged deposits needs further study before such techniques can | 
be applied. 
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OCEAN WAVE FORCES ON CIRCULAR CYLINDRICAL PILES® ( 


Closure by R. L. Wiegel, K. E. Beebe, and James Moon 


R. L. WIEGEL,! K. E. BEEBE,” and JAMES MOON.3—Anyone who has 
salt with large ocean waves would not be surprised at a large amount of 
satter in measurements of phenomena associated with these waves. No two 
-ean waves are identical, they only appear to be so due to the rather crude | 
ethods used to describe wave characteristics, considering that they are 
mn-uniform in three dimensions. This is especially true when both swell 
id locally generated wind waves are present. 

Mr. Stelson makes the statement that there was an area of deficient under- 
anding regarding inertia drag. There is much that is not known on inertia 
"ag, especially when associated with large velocities and the resulting phe- 
ymena; however this is not what he refers to. Mr. Stelson has either mis- 
iderstood the definition of the term Cy, which he has called the added mass 
ynstant—which it isn’t—, or he has misunderstood the problem in general. 

e are dealing with a fluid which accelerates about a cylinder at rest, not 
ith a cylinder accelerating through a fluid at rest—which is what he dis- 
isses. The force necessary to accelerate a mass through a fluid is associ- 
ed with the mass itself and what is termed the added mass (for a right cir- 
ilar cylinder the added mass constant is unity). The inertia force exerted 
ra fluid accelerating about a submerged body is associated with the mass 
the displaced fluid and the added mass; the coefficient of mass as given by 
amb for this case is Cy = 1+ 1 = 2 fora right circular cylinder.* This 
velopment of the coefficient of mass was used because of the ease in under- 
anding the physical phenomenon. Of far greater importance is the fact that 
acCamy and Fuchs(1) have solved the diffraction problem for a right circu- 
r cylinder with a diameter small compared with the wave length, using the 
near theory for waves in shallow water. When the resulting equation is put 
the writers’ terms it is 


Fay =. 2 eV ou 


1ich shows a value of 2 for what has been called the coefficient of mass. 


Proc. Paper 1199, April, 1957, by R. L. Wiegel, K. E. Beebe, and James 
~Moon. 

Associate Research Engr., Inst. of Eng. Research, Univ. of California, 
Berkeley, Calif. 

Graduate Research Engr., Inst. of Eng. Research, Univ. of California, 
Berkeley, Calif. 

Cons. Engr., Signal Oil & Gas Co., Los Angeles, Calif. 

See Reference 14 of original paper. 
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The scatter cannot be accounted for any one of the reasons stated by the 
writers; however, each of these reasons accounts for a portion of the scatte! 
Undoubtedly there are also other reasons. 

Some advances are being made in understanding the basic phenomenon. 
For example, data in the paper presented by John McNown and G. H. 
Keulegan, “Vortex Shedding and Resistance in Unsteady Flow” at the Octobe 
1957 meeting of ASCE showed a strong dependency of the coefficients of m 
upon the wake for a right circular cylinder in standing water gravity waves. 
Further, tests made by A. D. K. Laird at the University of California using 
right circular cylinder pair extending through an air-water free surface as 
the mass of a long pendulum oscillating in a direction approximately paralle 
to the free surface showed that variations in forces as great as 50 to 100 
percent could occur when large eddies in the wake of one cylinder hit the 
other cylinder in a certain manner. 
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THE ESTIMATION OF THE FREQUENCY OF RARE FLOODS? 


Closure by Benjamin A. Whisler and Charles J. Smith 


BENJAMIN A. WHISLER, | M. ASCE, and CHARLES J. SMITH,2 A.M. 
SCE.—The writers welcome the comments of Messrs. Beard, Benson, and 
filliams in their discussions of this paper. Several questions have been 
aised which the writers feel will be answered only after more data become 
vailable and after the suggested procedure and other procedures have been 
pplied to the data and the results compared. Since there appears to be no 
eneral approach to the problem of estimating the frequency of rare floods 
lat proves to be satisfactory under all conditions, the writers believe that 
ome completely new thinking, at present unknown to them, is needed. In 
me, with enough people working on the problem, a satisfactory approach 
lay be found. 


Proc. Paper 1200, April, 1957, by Benjamin A. Whisler and Charles J. 
Smith. 
Head, Dept. of Civ. Eng., Pennsylvania State Univ., Univ. Park, Pa. 
Associate Prof., Dept. of Civ. Eng., Pennsylvania State Univ., Univ. Park, 
Pa. 
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SYSTEMATIC CHANGES IN BEDS OF ALLUVIAL RIVERS@ 


Discussion by Vito A. Vanoni and Hsin-Kuan Liu 


VITO A. VANONI,! M. ASCE.—The authors are to be commended for 
winging their important observations of bed configuration of the Mississippi 
tiver to the attention of the profession. Such detailed observations, tedious 
nd time-consuming as they may be, are necessary to the advance of our un- 
lerstanding of the mechanics of alluvial streams. In this regard, the writer 
ot only agrees with the authors conclusion(1) that the study of river bed be- 
avior might throw new light on some theoretical aspects of river hydraulics, 
ut he feels certain that it will do so. 

The authors’ belief that sand waves in alluvial rivers form a major source 
f resistance to the flow is borne out by laboratory studies, as shown by the 
ata(A) in the following table. These data were obtained from experiments in 
yhich the bed of the flume was covered with loose sand with a geometric 
nean grain size of 0.152 mm and geometric standard deviation of sizes of 
-76. The flow rate was set for each experiment and the slope was varied to 
btain uniform flow. It will be seen that the friction factors with dune- 
overed beds are from two to five times as much as those for the flat-bed 
ondition when no dunes exist on the bed. It will also be noted that as the 
elocity and the sediment-discharge concentration increase the friction factor 
ecreases. Detailed observations of all dunes were not made in these experi- 
nents; however pictures and visual observations did not show any great 
hange in size or shape of dunes until the sand wave and smooth phases were 
pproached, i.e., runs 6, 7 and 8. 

The behavior of dunes in the flume seems to be opposite to that observed 
y the authors in the Mississippi River where dunes increase in size with 
tage. However, in the flume the increase in velocity was achieved by in- 
reasing the pump speed and changing the slope, while in the natural river 
1e slope remains essentially constant as the stage varies. Therefore, the 
WO cases, i.e., river and flume, are not exactly the same. Any quantitative 
ata that the authors have on the size of waves at high stages would be of 
onsiderable interest. 

It would be interesting to obtain simultaneous measurements of friction 
ictor and dune configuration for the river, and it is hoped that the authors 
ill have an opportunity to make such observations. It is also hoped that such 
bservations can be extended to the highest flood stages. In the absence of 
uch data one can only speculate on the effect of dune size and configuration 
n the friction factor. One is inclined to suppose that the friction factor 


, Proc. Paper 1331, August, 1957, by Walter C. Carey and M. Dean Keller. 
. Prof. of Hydraulics, California Inst. of Technology, Pasadena, Calif. 
Letters refer to bibliography at end of discussion. 
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would increase as the dune size increases. Since the highest dunes were Ob- 
served at the highest stages, this would mean that the friction factor in- i 
creases as the stage increases. However, measurements of the Manning 
roughness coefficient n, by Eden(B) for the Mississippi River in the vicinity 
of Memphis, Tennessee, showed that as the stage increased beyond about 30 ; ) 
feet, n stayed approximately constant or decreased slightly. As the stage | 
dropped below about 30 feet, Eden found that n increased rather rapidly and 4 
attained the highest value at the lowest stages. Eden also observed that the j 
largest dunes occurred at the highest stages. He observed dunes as high as 
34 feet, which agrees well with the maximum height reported by the authors 
for stages well below the flood stage. 

A careful observation of the profiles presented by the authors seems to 
show that more of the bottom was covered by small dunes at low stages than 
at high stages. It is possible that the small dunes actually offer more re- 
sistance to flow than the larger waves, which seem to have relatively smooth- 
er surfaces than the configurations of smaller height which occur at the lowe 
stages. 

The authors explain the loop in the rating curve in part by postulating a 
lag in adjustment of the bed configuration. It seems unlikely to the writer 
that such a lag exists since the rise and fall of the flood occur slowly, re- 
quiring days and weeks to take place. During this time the stages are rela- 
tively high, as are the velocity and sediment transporting capacity of the 
stream, and it seems to the writer that the stream should be able to wash 
away the larger dunes and readjust the bed with very little lag. The fact that 
the authors find that crossings adjust very rapidly lends support to the writ- if 
er’s belief that the lag is small. 1 

By following the idea advanced by Brooks, (C) that both flow rate and sedi- { 
ment transport rate are independent variables, one can explain the loop in 
the rating curve without resort to the idea of lag advanced by the authors. 

In line with the idea proposed by Brooks, one visualizes the dependent vari- 
ables, such as velocity and friction factor, as adjusting themselves to the in- 
dependent variables of flow rate and sediment transport rate. For instance, 
let us compare conditions in the stream for the same flow rate during the _ 
rising and falling stages of the flood. On the rising stage the sediment load 
is usually high, (D) so the stream must adjust itself to transport this load. 
This is done by smoothing out the bed and achieving a relatively high velocity, 
On the falling stage, the sediment load is usually smaller than on the rising 
stage, so the stream can transport the load with a lower velocity and larger 
depth. The stream then adjusts to this condition by forming rougher dunes 
than it did on the rising stage. 


a 
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Darcy-Weisbach Friction Factor f, for Sediment-Laden Flows 


0.241 ft Deep ina Flume 10.5 in. Wide 


U Cc f 
Run Average Sediment Friction Bed 
No. Velocity Discharge Factor Condition 
fps Conc. 
gr/L 
2b 0. 80: 0.30 0.124 Dunes 
2a 0.85 0.59 OFS Dunes 
2 0.91 0.82 OFS Dunes 
3a 0.98 aS 0.108 Dunes 
3 1.04 r.8 0-103 Dunes 
+ rT. 20 Faas: 0.075 Dunes 
5 1.38 3.4 0.050 Dunes 
6 ji EE z.9 0.038 Sand wave 
6 1.69 3.3 0.029 Sand wave 
8 i, 83 £9 7 0.024 Flat 
1 2.06 3.4 (Oral Flat 
S) 2.66 5.6 0.022 Flat 


HSIN-KUAN LIU* A.M. ASCE.—The authors are to be congratulated for 
resenting information on the sand waves of a natural river. Such informa- 
on is highly valuable to the profession. 

For several years, the writer has studied sand-wave motion under labora- 
ry conditions. The data collected from laboratory experiments are general- 
‘for steady, uniform flow in narrow straight channels. It is obvious that 
lis type of flow is different from what has been observed by the authors 
om the Mississippi River. Ordinarily the flow phenomenon of natural 
ivers is so complex that it is not possible to understand the hydraulics of 
vers based upon field investigations alone. On the other hand, conclusions 
ised upon investigation under simplified laboratory flow conditions should 
9t be applied to prototype flow conditions without due consideration of the 
ymplex problems involved in natural rivers. In the laboratory it is easier 
study analytically the primary forces which affect the movement of 


Asst. Prof. of Civ. Eng., Dept., Colorado State Univ., Fort Collins, Colo. 
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sediment. As a result, various concepts of mechanics of sediment transport — 
have been formulated. It is known that these primary forces such as tractive — 
force, lifting force, turbulent stress, vary with the flow conditions. The fl 
conditions in natural rivers are somewhat different from those in the labora~ — 
tories. For example, in natural rivers the flow is non-uniform and unsteady, — 
the scale of turbulence is large, and the geometric scale of a river is large; 
these are not so in the labroatory. Despite the fact of these differences, the 
authors’ conclusion based upon field investigations are shared to a large ex- 
tent by the writer who has been acquainted only with laboratory conditions. 

The authors are correct in stating that sediment transport in the bend re-— 
gion is different from that of a'straight reach. In the past, most of the labora- 
tory study of sediment transport has been conducted in straight channels. | 
Laboratory study of sediment transport in curved channels has not been don 
extensively. In view of the fact that river bends play an important role int 
river hydraulics, it is obvious that studies of sediment transport around the ~ 
bend should be emphasized both in the field and in the laboratory. ¢ 

The statement that sand waves occur generally on the bed of a river in a } 


systematic manner agrees with the laboratory data. This fact does not sup- } 
od 


port the hypothesis that sand waves and ripples are caused by turbulence as 
suggested by some hydraulic engineers. Needless to say that turbulence 
inter-reacts with the sand waves which in themselves cause turbulence. 

That the size of sand waves becomes larger with increases in discharge 
and smaller with decreases in discharge may be owing partly to the effect of 
change of depth, partly to the effect of change of velocity.(a) Studies(b) madell 
by W. H. Powless indicated that the sizes of the sand waves are also affected — 
by the rate of rise and fall of the stage. 4 

The writer agrees with the authors that ‘the major sources of resistance 
arises from systematic river-bed changes, which the river makes of its own 
volition, and modifies to meet changing discharge requirements.” In other 
words, the bed roughness may have correlation with the bed configuration. 
However, the bed configuration is not an independent variable, but it is de- 
pendent upon the river itself. 

No information regarding the rate of the sand-wave motion was given by 
the authors in their paper. According to laboratory data, the writer found 
that 


Gory (a) 


in which C is the velocity of the sand waves in ft per sec. 
V is the mean velocity of the flow in ft per sec 
N is the pure number 
K is an empirical coefficient, probably is a function of the flow and 
the bed material. 


According to laboratory data, the N-value is 5 and K is in the order of M X 
10-4, in which M varies from 1 to 10. According to data collected by Wolmar 
on the Du Noir Creek, Dubois, Wyoming, equation (a) still holds true as 

shown in Fig. A and N-value is also 5. For natural rivers there is very 

limited information regarding the rate of sand-wave motion. Table I was 
prepared by the writer by using data on the Mississippi at Fulton as reported 
by the Mississippi River Commission. 
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Table I. Rate of Travel of Sand Waves 
rate of travel 

pss ft/day ft/sec 

t 

22 to 24 rising 

18 to 22 rising 

13 to 18 rising 

16 to 15 falling 

15 to 13 falling 

13 to 12 falling 


Mean velocity 


Mate of Surveys of flow 


uly 7 to July 13 
uly 3 to July 7 
‘une 23 to July 3 
lune 11 to June 15 
lune 15 to June 19 
‘une 19 to June 23 


Jata in Table I are plotted on Fig. A. Those of falling stage are plotted 
eparately from those of rising stage. Data are insufficient for drawing 
lefinite conclusions. However, there are several points that can be dis- 
ussed, 


a) The K-value for the sand waves found in the Mississippi River is much 
maller than that of laboratory data. 


b) There are not enough data to determine the N-value. 


c) The mean velocity on a rising stage needs to be higher than that on a 
alling stage for a constant rate of sand wave motion. Regarding items c, W. 
i. Powless stated in his report(b) that the rate of sand wave motion is much 
reater on a rising than on a falling stage. Additional information of sand 
yave motion in Mississippi River is reported by E. W. Lane(C) and E. W. 
iden. (d) 

The Mississippi River data seem to indicate that the sand waves are in- 
luenced by the rising and falling stages of the flows. In order to understand 
his phenomenon, a laboratory study of sediment transport by introducing a 
ydrograph into the flow system will be very helpful. 

The writer had the opportunity to discuss the authors’ paper with Mr. E. 
¥. Lane. Mr. Lane thought that the authors’ information is highly valuable, 
owever, he pointed out that most of the data presented by the authors were 
ollected in a reach where the slope of the Mississippi River is appreciably 
ffected by the sea level. Furthermore, the Mississippi River in this reach 
3 deep and narrow. In order to understand fully the sand wave motion in the 
fississippi River data should be collected also in other reaches where the 
.Ow conditions are different. 

In conclusion it can be stated that, alluvial rivers, such as the Mississippi, 
arry not only large amount of water but also large amount of sediment. In 
rder to understand the nature of these rivers, investigation and study have 
9 be made not only of the motion of the water, but also of the motion of the 
ediment carried by the river. To a certain extent the study of the surface 
ater is easier because it can be seen from above the water surface, and ex- 
snsive studies have been made in an effort to understand its hydraulics. The 
tudy of sediment transported by an alluvial river is obscured by the fact that 
e movement is under water. For this reason, special effort has to be made 
1 order to understand the behavior of alluvial streams. Such understanding 
; essential in order to control and develop alluvial rivers for the welfare of 
jankind. 

The writer wishes to express his appreciation toward Messrs. M. G. 
folmanand E. W. Lane for their help in the collection of data and of refer - 
ices, respectively. 
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TURBULENCE IN A DIFFUSOR BOUNDARY LAYER2 
Discussion by Theodor S. Strelkoff 


THEODOR S. STRELKOFF.*—Before the problem of turbulent boundary 
yers can be solved completely, it seems apparent that a wealth of data will 
required, data which relate turbulence characteristics to mean flow 
rameters, and which are obtained from flows confined by various geomet- 
sal configurations. This fact makes the present paper most welcome. 
With reference to the terms in the turbulent energy balance, the produc- 
nm term, when describing conditions in flow with a pressure gradient, es- 
cially anywhere near separation, should include in addition to we at 

ist u' gu a variable which may be half the size of the shear stress term. 
though another applicable term, v' x » partially nullifies the effect of the 


ear stress. (1) 

(The expression for Wy, 4s given in the paper apparently suffered from a 
Sprint as it was being non-dimensionalized into its final form, because in- 
ad of U4 one reads U4.) 

It would be of some interest to know what method was used in placing “the 
ine of the hot wire at the wall of the diffusor.” As it is generally difficult 
locate the probe accurately, any information on this point would be helpful 
other researchers. In addition, one would wish to know what corrections 
any were applied near the wall to compensate for heat losses to the wall 
elf. 

If the true three-dimensional momentum thickness had been plotted, com- 
rison with Uram’s data(2) would still have been convenient, and in addition 
. Karmdn equation, modified to include normal stress terms, could have 

2n employed for calculation of the shear stress at the wall to be compared 
sh the value obtained at station 12 with the pitot tube method. 

It would also have been helpful if a description of the “history theory,” 
erred to in Fig. 6, had been supplied, together with indications of how the 
‘culations were made and what the various symbols stood for. Schubauer 
i Klebanoff,(1) basing their reasoning on a simplified form of the Reynolds 


lation, which neglects viscous stresses, derived the relation as oye 


; ; oR 
ere T is the turbulent shear stress. It is interesting that at station 12 a 


1 be computed as 2.2 psf/ft, in fair agreement with oP = 2.56 psf/ft. The 


Proc. Paper 1393, October, 1957, by J. M. Robertson and G. L. Calehuff. 
Research Associate, Iowa Inst. of Hyd. Res., State Univ. of Iowa, Iowa 
City, Iowa. 
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value of 2.2 psf/ft was obtained from the stated value of Rg = 25, 000, by as- 
suming an air temperature of about 80° and taking the slope of the shear 
stress coefficient function at y = 0 as 0.016 per inch. Actually, this computa- 
tion should be carried out with variables whose values are directly known, as 
an error of 20° in air temperature, for example, would cause a 10% error in 
Or 

oy 

If the graphs of u' had been normalized through division by the velocity Uy 
instead of U, the peak, noted by other researchers studying adverse pres- 
sure gradient boundary layers, would have appeared at each station. The 
data of Schubauer and Klebanoff(1) peaked in this manner; that of Sandborn 
and Slogar(3) exhibited a somewhat similar peak at the station farthest down- 
stream. Both sets of data applied to two-dimensional boundary layers. 

It would have been informative to see a graph similar to figure 8 with data 
from all stations superimposed as an indication of the applicability of the 
inner-law to turbulent fluctuations. 

Further studies with the same apparatus seem to be definitely in order as 
time permits. In particular, careful investigation of U, u', and uv very near 
the wall should prove profitable. 
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100 FREQUENCY CURVES OF NORTH AMERICAN RIVERS2 


Discussion by G. N. Alexander, Max A. Kohler and Ralph W. Powell 


G. N. ALEXANDER* A.M. ASCE.—The first conclusion made by the author 
in his synopsis reads - 


“1. That frequency curves can generally be drawn as straight lines.” In 
the light of the paper this statement is interpreted as meaning that the most 
Suitable distribution for flood flows is the log-normal. The author carried 
out a sampling (Monte Carlo) experiment where he assumed a log-normal 
population of 1000 floods and drew twenty random samples of size fifty. He 
noted that two of them were far from linear and rightly states that this “is no 
conclusive proof that frequency lines have to be straight.” In fact it is just 
the opposite. From this experiment it is clear that samples, even if drawn 
from the ‘true’ population may be far from linear. What is more, even if a 
Sample is linear it may still be a poor representation of the population. Even 
in a record of 1000 years the order statistics will show a pattern of scatter- 
ing similar to a small sample. Presumably, one would obtain an even greater 
scattering than that found by the author if the samples of fifty were drawn 
from a ‘sample’ of 1000 from a very large population, rather than the ‘true’ 
population of 1000. 

Statisticians, and engineers familiar with quality control, have recognized 
these sampling effects, but many hydrologists are still not sufficiently aware 
of the effects of sampling and hence apt to ascribe them to other causes. The 
problem is to differentiate between sampling effects and characteristic ef- 
fects, and the writer believes that is best done if we first have some concept 
of the effects of sampling. 

These effects can be shown by using a ‘band’ instead of a line to represent 
the distribution. The width of the band would approximate to the interval be- 
tween the outer pair of samples obtained by the author from his 20 samples 
if we wish to cover about 95% of the samples from his population of 1000. 
Hence if we were to use a 5% significance test one would ascribe fluctua- 
tions within the band to sampling effects rather than seek natural causes. 

The fourth conclusion in the author’s synopsis reads - 


“4, That an increase in the size of the drainage area tends to flatten the 
frequency line.” 

The author has tried to prove this by an a priori argument comparing the 
slope of the frequency curve of annual floods obtained for stations a or b on 
tributaries identical both physiographically and climatologically with that for 
€ just below their confluence. In the standardized log-normal distribution as 


a. Proc. Paper 1395, October, 1957, by E. Kuiper. 
* Hydr. Engr., U. S. Geological Survey, Washington, D. C. 
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plotted by the author, the slope may be measured by the ratio of the standard 
deviation o to the mean , i.e. the coefficient of variation C.V. 

Now under the above assumptions the mean flood at c, mu (c) say, will be 
twice the mean flood at a, u (a) or 


H (c) = 24 (a). (1) 


If in the first instance we assume the floods at c are a random conbination of 
the annual floods at a and b then for uncorrelated distributions, the variances; 
(square of the standard deviations) may be taken as additive (this ignores 
their non-normality) 


Le. | Vara c)== 2aVareta) 
(2) 
Or ta Cet cera) 


or from (1) & (2) 


Cove) ea WF 
C.V. (a) V2 
Hence, in this case slope of c is flatter than slope of a. 
However, since we have tributary streams there will be some correlation 


p say between their flows and the relation between a, b, & c may be written 
as 


Var(c) = Var(a) + Var(b) + 2p. [Var(a).Var(b)]2 
For identical tributaries 
Var(a) = Var(b), 
Hence 
Var(c) = 2(1+P). Var(a). (3) 
Hence, using (1) and (3), 
oa ation ee el 
giving as before = if P=0. If P = 1., i.e. perfect correlation the slopes are 


identical. 

However, the floods at c will not generally be a combination of the “annual 
floods at a and b as assumed above. We could divide each annual flood at c 
into its a ‘and b “components” and hence obtain the distributions of floods in 
aandb causing the annual floods at c. The distributions of these components 
at a will usually have a wider range ‘and hence greater variance than the an- 
nual floods at a, and will also have a lower mean. Hence the “component” 
distribution at a will have greater C.V. than the “annual” distribution at a. 
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Thus although the combination of the components will result in annual floods 
at c with a lesser C.V. than the component floods in a or b, this C.V. may or 
may not be less than the C.V. of the annual floods at a or b. 

Actual analyses are therefore, required to justify the general thesis that 
the slope decreases with the drainage area: Such analyses now being under- 
taken by the U.S. Geological Survey do support this thesis. 

The author raises the question of suitable parameters—stating “it would 
be desirable if the inclination of the frequency line could be expressed in 
Some numerical parameter.” He has selected the ratio of the 1% flood to the 
mean annual flood. In dealing with the log-normal distribution a decision has 
to be made concerning the stage, if any, in the analysis at which the original 
data will be transformed to log-units. The statistical analysis is consider- 
ably simplified and significance tests depending on normality may be used 
more readily if the transformation to log-units is made at the outset. In this 
case it is logical to use the independent and efficient parameters, namely, 
the mean ? and variance o2 of the common logarithms as the parameters to 
describe the distribution. It is readily shown (1, p. 161) that in a “true” log- 
normal distribution the mean of the logarithms is also the geometric mean 
of the original variate. This value is also the 50% or median value, with a 
mean return period of two years, irrespective of whether we use the loga- 
rithms or the original variate. It therefore seems more logical to measure 
“central tendency” by the median or geometric mean rather than by the 
arithmetic mean of the original observations used by the author. The rela- 
tion between a, Y & a is given by (Ref. 1, p. 161) 


108, 5 % = ct Lisa". 


As a measure of ‘dispersion’ the standard deviation, in some form or 
other, is almost universally used. For some curious reason, Statisticians 
who derided as old fashioned the engineer’s use of the “probable error” with 
its 50% interval were apparently quite happy to measure in percentage sig- 
nificance levels, rather than using the natural unit, viz.; the standard devia- 
tion. As probability paper is not normally printed to facilitate measurement 
in standard deviations, and statistical tables are usually constructed with 
Significant levels shown percentage wise, it is generally more convenient to 
use percentages as the author has done. 


REFERENCE 
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MAX A. KOHLER,* M. ASCE.—The study reported by Mr. Kuiper was di- 
rected toward the development of a technique for using discharge records at 
numerous points in an area when developing the frequency curve for a short- 
record station. This generalized approach to frequency analysis may be ex- 
pected to yield improved results (over single-station analysis), and is perhaps 
to be recommended when the nature of the problem does not justify detailed 
study of rainfall-runoff relations,(1) etc. The approach is in some respects 
similar to that widely employed by the U. S. Geological Survey,(2) but 


* Chief Research Hydrologist, U. S. Weather Bureau, Washington, D. C. 
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Mr. Kuiper attempts to introduce “climate” as a factor, whereas the Survey : 
restricts each regional study to an area of climatic homogeneity. 

Although much is said of climate in the paper, the discussion under “Pos- 
sible Application” would indicate that the frequency curve is derived from 
only two independent factors—mean-annual flood and drainage area. Dividing 
the 100-year flood by the mean-annual and the mean-annual flood by the 
drainage area only tends to obscure the high correlation which exists in the © 
basic data. 

The upper plotting of the accompanying chart (derived from data in Table — 
7) shows that the 100-year and the mean-annual floods are highly correlated 
and that the residual error is not appreciably correlated with drainage area. 
Although not immediately evident from the plotting, study of the data shows 
the relation is largely independent of physiography or climate. On the other © 
hand, the plotting in the lower portion of the accompanying chart shows clear- | 
ly that the mean-annual flood is highly correlated with both drainage area and 
physiography (or climate). 

Treating the entire 100 stations as constituting a single population may 
enhance the reliability of the relation between the 100-year and mean-annual 
floods, but the data must be classified by basin, or region, in deriving the re- 
lation between drainage area and mean-annual flood. At least this is the case 
until such time as we are able to explain the basin-by-basin variation in terms 
of measurable factors. Even so, the mean basin-curves in the lower part of 
the chart present an interesting array. Without belaboring the issue, it can 
be said that: 


1. As would be expected, the mean-annual flood increases with increasing 
annual precipitation and runoff (compare Kansas and Ohio River Basins). 


2. Proceeding downstream, the rate of change of mean-annual flood with 
drainage area changes abruptly when a major tributary of different 
climatic and physiographic characteristics enters (see points 27, 42 
and 74, which are for Keokuk, Kansas City and St. Louis, respectively). 
In other words, the relation between drainage area and mean-annual 
flood for mainstem stations may differ appreciably from that for sta- 
tions on adjoining tributaries. 


3. Climate and physiographic characteristics affect the exponent as well © 
as the coefficient in the equation relating area and mean flood. 
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RALPH W. POWELL,! M. ASCE.—Mr. Kuiper is to be thanked for as- 
sembling this large mass of data and giving the profession the benefit of his — 
one hundred frequency diagrams. However, his treatment of the data raises 


1. Visiting Prof., Applied Mechanics, Univ. of Kansas, Lawrence, Kansas. 


DISCUSSION 1616-23 


ASCE 


1000 


S500 


200 


192190 


20 
— 


10 


MEAN ANNUAL FLOOD (X) IN O00 cfs 


2000 


1000 


500 


om Ss 5 IW'OS OOO! Ni(¥) VaNV SOVNIVYG 
ew .O°r fo} fo) [e} 
2 ® s S 8 g ° 3 7 
== T Uke aL Sin eae ee al ieee T T | Lea a iw T 7 
ry 3 


z 
2 
AEA 
BASIN NUMBER 


s 
Yy 


DRAINAGE AREA IN 1000 SQ. MI. 


° © a = 


$19 000'I NI (°°'X) G0014 NVaA-001 


200 
100 
20 


1000 


16 n 


pie | eae 
100 200 


50 


20 


0 5 
MEAN ANNUAL FLOOD(X) IN 1,000cfs 


1616-24 HY 2 April, 1958 


several questions. In the first place, would it not have been better to use the — 
median annual flood, rather than the mean annual flood for the denominator 

of the ratio represented by R? It would seem that if the sample is unbal- 
anced by having too large or too small a number of large floods, it would 
produce more error in the mean than in the median. 

The second question regards using the eight-tenths power of the drainage 
area as denominator of the ratio represented by F. This evidently came 
from Fuller’s formula, but it is doubtful if it applies here. In fact, by 
Fuller’s formula all the R’s would be 1+ 0.8 log 100 = 2.6. The writer 
plotted F against R on semi-logarithmic paper and drew a regression line by 
eye with the equation log F = 2.5 - 0.5 R. There was a great deal of scatter 
and the larger drainage basins fell below the line and the smaller ones above, 
Stations with approximately the same R were then grouped together and for 
each group the mean annual flood plotted against drainage area. It was found 
that the mean annual flood tended to vary as about the two-thirds power of the 
drainage area. Values of K' = = Q/A2/ 3 were then tabulated, where Q = mean 
annual flood in c.f.s. and A = drainage area in sq. mi. This is easily done by 
using tables of two-thirds and three-halves powers from a hydraulics hand- 
book. (Moving the decimal point three or six places in A moves it two or four 
places in A2/3), F* was then plotted against R and the regression line drawn 
by eye has the equation log F' = 3.30 - 0.55 R. The scatter was reduced only 
a little, but there no longer seemed to be any systematic variation with drain- | 
age area. As mentioned above, it is possible that using the median flood 
rather than the mean might reduce the scatter somewhat. However, it is 
quite evident that there are other factors determining R besides the relation 
of the mean or median annual flood to the drainage area. 

As the author points out, it is probable that one of these is the climate as 
represented by the amount of the mean annual rainfall. Perhaps some mea- 
sure of the variability of the rainfall would be a still better criterion. An- 
other factor mentioned by the author is the average permeability of the soil. 
But to try to develop a formula involving these and any other factors would 
require a vast amount of research. It raises the question as to whether it 
might not be just as feasible to try to develop a formula for estimating the 
one per cent flood directly. Attention is called to the fact that Pettis(1) pro- 
posed such a formula many years ago. It would be interesting if some one 
had the time to test this formula against the data the author has assembled. 
More recently Potter(2) proposed a formula for the ten year (ten per cent) 
flood for small drainage areas. Perhaps some day some one can investigate 
whether it would apply to larger drainage areas, or how it would have to be 
modified to do so. 
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THE HYDRAULIC DESIGN OF STILLING BASINS: 
HYDRAULIC JUMPS ON A HORIZONTAL APRON (BASIN I)a 


Discussion by E. A. Elevatorski and A. Rylands Thomas 


E. A. ELEVATORSKI,* J.M. ASCE.—The authors have presented useful 
information on the elements of the hydraulic jump formed in a rectangular 
channel on a horizontal apron. To the designer, the height and length of the 
hydraulic jump are of particular significance in determining the necessary 
dimensions of a stilling basin or channel transition. 

D 
In Figure 5, the authors’ relationship of = versus Fj verifies, without 
1 
doubt, the applicability of the pressure-momentum formula to the hydraulic 
jump. However, the extrapolation of the lower curve seems to be in error, 
since the undular-type jump could not be formed unless the Froude number 
exceeds 2, regardless of the tailwater depth. 

As an alternative to the employment of Figure 5 for finding the downstream 
depth, Dg; a dimensionless relationship between Do, Dj, and the critical 
depth, Dc can be found. This seems to have merit since the Froude number 
need not be computed for each discharge. The following analysis presents 

D D 
the relationship between = and D for a given discharge. 
(B c 

For a rectangular channel, the difference in hydrostatic pressures at sec- 
tions (1) and (2), referring to the beginning and end of the jump, respectively, 
must be equal to the change of momentum per second of the mass flow. 
Equating and transforming, Dj and Dg in terms of q, the discharge per foot of 
channel width becomes 


D, Dy (Dy +D,) = aa (1) 


Since D x = q2 for a rectangular channel, equation (1) becomes 
& 


2 3 
D, Dy (Dy +D,) = 2D, (2) 


To obtain dimensionless ratios, divide each term of equation (2) by Dees 


a. Proc. Paper 1401, October, 1957, by J. N. Bradley and A. J. Peterka. 
* Asst. Prof. of Civ. Eng., Univ. of Arizona, Tucson, Ariz. 
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i 
Therefore D, | Do-) Dy dan, 22 ae (3) 
De De De De : 

nt 

To assist the designer, Figure 1 presents the relationship between D- 

c 

D 


and s for the range normally encountered by the designer. For a given dis- 
c 
charge q(c.f.s/ft. of basin width), the critical depth can be computed from 


equation (4) 


3 
D = ra (4) 


c iy SS 


& 


The longitudinal element of the jump is, without doubt, the most difficult 
element to measure. Certainly this is partially due to the differences in 
opinion of exactly where the terminus of the jump lies. The writer prefers 
the following definition: 


“The origin of the jump is the point where the roller turbulence begins © 
and the water becomes white and foamy, due to the large amount of en- — 
trained air. Downstream, the terminus of the jump is defined as the 
point where no return flow is observable.” 


This seems to be in harmony with the definition of the authors who judged 
the length of the jump as it would best represent the terminus of a conven- 
tional stilling basin. 

Contrary to the conclusions of the authors, the writer is of the opinion that 
the length of the jump is independent of the Froude number. Past studies 

clearly show that the length of the jump, Lj, bears a definite relation to the 
jump height, Dj- To support this, the writer has taken the liberty of using the 
authors’ published data to plot Lj versus Dj = Dg - Dy in Figure 2. Despite 
the fact that the Froude number varied from 2 to 20, a good correlation was 
obtained. From the plot, the correlation can be summarized by equation (5). i 


L, = 6.6 (Dp - Dy) (5) 


Equation (5) compares favorably with the results obtained by other investi- 
gators and listed below: 


Investigator Length of Jump Formula 
Smetana L, —aG (D, - D,) 

wu L; = 6.8 (Dy - Dy) 
Safranez® 10, Sb Hy (D, - D,) 
Page? L, = 5.6 Dy 


The writer wishes to emphasize that equations (3) and (5) are only 
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applicable for hydraulic jumps formed in rectangular channels on horizontal 
L, 


aprons. For trapezoidal basins and channels, the ratio ft is much greater 


J 
than the corresponding ratio for a rectangular channel. Experiments by 


Hsing(1) indicate that a increases as the side slope of the channel decreases. 
j 


It is strongly recommended that model studies be conducted to ascertain the 
elements of the hydraulic jump in a trapezoidal basin before the particular 
design is adopted. 
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A. RYLANDS THOMAS,@ M. ASCE.—The authors deserve great credit for 
their comprehensive study, supported by extensive experimental data, pre- 
sented in the six papers 1401 to 1406. 

The length of the hydraulic jump has been the subject of much discussion 
in the past and, as shown by Fig. 7, considerable differences exist between 
the various series of observations. One reason for this is differences in 
definition of the position of the end of the jump. 

The hydraulic jump is essentially a transition from a high-velocity (super- 
critical) stream to a lesser-velocity (subcritical) stream with a higher sur- 
face level. This involves an upward expansion of flow in which transfer of 
momentum plays a major part, and it is the rate of net upward momentum 
transfer which determines the length of jump. The process of momentum 
transfer is not complete at the end of the surface roller but continues down- 
stream, the rate of upward transfer gradually diminishing as the velocity 
distribution becomes adjusted until the loss of momentum due to bed drag re- 
verses the process into one of net downward momentum transfer. Under such 
conditions of gradual change the difficulty of defining the end of the jump is 
manifest. 

Bakhmeteff and Matzke conceive the end of the jump as the section at 


a. Consultant, London, England. 
b. Bakhmeteff, B. A. and Matzke, A. E., “The Hydraulic Jump in Terms of 
Dynamic Similarity.” Transactions ASCE., Vol. 101, p. 640, 1936. 
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which the surface level reaches its maximum height. This has the drawback 
that it is affected by the surface roughness of the bed and depends on the bed” 
slope—a gently sloping bed might increase the apparent length of jump in- 
definitely. The authors (page 14) have taken the jump as ending at a point ~~ 
downstream either where the high-velocity jet began to leave the floor or on 
the surface immediately downstream of the roller, though which point was 
used in each case is not stated and the former point must have depended on 
the roughness of the bed. [ 

The writer feels that rather than assume a length based on arbitrary con~ 
ditions, on which experimentors are not unanimous, it is preferable to record 
the length of roller, as done by Kindsvater,© a length which is reasonably 
clearly defined, easily observed and a characteristic parameter of the jump. 

The length of jump may be said to depend on the rate of upward momen- 
tum transfer and on the degree or extent of expansion. These may be com- 
plex functions of Froude Number and depth, that is of Fy and Dj and Dg or Dj 
(height of jump = Dg - Dj). Lj is often expressed as a ratio of Dj and, when © 
plotted against Fj, this has the advantage of giving a flatter due in the range! 
of Fy between 2 and 6 than does the ratio Lj/D2 which as may be seen from 
Fig. 7 varies sharply in this range. In fact the Authors’ data when plotted in 
the form Lj/D; : F, indicates a linear relationship over the whole range, 

L; j/Dj reducing from approximately 7.4 with Fj = 2 to 6.0 with F, = 16. 

The Authors’ classification of jumps (Figs. 7 and 9) is of much interest 
and is to some extent supported by experience with canal falls in India. These 
consist of a free overfall followed by a stilling basin to contain the hydraulic 
jump, generally with side expansions. Bank scour downstream has been more 
troublesome when the fall was relatively low (low Fj) than when it was high 
(high Fj). In the former case there is less turbulence to stabilize the flow 
and it is possible that jets might oscillate not only in a vertical plane, as F 
illustrated by the Authors in Fig. 9B, but also in a horizontal plane, resulting 
in the formation of large eddies. Conditions could no doubt be improved by 
adopting a greater width, resulting in reduced discharge concentration and 
higher F1, but width was generally determined by other considerations, for 
example degree of contraction required to give a desired stage : discharge 
relation upstream. 

Insufficient evidence has been produced for the conclusion at the foot of 
page 22 that baffle piers are of little value when Fy lies between 2.5 and 4.5. 
When Fj is low, baffle piers may be less effective in providing back pressure 
but they can be valuable in this range in stabilizing the flow and reducing nes 
length of basin required. 

It would be of value if tests similar to those described in this paper were 
carried out on jumps in channels with side expansions to determine the basic 
relationships in regard to position and length of jump in laterally expanding 
flow. 


c. Kindsvater, C. E., “The Hydraulic Jump in Sloping Channels,” Transactions 
ASCE., Vol. 109. p. 1112, 1944. 
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HYDRAULIC DESIGN OF STILLING BASINS: 
HIGH DAMS, EARTH DAMS AND LARGE CANAL STRUCTURES (BASIN II) 


Discussion by A. Rylands Thomas 


A. RYLANDS THOMAS,» M. ASCE.—The titles of Papers 1402 and 1403 
ind the conclusions therein indicate that stilling basins with chute blocks and 
ond sill (Type II) are recommended for high dams, earth dams and large 
anal structures and shorter basins with chute blocks and end sill but with 
loor blocks in addition (Type III) are recommended for other canal struc- 
ures, Small outlet works and small spillways. In the case of high spillways 
he use of floor blocks might be undesirable because of cavitation but it is not 
slear why in all other cases advantage could not be taken of the benefits of 
loor blocks as evidenced in the Type III design, except in the case of such 
small works that the saving in cost would be immaterial. 

Where cost is the main criterion for the adoption of a Type II design, and 
he choice lies between chute blocks and floor blocks, it may be asked 
vhether a single row of floor blocks might not be more effective than chute 
locks. Comparative effect would be indicated by scour tests in a model with 
nobile bed downstream. 

The observations of jump length are of value but the real criterion of 
ninimum safe length of basin is protection against scour downstream. Ex- 
keriments have shown that depth of scour is reduced as length of basin is in- 
sreased to as much as 12 times the downstream depth. Such length would in 
nost cases be far greater than is required, but the result indicates that 
asin length is a matter of degree of protection required, and this depends to 
;ome extent on the case considered. Model tests with mobile beds give a 
ood indication of scour pattern and comparative depths, but there have been 
ases of deeper scour in the full-scale works than indicated by the models. 

The design shown in Fig. 10 with dentated sill in the middle of the floor 
ppears to be an unnatural compromise between baffle and end sill, not likely 
O be very effective as either. To reduce scour at the end of the floor it is 
nore important to provide an end sill and the Authors are certainly right in 
ecommending that the sill should be at the end (Fig. 14). The function of a 
affle is quite distinct and if a baffle effect is required a Type III design 
hould be used. 

In Fig. 15 the curves are carried to great heights of fall. Where the up- 
tream head at the jump is high a deep and long basin is required. The pos- 
ibility of dissipation of energy in the chute, for example by “splitters,” 


. Proc. Paper 1402, October, 1957, by J. N. Bradley and A. J. Peterka. 


. Consultant, London, England. : 
. Blench, T.: Punjab Irrigation Research Institute, Government of Pakistan. 


Reports for 1946 (p. 52) and 1947 (p. 21). 
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might profitably be explored in all cases when F, exceeds 10 or where it is 
found necessary to excavate the basin appreciably below ground level. 

If Fig. 11 is used for determining the basin floor level, the process is o1 
of trials when Fy, Dy and D2 are not known until the floor level is deter- 
mined. The level may be directly determined by the use of the diagram 
shown in Fig. I of the Writer’s discussion on Paper 1403. 
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HYDRAULIC DESIGN OF STILLING BASINS: SHORT STILLING BASIN FOR 
CANAL STRUCTURES, SMALL OUTLET WORKS, 
AND SMALL SPILLWAYS (BASIN III) 


Discussion by A. Rylands Thomas 


A. RYLANDS THOMAS,» M. ASCE.—The hydraulic design of a stilling 
basin may be undertaken in two stages: the first producing a two-dimensional 
design in which the floor level and length and the arrangement of baffles are 
determined for the worst combination of discharge and tailwater level and 
then tested for other combinations which may occur, and the second expanding 
this design into a three-dimensional basin with wing walls and bank protection 
downstream, taking account of conditions which may produce asymmetry of 
flow. 

The authors’ experimental work and recommendations relate to the first 
stage, in which the designer has to select from a variety of shapes, sizes and 
positions of baffles and end sills. They are to be complimented on their 
methodical attack on the problem, covering a wide range of Froude Number, 
and their results will be of value to designers not only of the works mentioned 
in the title but also of larger river works where surplus energy is to be dis- 
Sipated without danger to the structure. In the following discussion the 
Writer has in view this wider field of application. 

Given a set of conditions the designer’s first step is to determine the basin 
floor level without baffle blocks. Probably in most cases in practice the 
quantities known are the discharge, upstream total energy level, tailwater 
level and total energy level, but not V1, Dj, Dg or F1 which depend on the 
floor level to be determined. If Fig. 11 is to be used a process of trials is 
involved. Curves have, however, been devised for direct determination of 
floor level, for example those of the Central Board of Irrigation and Power, 
India.(1) The Writer uses curves of the type shown in Fig. I; Ey, the differ- 
ence between upstream and downstream total energy levels is known; the 
critical depth Dc is calculated from the discharge per unit width q, and the 
floor level is found by reading E2/D, and deducting E2 from the downstream 
total energy level. In the same diagram curves are provided for E1/D¢e and 
Dj/De, the latter representing height of jump which is of use in estimating the 
length of basin required. D4 may be calculated from Fj and Dc as shown and 
D2 from EQ if required. 

These curves, like those in the Paper, relate to jumps on horizontal floors. 
In many cases a floor with a gentle slope in the direction of flow is advan- 
tageous, in particular where the flow is controlled by gates; should the jump 
begin to move away from the toe of the chute there is more time to take 


a. Proc. Paper 1403, October, 1957, by J. N. Bradley and A. J. Peterka. 
b. Consultant, London, England. 
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remedial action before the jump moves out of the basin and greater protec- 
tion against scour at the lower tailwater levels if the floor is sloping. Ex- 
periments by the Writer(2) indicate that the same curves apply with sufficient 
accuracy to gently sloping floors (the slope tested was 1 in 43) when the floor 
level is taken at its mid-point. 

The provision of baffle blocks or piers in the stilling basin brings three 
benefits: - 


i) Because the blocks create back pressure, the floor may be higher 
with equal margin of safety in stabilizing the jump. 


ii) Because jet expansion is more rapid, the basin may be shorter with- 
out risk of deeper scour downstream. 


iii) Because of the additional turbulence generated in the jump, the down- 
stream flow is more stable, there being less tendency for asymmetri- 
cal flow or large vortices to develop. 


(i) and (ii) bring economies in construction costs. (iii) is of importance in 
the case of canal works in the prevention of bank erosion downstream and in 
the case of large river works in reducing the risk of deep local scour which 
might endanger the structure. 

The value of baffle blocks in enabling a higher floor level to be adopted 
may be seen from Fig. II in which the Authors’ data (Table 4) are shown to- 
gether with those of the Writer.(2) These data are not strictly comparable 
because of differences of relative areas of the blocks and distances from the 
toe of chute. In the case of the Writer’s results the floor had a slope of 1 in 
43 and no chute blocks were provided; the data represent the condition when 
the jump was just held at the toe of the chute. In the Authors’ results the 
floor blocks were nearer the toe of the chute and the tailwater had been 
lowered to the stage when the beginning of the jump had receded from the toe 
and reached the blocks (see Fig. 19); prolongued flow under this condition 
might lead to cavitation damage and dangerous scour downstream. In gener- 
al, however, it can be said that, if the dimensions and location of the blocks 
are reasonably near to those shown in Fig. 17, a reduction of 15 to 20 per 
cent in depth of basin, compared with depth required with no baffles, is 
permissible. 

The Authors would reserve this margin as a factor of safety and design 
the floor level as though there were no blocks. This amount of margin may 
Suit many cases but the Writer prefers to adopt a margin in each case suit- 
able to its needs. For example, where ample safety margin has already been 
allowed in the maximum discharge and downstream retrogression of levels 
advantage may well be taken of the reduced depth of basin required when 
floor blocks are provided. A 10 per cent reduction is represented by the 90 
per cent curve in Fig. II and this may be used for design purposes in appro- 
priate cases. 

The experimental data plotted in Fig. 12 support the conclusion that the 
length of jump is considerably reduced when floor blocks are provided. The 
length of jump is perhaps more directly a function of height of jump, Dj = D2 - 
Dj, than of Dg. A suitable length of basin might be 3 Dj which is easily cal- 
culated from the curve in Fig. I. The minimum length of basin required, 
however, is not necessarily the same as the length of jump as measured by 
the Authors. It is the minimum length with which dangerous scour will not 
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occur downstream under any possible flow conditions. For this aspect of de- 
Sign scour tests are required. 

The Authors’ finding that the most effective floor blocks are those with 
vertical upstream faces confirms the same conclusion reached by the Writer 
as a result of model tests with blocks of various shapes and sizes.(2) Blocks 
of total height equivalent to 2m. of which 1.5m. was vertical face and the top 
0.5m. was sloped back at 33° were found to be no better than rectangular 
blocks 1.5m. high without the sloping top. There are, however, other con- 
siderations which may decide the shape of block in a given case. Where high 
velocities may be encountered precautions must be taken against cavitation 
(Fig. 4 of Berryhill’s paper(3) is surely a warning) and streamlining may be 
necessary. If fish are to be protected from damage the corners of blocks 
should be chamfered or rounded. In such cases the effective area of the 
blocks will be less than the projected area and generally reliance can only 
be placed on the area of vertical face. 

The position of the blocks on the floor is a matter on which guidance is 
needed. The nearer the blocks are to the toe of the chute the higher the 
velocity and hence the greater the back pressure created, but if they are too 
near the jet “hurdles” over the blocks which are then ineffective. This was 
found by Inglis and Joglekar(4) who used a continuous baffle and the same 
effect has been observed by the Writer with rectangular blocks and by Harle- 
man(5) with streamlined blocks. In none of these tests, however, was the 
range of Fj sufficient to establish a general rule for design, such as that of 
the present Authors that the distance from the toe of chute to floor blocks 
should be 0.8 Dg. As dimensions depend on the size of the work it is con- 
venient to express all linear dimensions as ratios of the critical depth De 
which is a suitable linear parameter of the discharge per unit width. The 
optimum distance ratio X/Dc (where X is the distance from toe of chute to 
blocks) is clearly a function of Fj and probably also of h3/Dc. It does not 
depend on Dg though it is possible that if expressed as a ratio of Dg it is con- 
venient for design purposes and sufficiently near the optimum value for 
practical purposes. This, however, cannot be judged from the results pre- 
sented. A curve representing the Authors’ recommendation is shown in 
Fig. Ill together with the ratios pertaining to the basins tested, vide Table 4. 
These are divided into three classes according to the results obtained and it 
will be seen that the best results were given by blocks relatively nearest the 
toe, though no other decided pattern emerges. 

With high velocities, cavitation imposes a limit, so that in a large work 
blocks must be placed at a greater relative distance from the toe than in 
smaller works with equal Fy. Some flexibility in design rules is therefore 
required and it would be of value if the Authors gave the results of tests on 
which they based recommendations, in particular tests made with blocks at 
various distances from the toe of chute, other conditions being the same. 

The optimum height of the floor blocks is related to the depth of jet at 
their position. It therefore depends on their distance from the toe of chute, 
the greater the distance the higher the blocks required for the same effect, 
but increasing the height beyond the optimum is not likely to increase their 
effectiveness and may result in undesirable surface disturbance at lower 
stages. Though the Authors state that the height of the blocks should vary 
with Fy and in Fig. 18 the ratio hg/D is shown to increase with F1, the true 
relative height of blocks in their tests was almost constant, as may be seen 
from Fig. I. The Authors’ rule for block height could be stated in the form 
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hg = 0.5 De. This, however, applies only to the authors’ standard design 
where X = 0.8 Dg. If X is varied hg should be modified, within limits. 

In the absence of data in regard to scour downstream there is insufficient 
evidence for the conclusion that a second row of floor blocks is unnecessary 
and the Writer feels that the second row should not be so lightly written off. 
Besides constituting a reserve in case any blocks in the front row sustain 
damage, the second row may well reduce bed scour and have the effect of 
creating turbulence which helps to stabilize the flow downstream. Where 
streamlined blocks are used to avoid cavitation, a second row may be re- 
quired to provide the desired back pressure. 

The provision of chute blocks is recommended by the Authors but no data 
are produced to indicate their value where floor blocks also are provided. If 
any tests were carried out the results would be of interest. 

The end sill shown in Fig. 17 is an almost essential feature to reduce 
scour at the end of the floor. Deep scour may still occur further downstream, 
where however it is less dangerous. In model experiments by the Writer (re- 
sults unpublished) a continuous sill with two rows of rectangular blocks just 
upstream (these blocks-were of the same height as the sill and were distinct 
from the baffle blocks nearer the chute) were found to be more effective in 
reducing scour than end blocks alone or a Rehbock type sill alone. But the 
continuous sill alone was almost as effective as the sill-cum-blocks in re- 
ducing scour at the end of the floor and, though less effective in respect of 
scour further downstream, is probably sufficient for most cases. 

The height of sill is shown in Fig. 18 as a ratio of Dj. It should more 
properly be considered in relation to Dg (as in Fig. 14) because the height 
required is determined by the flow at the end of the jump irrespective of Fj 
or Dj which represent conditions at the beginning of the jump. In India, 
where the end sill is also known as a deflector,(4) the height is generally 
1/12 to 1/10 of Dg and this has proved satisfactory. The Authors’ sills are 
in some cases within this range and in others somewhat higher. 

As the purpose of a stilling basin is to afford protection against scour, 
the value of the results would be greatly enhanced if supplemented by tests 
with a mobile bed downstream to indicate comparative depths of scour. In 
this connection mention may be made of the research into the design of canal 
falls with stilling basins which has been in progress for several years in 
India. Observations have been made both in models with mobile beds and on 
full-scale structures and accounts of these have from time to time been 
published. (6) 

In this paper a bold attempt has been made to reduce the design of the 
Type II stilling basin to a routine procedure of determining dimensions 
based on two main parameters, Dg and Fj. In the design of complete struc- 
tures it is not always desirable to adhere rigidly to a standard design and 
flexibility is required to meet local conditions. With so many features to be 
determined it is understandable that the detailed results which led to the 
adoption of the recommended dimensions have not been presented, yet with- 
out these it is impossible to see whether in fact the proportions are those 
most suited to requirements or to judge the effect of deviation. Perhaps the — 
Authors will be able to include some of the more important supporting data 
in their reply to the discussion. Even with the data already given, however, 
this Paper will enable an effective initial design to be produced and much 
time to be saved in the experimental work. 
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HYDRAULIC DESIGN OF STILLING BASINS: STILLING BASIN WITH 
SLOPING APRON (BASIN v)@ 


Discussion by A. Rylands Thomas and Earl J. Beck 


A. RYLANDS THOMAS,» M. ASCE.—Sloping aprons offer advantages over 
horizontal aprons in operational control and the numerous experimental re- 
sults presented in this paper are a welcome addition to the rather few series 
of data previously on record. 

The methods of design advocated by the present and previous authors are 
concerned with the determination of tailwater depth at the end of the jump and 
this requires an assumption of jump length, a dimension which is not pre- 
cisely known. A more direct method, not involving any assumption of jump 
length, is illustrated in Fig. A. Where the initial information includes the 
discharge per unit width, q, and the upstream and downstream total energy 
head levels, and therefore Ej, the loss of head in the jump, it is a simple 
matter to calculate Ey,/Dc, where D, is the critical depth, and read E9/Dc 
from the appropriate curve for horizontal or sloping floor in the diagram. 

In this diagram, E2 is the downstream total head above the point on the floor 
vertically below the beginning of the jump, and the required floor level at this 
point is therefore easily determined. 

The curve for the horizontal floor is calculated by the well established 
formulas based on the pressure—momentum theorem. The curves for slop- 
ing floors are drawn by eye from the Authors’ data (omitting those of Flume 
A which had expanding sides) and the Bakhmeteff and Matzke data relating to 
the slopes of 0.046 and 0.070. 

Apart from the experimental scatter of points within each series of data, 
the results of the series are not completely consistent among themselves. 
The points relating to Flume D lie above those of Flume B for the same 
Slope (0.10) and the latter lie rather higher than would be expected from the 
positions of the curves for 0.05 and 0.20. Also the points for the range of 
slopes between 0.16 and 0.263 lie so close together that one curve (0.20) has 
been drawn for all these. This may indicate an optimum slope of about 0.16 
beyond which less advantage is obtained in depth reduction, or it may be a 
result of a particular characteristic of the channel or experimental technique. 
There does seem to be a limit beyond which the flow pattern changes and 
this may occur between the slopes of 0.20 and 0.35, depending on Fj. There 
is considerably more scatter in this range, as evidenced by the data of 
Hickox for a slope of 0.336, possibly indicating instability of flow pattern. 


a. Proc. Paper 1405, October, 1957, by J. N. Bradley and A. J. Peterka. 


b. Consultant, London, England. 
c. Hickox, G. H., Discussion on “The Hydraulic Jump in Sloping Channels, 


by C. E. Kindsvater, Transactions ASCE., Vol. 109, p. 1141, 1944. 
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Whether or not future experimental work will confirm that the curve 
shown for the slope of 0.20 is too conservative and is more applicable to a 
slope of 0.16, it is suggested that in the meantime the results for slopes up 
to 0.20 are sufficiently consistent for the curves to be used for practical de- 
sign. In view, however, of the scatter, a safety margin of 10 per cent of E2 
is desirable in determining the floor level, i.e. the floor should be this 
amount below the level indicated by the diagram. 

Air entrainment can have significant influence on the jump profile and 
tailwater level. It introduces a scale effect which may be of importance 
where the upstream velocity is high.. It would be of interest if the Authors 
would state whether appreciable air entrainment was observed in any of their 
tests, and its effect. 

In discussing the comparative benefits of sloping and horizontal aprons 
(page 28) the Authors suggest that instead of adopting a sloping apron a 
slightly lower horizontal apron be used. This would increase the cost, but 
also the chief advantage of a sloping apron would not be realised, that of a 
visual warning during gate operation. With even a gently sloping apron, when 
it is seen that the jump has begun to recede from the toe of the chute there is 
still time to control the jump before it moves off the apron. With a horizon- 
tal apron this is not so, and damage may occur due to scour and the tail end 
of the apron before the jump can be controlled. 

Case B is more complex than Case D in that an additional variable is in- 
troduced. If points representing the data of Table 9 were plotted in Fig. A 
they would, however, lie between the curve for a horizontal floor and the ap- 
propriate curve for the sloping floor, in positions approximately consistent 
with the relative proportion of jump length on the sloping floor. This appears 
to provide a suitable basis for design, points being interpolated according to 
the ratio of the jump length on the slope to 6E9. 


The authors have not suggested the use of baffle blocks on sloping aprons. 
Such blocks are, however, as useful in this case as on horizontal aprons in 
reducing the depth and length of stilling basin required. Provided the blocks 
are safe against cavitation damage, the most economical design of stilling 
basin in most cases with moderate upstream head is probably one with a 
gentle slope to the floor and provided with floor blocks and end sill. To de- 
termine the length of basin required to afford adequate protection against 
scour, tests are needed in models with mobile beds, and it is hoped that 
future work on these lines will include scour tests with sloping aprons. 


EARL J. BECK, 1! A.M. ASCE.—Application of their experimental data to 
existing apron spillway designs helped the authors substantiate the practical 
value of their recommended design method. 

To extend the range of their comparisons, to spillways of lower head and 
lower Froude numbers, the three existing spillways shown in Fig. 1 were 
re-examined by this same method. These spillways, each of similar sloping 
apron design and founded on sand, were designed with the aid of hydraulic 
model studies.2;3 Although they have low Froude numbers, they have been 


1. Harza Eng. Co., Chicago, IL. 

2. Model Study of the Spillway for the Santee River Dam—US Waterways, Ex-- 
periment Station, Vicksburg, Mississippi, 1944. 

3. Model Tests of Spillway for Petenwell and Castle Rock Dams—Arno T. 
Lenz—Engineering Experiment Station, University of Wisconsin, 1949. 
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operating successfully for many years. The pertinent data for these spill- 
ways is given in Table 1. The results plotted on Fig. 2, further substantiate 
the authors findings. 

The experimental information furnished by the authors would have been 
very helpful in each of these cases for establishing the most advantageous 
apron arrangement and thereby eliminating tests of many less satisfactory 
arrangements. However, model test would still have been desirable if not 
mandatory for furnishing additional information which at best the ex- 
perienced designer can only approximate. For example, model tests of the 
spillways mentioned above were used to determine the affect and critical 
limit of tailwater recession due to downstream degradation, to determine 
length of training walls downstream of apron required to prevent erosion of 
adjacent fill sections, to establish erosion patterns for determining safe 
depth of footings for training walls, to determine maximum uplift pressures 
on apron and walls for varying discharge and tailwater, and to investigate 
special operating conditions. 

These spillways were found to be a very economical design for the sand 
foundations. However, due to the erodability of the foundation material it is 
imperative that such spillways be operated with symmetrical flow. For this 
reason a separate regulating bay was provided in each case to permit small 
releases without operating the entire spillway. This has the added advantage 
of allowing tailwater to rise before opening the main gates. The regulating 
bay was provided with a longer apron and longer training walls. 

Many structures have special problems which can be solved best with the 
aid of a model. In these cases the paper will be a valuable aid to the de- 
signer in establishing preliminary designs for verification by model tests. 

For cases where model tests will not be made this paper will be invaluable 
to the designer. 
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THE HYDRAULIC DESIGN OF STILLING BASINS# 


Discussion by W. H. R. Nimmo and Walter Rand 


W. H.R. NIMMO,! M. ASCE.—This comprehensive group of papers may 
well serve as a manual of design wherever conditions are suitable for one of 
the type of stilling basins proposed by the authors. In the case of two dams 
in Queensland, designed under the direction of the writer, however, such 
conditions did not exist. 

Somerset Dam on the Stanley River is a straight gravity structure for 
water supply and flood mitigation. A central spillway, having a crest 130 ft. 
above the horizontal apron, is divided by piers into eight openings, each 26 ft. 
long controlled by sector gates 23 ft. high. The downstream face of the spill- 
way is joined to the apron by a curve. There are no chute blocks. Overall 
length of spillway including piers, is 253 ft. 4 in. Below each bay of the spill- 
way is a 12 ft. x 8 ft. low level sluice, which has a curved profile and dis- 
charges on to the apron close to the foot of the spillway. At either end of 
the spillway are valve chambers, each containing two 7 ft. 6 in. cone type dis- 
persing valves. The width between the 40 ft. high parallel sidewalls of the 
stilling basin is 340 ft. The sluices and valves are required to rapidly lower 
the reservoir after floods which temporarily submerge a railway bridge. To 
avoid excavation of sound rock—the river bed was originally only 150 ft. 
wide —the sluices were placed below the controlled spillway although when 
construction commenced in 1935 there appeared to be no precedent —at least 
not in this country—for such an arrangement. The maximum discharging ca- 
pacity, with a head of 23 ft. on the spillway, is spillway 100,000; sluices 
50,000; and valves 10,000 second feet respectively. At low flows, tailwater 
level is insufficient to induce a satisfactory jump. At high flows, tailwater 
may exceed 50 ft. and a stilling basin then becomes almost unnecessary. The 
problem was to dissipate the energy of the 70 ft-sec. jets from the 12 ft. x 
8 ft. sluices acting singly, in groups, or all simultaneously. 

Model experiments showed that it was essential to install a continuous 
barrier sufficiently massive to not be displaced by the jets. As constructed, 
the stilling basin has a length of 196 ft. There is a baffle wall 10 ft. high by 
18 ft. wide, 106 ft. from the foot of the spillway slope and a 4 ft. x 4 ft. con- 
tinuous sill at the downstream end. The model experiments showed that under 
the worst conditions the velocity of water leaving the basin will not exceed 14 
ft-sec. which is regarded as satisfactory for the rock in the river bed. 

Assuming that the combined maximum discharge of the spillway and 
sluices, 150,000 sec. ft., is spread uniformly over the 253 ft. 4 in. length of 
the spillway, the discharge per foot of width would be 593 ft. From the charts 


a. Proc. Papers 1401-1406, October, 1957, by J. N. Bradley and A. J. Peterka. 
1. Civil and Hydraulic Engr., Brisbane, N. 8, Queensland, Australia 
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given in the papers, it is estimated that the length of a type I basin would be 
219 ft. If the action of the Somerset Dam basin is less elegant than that of a 
Basin II, it is effective and costs are about the same. 

Tinaroo Falls Dam nearing completion on the Barron River is a gravity 
structure with spillway crest 131 ft. above the level surface of the apron. 
Distance between the side walls of the spillway is 250 ft. With a head of 12 
ft. on the crest, the maximum outflow will be 41,000 sec. ft. Model experi- 
ments showed that at no rate of outflow was tailwater high enought to de- 
velop a jump even with dentated sills. This could have been corrected by 
lowering the floor of the basin but this would involve removing sound rock 
below the elevation of the foundation of the dam itself and was considered to 
be inadvisable. The stilling basin as constructed includes a continuous baf- 
fle wall 8 ft. high by 12 ft. wide, 68 ft. from the intersection of spillway 
slope with the apron and a continuous sill, 38 ft. from the baffle wall. The 
model showed that the basin, though turbulent in action, functioned satisfac- 
torily. It is no longer than a type I basin proportioned by the rules given by 
the authors. 


WALTER RAND, 1 A.M. ASCE.—The effort of Messrs. J. N. Bradley and 
A. J.. Peterka to supplement the 23-year experience of the U. S. Bureau of 
Reclamation by numerous new experiments, and to publish the results with 
the idea to *. .. . emphasize practical design... . so that dependable stilling 
basins may be designed without the need for exceptional judgment or exten- 
sive experience.”, as declared in the foreword of the paper, should be appre- 
ciated by any engineer who deals with the design of stilling basins. These 
papers certainly belong among the best that have been written on this subject. 
However, because of the approach, quoted above, what is offered is not so 
much a generalization of the design of stilling basins, as intended by the 
authors (Paper 1401, page 2), but a standardization of the design. In stress- 
ing structural standards, the hydraulic features of the design have been less 
emphasized, even as much hydraulic data have been given in the form of 
summaries and tables. Such presentation makes it difficult to see why cer- 
tain dimensions of a proposed basin are the most appropriate, and how they 
depend on the flow pattern in the given basin. Some of the proposed basins of 
finished design may be confidently adopted without much judgment or ex- 
perience, and without fully understanding their hydraulic background. If, 
however, due to particular reasons, a somewhat different design is contem- 
plated, it would be difficult to predict the resulting changes in the flow pat- 
tern, and to determine which structural changes may be acceptable. 

In the writer’s opinion, a general design method should consist of a set of 
hydraulic and structural design principles on the basis of which an individual 
design solution can be found for a particular stilling basin. Ordinarily, de- 
sign does not mean that complete structural solutions should be copied. How- 
ever, the hydraulics of stilling basins is not yet so far developed, and the 
typical design solutions, offered by the authors, are of great immediate value. 

As, in any case, it is an advantage to understand better the hydraulic back- 
ground of the proposed stilling basins, it may now be propitious to discuss 
certain aspects of it, One of the basins, for example, is provided with a 
dentated end sill (Basin II). The length of the forced hydraulic jump, pro- 
duced by this sill, is considerably shorter than the natural hydraulic jump. 


1. Asst. Prof. of Civ. Eng., The City College, New York, N. Y. 
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is it exactly of the same length as the basin? Could this forced jump, and 
sonsequently the basin, be even shorter? How long is the basin relative to 
he shortest possible forced hydraulic jump? Why is the chosen basin 

length the best? There is an even shorter basin (Basin III). How are the 
lentated sill of Basin II and the baffle piers of Basin III in their action hy- 
jraulically related? In connection with a general study of stilling basins, the 
writer has investigated the properties of a forced hydraulic jump, that was 
aroduced by a continuous, vertical end sill.“ The results seem to agree 
qualitatively with the effect of a dentated sill or baffle piers, and some of the 
juestions above may be tentatively answered. 

Fig. A summarizes observations on lengths of basins. The ratio of the 
jistance Lg of the sill from the entrance of the stilling basin to the entrance 
Jepth Dj is given as a function of the Froude number F. The shortest pos- 
sible basin length, given by min.Ls/Dj, corresponds to the shortest physical- 
ly possible forced hydraulic jump. The maximum relative distance of the 
sili from the entrance is given by max.Lg/D . In this case, the natural hy- 
iraulic jump of the relative length L/D, occurs in the basin. Between 
min.Lg and max.Lsg any particular sill of a position Lg and height hg should 
ye accompanied by a particular tailwater depth to avoid a sweepout of the 
ump. Any distance L,/D, can be expressed by a coefficient K with a value 
Jetween zero and one. The forced jump that corresponds to K = 0 to K = 0.2 
iS an effective energy dissipator. However, the flow conditions downstream 
Mf the sill would not be acceptable for an unprotected channel bottom. At 
K<-- 0.2 to K = 0.4 an effective forced jump is accompanied by satisfactory 
low conditions downstream of the sill. If the basin length is extended over 
what corresponds to K = 0.5, there would be only little improvement in flow 
sonditions, and the corresponding stilling basin would be of considerable 
ength. At K = 0.8 the basin length is equal to the natural hydraulic jump. 
iowever, as there is an end sill, the hydraulic jump that actually occurs is 
10t a natural one. Only at K = 1 is the basin of sufficient length to produce a 
latural hydraulic jump upstream of the sill. 

Fig. A has been completed by the drawing of lines, representing the dis- 
ance from the entrance of the basin to the end sill or baffle piers, corres- 
Jonding to the discussed stilling basins of type II, I, and IV, and to the 
shortest existing stilling basins of type I, and to the basin length proposed 
yy Forster and Skrinde.3 The plot impresses one with its diversity. All 
hese are acceptable design solutions for basin length, provided the height of 
he tailwater and of the sill are adjusted to the given length. These solutions 
ire not contradictory as an analysis of flow conditions will show. Basin II, 
hat corresponds to K = 0.4 to K = 0.5 makes use of an effective forced hy- 
lraulic jump accompanied by a relatively tranquil flow pattern downstream 
f the sill. This is an excellent design solution. The shortest existing type I 
tilling basins correspond approximately to K = 0.1 to K = 0.2. Acceptable 
lesign solutions are possible also in this range. The corresponding forced 
ydraulic jump is an effective energy dissipator, and if the height of the sill, 
ind the given tailwater depth, produce acceptable flow conditions downstream 
f the sill, depending on the erodability of the channel, the basin could still be 


satisfactory. 


. “An Approach to a Generalized Design of Stilling Basins” by Walter Rand, 
The N. Y. Academy of Sciences Transact. Vol. 20, No. 2, 1957. 

. “Control of the Hydraulic Jump by Sills” by J. W. Forster, and R. A. 
Skrinde, ASCE. Transactions, Vol. 115, p. 973, 1950. 
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The baffle piers of Basin II, corresponding to K = 0.0 to K = 0.05, pro- 
duce a jump that is practically the shortest physically possible. These piers” 
are close to a critical position (min. Lg), and also their height is close toa ~ 
critical value as shown later. The notion of the authors that there is no 
Sat walere intention to give the impression that the position or height of 
thé baffle blocks are critical” (Paper 1403, p. 16) results from their different 
experimental approach. Of course, the basin has to be extended beyond the 
baffle blocks to protect the channel against harmful erosion. Because the 
forced jump is produced by the baffles, the end sill serves only to improve 
flow conditions downstream of the piers, and it can be smaller. 

Fig. A also shows that Basins II and II become relatively short (smaller 
K values) at smaller Froude numbers, and no solution is therefore possible 
for F< 4. This leads to Basin IV. Would it not be possible to have a con- 
tinuous variation from the length of Basin II to the length of Basin IV? 
Forster and Skrinde apply such a trend. However, their basin is relatively 
long at larger Froude numbers. ; 

This analysis shows that the range and properties of the forced hydraulic 
jump can be used to compare the relative length of different stilling basins. 

It also shows that there is no basin length that is the only correct one. Many 
choices are possible in a certain acceptable range, provided the height of the 
sill and tailwater, and the erodability of the channel are correctly related. 

Some observations will be offered about these relationships. ' 

The concerted action of the sill and tailwater keep the jump from being 
swept out from the stilling basin. If the problem of erosion downstream of 
the sill is provisorily left out of consideration, then a higher sill could be 
used with a lower tailwater to stabilize the jump in the basin. It is possible 
to define the height of the sill as a function of tailwater, and to determine the 
smallest sill that is sufficient to stabilize the jump without any collaboration 
from the tailwater (in the case of very low tailwater). With this critical sill, 
the required depth of the tailwater depends only on the degree of erodability 
of the channel. 

The authors have another approach. They recommend a tailwater of con- 
jugate depth Dg, that is sufficient to stabilize the jump without any action of 
the sill. This recommendation is partially based on a statistical average of 
the tailwater depth at existing structures, equal to 0.99 Dg. As sometimes 
an unnecessarily high tailwater is naturally present, this average may not 
indicate what is actually needed. 

The height of the sill of Basin II is equal to 0.2 Dg (Paper 1403, p. 7). It 
is relatively low, especially at larger Froude numbers, if compared to the 
critical height of the sill, that would stabilize the jump without any action of - 
the tailwater. Therefore the lowest tailwater depth that prevents the sweep- 
out is close to D2, and relatively closer at larger F values (Paper 1402, 

Fig. 11). The effect of a high sill is well illustrated by the baffle piers of 
Basin II]. Their height is close to the critical, and a lower tailwater stabi- 
lizes the jump. | 

By the limitation of the height of sill to 0.2 Dg, which is an average value 
based on the existing stilling basins, and in no shown relationship to the basin 
length or tailwater, a rather rigid standard has been established. Could not . 
more flexibility be provided? | 

The writer hopes that these observations further illustrate the hydraulic — 
aspects on which the design of stilling basins depends. 
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FRED W. BLAISDELL,* M. ASCE.—Messrs. Bradley and Peterka have 
performed a most useful service to the designer in presenting a number of 
generalized stilling basin designs in this series of papers. For the smaller 
structure, where a special model study is usually considered too costly, these 
papers make available a number of carefully tested designs. For the larger 
structure, where an individual model may be justified, these designs will pro- 
vide valuable basic information if they do not actually eliminate the need for 
the special study. 

The generalization of the results obtained from model studies on specific 
structures, supplemented by additional laboratory research to fill in the gaps, 
is something that has been long talked about but on which little has been ac- 
complished. From these examples, it is hoped that others will see the result- 
ing significant benefits and will make similar analyses of additional struc- 
tures. 

The authors show an equation for the curve representing E,/D1 in Fig. 8. 
They do not mention that this is a theoretical equation derived through the use 
of the energy equations 


v, 
E, = 2g ots D, 
and 
5 
E, = 2” + Do ? 


the equation of continuity, and Eq. 3. The theoretical equation representing 
the curve of Ey,/E1, derived from the same basic equations, can take the form 


Sag NS EE Vane caterer a 2(-la \SEoh 1)|2 
1 F +2 


The curve of Ey,/E, in Fig. 8 does not represent Eq. A precisely, so the 
writer has shown it in Fig. 8A. The very good agreement between the experi- 
mental data and the theoretical curves of Fig. 8a is not unexpected in view of 
the good agreement between the theoretical and experimental conjugate depths 
shown in Fig. 5. 

The proportion of the energy entering the hydraulic jump which is lost in 
the jump is seen in Fig. 8 to increase with the Froude number. The advantage 
of high Froude numbers is apparent if a high degree of energy dissipation is 
desired. This is true not only for the hydraulic jump stilling basin but also 
for those basins where the “jump” is assisted by means of blocks, piers or 
sills. 

It is sometimes possible to exercise a degree of control on the Froude 
number. A transition between a chute or a culvert and a stilling basin may be 


a. Project Supervisor, U. S. Department of Agriculture, Agricultural Re- 
search Service, Soil and Water Conservation Research Division, St. 
Anthony Falls Hydraulic Laboratory, Minneapolis, Minnesota. 
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ised to spread the flow into a thinner sheet. This may cause some change in 
V1 but the major change will be a decrease in Dy and a resulting increase in 
Fy (see Eq. 1). Based on certain assumptions, the writer has shown(!) that 
reductions in the concrete and excavation volumes for chute or culvert outlets 
us high as 62 percent and 41 percent respectively can be achieved. Cost re- 
juctions might be even greater since the net reduction in concrete volume is 
obtained by a large reduction in the more expensive wall concrete and an in- 
srease in the less expensive floor concrete while the net reduction in excava- 
‘ion is achieved by widening the excavation and decreasing its depth, the 
Jeeper excavation being more expensive in the usual case. The use of a 
Yransition to increase the Froude number will frequently be economical for 
sulverts, outlet works and chutes, but will probably not be practical for rela- 
ively wide and short structures such as dams. Further advantages resulting 
rom an increase in the Froude number are that a larger percentage of the 
energy is dissipated, as shown in Fig. 8, and it may be possible to use the 
shorter Basin I instead of Basin IV or Basin VI. 

When discussing Fig. 9 the authors mention the pulsating action observed 
‘or Froude numbers between 2.5 and 4.5. The writer made a similar observa- 
ion during the studies leading to the development of the SAF stilling basin(2) 
‘or a range of the Froude number from 1.73 to 17.3. During these studies it 
was noted that the best performance of the stilling basin was achieved at 
Froude numbers greater than about 4.5. Under the heading, “Practical Con- 
siderations” in Paper 1401, the authors suggest that the 2.5 to 4.5 Froude 
lumber range sometimes can be avoided “by altering the dimensions of the 
structure.” The writer agrees and suggests a transition as one means of 
naking this alteration. 

In Paper 1402 the authors state that wingwalls are not amenable to generali- 
zation. The writer is not in agreement with this broad statement. Three 
feneralized stilling basins have been studied under the direction of the 
yriter(2,3,4) and in each basin wingwalls having a sloping top and flaring in 
lan have been found to be superior to other types of wingwall in protecting the 
lam fill from scour. They, in combination with a good stilling basin, have 
een found so effective that the writer ordinarily suggests that no riprap be 
ised in the downstream channel to prevent bed and bank scour, except possibly 
or that bank scour and beaching caused by surface waves. Experiments(5) 
onducted by T. M. Brett for the Soil Conservation Authority of Victoria, 
Australia verify the effectiveness of the sloping top, flaring wingwall. 

Basin III described in Paper 1403 and the SAF stilling basin are similar in 
ppearance in that chute blocks, baffle piers and an end sill are used. The 
asin dimensions, however, are considerably different in some particulars. 
fhe similarity may be due in part to the fact that the principal features of the 
IAF stilling basin were obtained from a rectangular stilling basin design re- 
jorted by Jacob E. Warnock, M. ASCE,(6) a former associate of the authors 
t the Bureau of Reclamation Hydraulic Laboratory. The differences perhaps 
an be ascribed to differences in philosophy of those who developed the pro- 
jortioning rules. The authors apparently proportioned Basin III to insure that 
he energy dissipation took place entirely within the stilling basin. From the 
hotographs of Flumes B, C, D and E in Figs. 1, 2 and 3, it appears that the 
uthors’ tests were performed in fixed-bed channels. The SAF stilling basin 
ests were performed in moveable bed channels where the bed scour was used 
S a measure of the stilling basin performance. No attempt was made to ac- 
omplish the entire energy dissipation within the stilling basin if the scour or, 
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more properly, the lack of scour was such as to not endanger the structure. z 
That some energy dissipation takes place outside the SAF stilling basin is 
shown in Fig. A where it can be observed that a standing wave may occur 
downstream of the basin exit. Since this standing wave is on the surface, the 
end sill protects the stream bed, and the basin performance is satisfactory, 
it is not felt necessary to lengthen the basin and thereby increase its cost. 

This difference in philosophy and experimental methods has resulted in 
the SAF stilling basin being much shorter than Basin II]. The comparison is 
made in Fig. 12A. Relatively longer SAF stilling basins are required at low — 
Froude numbers. This is in contrast to the relatively shorter lengths re- 
quired at low Froude numbers for Basins I, I and IL The writer feels that 
the relatively shorter basin lengths required for the SAF basin as the Froude 
number increases reflects the greater percentage of energy dissipation that 
occurs naturally in the hydraulic jump with increasing Froude numbers, as is 
shown in Fig. 8A. 

The authors locate the Basin III baffle piers 0.8D2 downstream from the 
chute blocks. The basin length is given in Fig. 12 as from 2.2D9 to 2.75D9. 
The baffle piers are thus located from 1/2.75 to 1/3.34 times the basin length 
downstream from the chute blocks. The corresponding distance recommended 
for the SAF stilling basin is 1/3. The agreement is close. i 

A difference is that the baffle piers in Basin III do not have to be staggered 
with the chute blocks. Staggering is recommended for best performance of 
the shorter SAF basin. With this one exception, the writer emphatically en- 
dorses the authors comments concerning baffle piers. 

The method used to determine the elevation of the stilling basin floor illus- 
trated in the “Application of Results (Example 3)” deserves emphasis. It has 
been the writer’s experience that many designers consider that the basin 
floor or end sill elevation should be at stream bed level. This is not neces- 
sarily the correct elevation. The correct elevation is determined by sub- 
tracting the required basin tailwater depth from the elevation of the natural 
tailwater downstream from the stilling basin. The authors have correctly 
shown that the floor elevation should be checked at maximum and lesser dis- 
charges. 

The recommended width and spacing of the baffle piers is shown in Exam- 
ple 3 to be 25.5 in. The basin width is computed as 50 ft. from the total dis- 
charge and the discharge per unit width. This results in about 23-9/17 
spaces. An even number of spaces is required if a half space is to be located 
next to each sidewall. There will be 24 spaces if the baffle pier width and 
spacing is 25 in. A baffle pier 25 in. wide can then be centered on each odd 
numbered 25-in. distance from the basin wall. 

In the case of the SAF stilling basin, the baffle piers are specified to occu- 
py from 40 percent to 55 percent of the basin width. One reason for giving 
this latitude was to permit the designer to select sizes that would eliminate 
fractions of inches. Presumably some such latitude is possible in Basin Il, — 
but the writer agrees that, for best performance, the piers should occupy 
about 50 percent of the basin width. é 

In the section headed “Effect of Slope of Chute” in Paper 1405, mention is — 
made that “The baffle piers on the floor [ of the Type II basin] tend to alter 
the asymmetrical jet, resulting in an overall improvement in operation.” The 
writer observed similar results for the SAF stilling basin where the depths at 
the basin entrance were non-uniform. In one case, the circular pipe outlet 
discharged onto a diverging chute. At the basin entrance the jet was high in © 
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ie center and there were standing waves at each sidewall. In spite of this 
sry poor depth distribution, the scour pattern at the basin exit was complete- 
‘normal. This is ascribed to the action of the baffle piers in correcting for 
ie non-uniform depths at the basin entrance. Attempts should be made to 
2cure uniform depths and velocities in the approach to the stilling basin, but, 
they are not uniform, the baffle piers of Basin III or other similar basins 
ill help correct the adverse condition. 

No single stilling basin can be expected to best fit all the many combina- 
ons of discharges, velocities, tailwater conditions, bed materials, sites, 
isks to life and property, and other tangible and intangible conditions that 
just be considered by the designer. The selection of the most suitable 
‘illing basin for detailed investigation from those types presented by the 
uthors is facilitated by the comparison of the important features of each 
asin summarized in Fig. 47. There are, of course, other types of stilling 
asins that were not studied by the authors, but this in no way detracts from 
1e value of their contribution. 
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HYDRAULIC DESIGN OF STILLING BASINS: HYDRAULIC JUMP ON A 
HORIZONTAL APRON AND HYDRAULIC JUMP ON A SLOPING APRON2 


Discussion by Solano Vega-Vischi 


SOLANO VEGA-VISCHI,! A.M. ASCE.—In analyzing the data obtained by 
the U.S.B.R., instead of using the Froude Number as the parameter for relat- 
ing the data, a different approach was attempted, making use of the critical 
depth as a factor in defining the dimensionless relationship of the different 


features of the phenomenon. This relationship is expressed as the ratio a L 
C 
where do is the depth of the flow entering the hydraulic jump and d, is the 

critical depth for that discharge. 


A correlation between the Froude Number and the parameter a can be 
C 
established. In the following equations this parameter will be noted as Xp. 
The Froude Number is defined as: 


where F is a dimensionless parameter, V and d are velocity and depth of 
flow respectively entering the jump, and g is the acceleration of gravity. 
The correlation of both parameters can further be demonstrated as 


follows: V d V 
ee Nee 
Vgdg’ c 


fo) 


1. Proc. Papers 1401 and 1405, October, 1957, by J. N. Bradley and A. J. 
Peterka. 
|. Cons. Eng., Santiago, Chile. 
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In this form it is evident that the parameter X, still retains the same ad- 
vantages as the Froude Number with respect the dynamic similarity of the 
models, and further allows the use of simple analitical expressions in the 
foundamental equations. 

The theoretical aspects for the derivation of the theory using the critical 
depth criteria has been taken from the treatment given by F. J. Dominguez. 2 
This derivation is explained in the following paragraphs. 

Referring to figure 1 which represents a rectangular channel with a slope 
o@, the theorem governing the variation of the linear momentum is applied 
over a projected length above the horizontal axis Y-Y, over which the rate of 
change of momentum and the forces producing that change are equated. 


Notation of symbols used in the following equations (reference to figure 1) 


Vo= velocity in the initial section of the hydraulic jump 
V,= velocity in the final section of the hydraulic jump 
d = depth of the supercritical flow at The initial section 


of the jump 
d, = depth of the subcritical flow at the final section of 
the jump. 
a = angle of the bottom grade with a horizontal plane 
L = length of the jump measured on the bottom grade 
L,= length of the jump measured ona horizontal plane 
g = acceleration of gravity 
q = flow of water per unit width of channel, in the same 


units of measurement used ing 
w = unit weight of water. 


The increment of linear momentum projected on Y-Y is: 


g4(i-Vo) cosa 


The difference of hydrostatic pressure in the initial and final sections 
over Y-Y is: 


W cos?a (d5- d?) 


In addition to the external horizontal hydrostatic pressure forces acting . 
on the water enclosed by the initial and finalsections, the horizontal q 


2. F. J. Dominguez.—Curso de Hidraulica.—Santiago de Chile, 1945. 
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Omponent of the reaction caused by the water’s weight on the bottom of the 
hannel, should also be included. This component force, projected onto Y-Y 


iE 
4 wf hdL sina 
0 


This integral can be simplified by substituting an average value h, for the 


nean integrated values of h. Substituting Li for L, the value of the 


$Q 


ntegral becomes: 


wheL, sina 


en wheL, tana 


The equation for the hydraulic jump can then be written: 


1 7 4Vi- Vo) cosa = 5 wcos-a(do-d;) + wheL, tana 
2 2 
2) (g,~ a) 0080 ~2S 2 (eE-dF) + ht, tana 


2 
but + = 2 in which dg is the critical depth 


2 
Dividing equation (2) by d C 


5 (- )cosc ie cg oe + Tana ; 7 


eed. d- de 
Noting: 
di Ne. 
(i Reeih dk H 
do | Le 
ioe ae pee 


squation (3) can be written: 


ates - + costa (X2- X) + HL, tana 
4) (x. x,) cosa > COS*A(X- Aj c 
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If the average value for h, is represented by the arithmetic mean 


5 (dot d,) , equation (4) becomes: : 


a Calan ae 2 2 mye a8 
(x x )cosa 2 COSA (Xo- Xi) + (Xo of XC ran 


or ; 


cos Q 
(5) x, 


ae +(X,cosa) + >Xol-tane 


+ £(X cosa’ - > X,Letana = 


Since the influence of cos for slopes up to 0.200 is small, the use of 
equation (5) can be restricted to slopes above 0.200. 

The following simpler equation can be appiied for slopes between 0.000 
and 0.200: 


| | v2 | | ly2 / 
—— —_— — — °° _ — + — 

In applying equation (5) to slopes above 0.200, and considering that the 
tailwater level was maintained horizontal for the test conditions, (reference 
to figure 2), equation (5) becomes: 


(7) x t2X-% L.tana = Yaa £(X,cosa)” he > Xol¢ tana 


Length of jump. 


It is stated in Proc. Paper 1401, “the length of jump was the most difficult 
measurement to determine. Where chutes or overfalls were used the front of 
the jump was held at the intersection of the chute and the horizontal floor. 
The length of jump was measured from this point to a point downstream 
where either the high velocity jet began to leave the floor or to a point on the 
surface immediately downstream from the roller, whichever was the longer. 
In the case where the flow discharged from a gate onto a horizontal floor, 
the front of the jump was maintained just downstream from the completed 
contraction of the entering jet. In both cases the point at which the high , 
velocity jet begins to rise from the floor is not fixed, but tends to shift up- 
stream and downstream. This is also true of the roller on the surface. It 
was at first difficult to repeat length observations within 5 percent by either 
criterion, but with practice satisfactory measurements became possible. It 
was the intention to judge the length of the jump from a practical standpoint; 
in other words, the end of the jump, as chosen, would represent the end of the 
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concrete floor and side walls of a conventional stilling basin.” 

Table I includes 94 U.S.B.R. test results for slopes between 0.050 and 
0.280. (Table 8 Proc. Paper 1405). The experimental values for do, dy, Ly, 
and q were introduced into each side of equations (6) and (7) to check the 
correlation with the equations, and to determine the acceptance of the hypo- 


thesis used for solving W hd L . The computed difference between the two 


sides of the equation for any given set of experimental values was computed 
as a percentage referred to the average value of the two sides of the equa- 
tion (table I col. 14). The average variation for all the tests was 4.45%, in- 
dicating it to be within the same range of variance usually expected in the 
experimental observations. It appears, then, that the correlations are suf- 
ficiently close to substantiate the mathematical analogy. 

In an attempt to correlate the length of the jump with the degree of slope, 
plottings were made using various parameters and different coordinates sys- 
tems. The plottings where the parameter Xp was used in relation to Lc and 
the slope tengent showed some apparent consistency for values of Xp be- 
tween 0.200 and 0.500, but became erratic for values of Xp outside that range. 
Plottings on natural and semi-log scales and the use of others parameters 
still failed to show a consistent relationship. 

When the test values were plotted on log-log paper, using the parameter 
Lc as ordinates and X, as abscissas, the correlation proved to keep some 
jefinite relationship (figure 4). The test slopes varied from 0.05 to 0.280; 
and the plotting appeared to follow a definite pattern, showing a consistent 
relationship between Lc and Xo, irrespective of the slope. The arc of an 
sllipse very suitable approximated the locus of the plotted points and offered 
a Simple analytical medium of definition in the form of a logarithmic ex- 
pression. 

The percentage of error of the curve in relation to the actual values mea- 
sured did not exceed about 5.07%, still within the range of variance expected 
in experimental observations. The expression for the assumed curve is: 


Log Le = 1.36VI-(Log 10X,)° (8) 


(Log Le and Log 10Xo are common logarithms) 

This relationship is to be considered valid between Xp = 0.100 and 0.700, 
which ordinarily includes the practical accepted range for definable hydraulic 
umps. 

From an academic point of view the function linking Le and Xo might be 
‘riticised, since it does not fulfill the condition of having a value of zero for 
he extreme case Xo = 1 (instead of that, Log Le = 0). However, considering 
hat this relationship has been assumed to be acceptable only in the range of 
he experiments, namely, for values of Xp between 0.100 and 0.700, and that 
vhen Xp values approach to 0.700 (F = 1.7) the jump will lose its stability and 
ts length will become indefinable, the above mentioned criticism can be over- 
ooked, and the function accepted within the aforementioned range. 

Taking the data from the model tests and substituting the length of jump L' 
omputed from the curve whose expression is equation (8) (figure 6), new 
ralues for each side of equation (5) and (6) were computed. In determining 
he value for Xj, the tail water depth dj also was modified to reflect the new 


1616-64 HY 2 April, 1958 


length, this new value becoming: : 
d}=d,+(L-L) tana (igure 2) 


Similarly to the comparison of the original values, the computed differ - 
ences between the two sides of the equation were compared to the average 
values of the two sides of the equation to arrive at a percentage variation. 
This percentage variation was slightly less than the original values, showing” 
an average error of 4.29% (Table I, col. 16). 

As can be noted in the expression for Le given in equation (8), the length © 
of the jump is independent of the slope, Le being a function only of the depth © 
of flow at the beginning of the jump related to the critical depth of flow. 

In order to substantiate the hypothesis that L, and Xo are interrelated in- 
dependently of the slope, the experiments performed on horizontal grades 
were plotted in the same manner as those for the sloping channels. (figure 
5). Superimposed on these plottings is the curve defined by equation (8). 

The consistency of the curve and plotted points are noted. The test data for — 
horizontal grades are tabulated in Table II (from Table I Proc. Paper 1401). 
The difference between the recorded and computed lengths by equation (8) for 
117 tests shows an average difference of 5.86% (Table II, col. 15). 

On the basis of using the value of L, in equation (8), equations (6) and (1) 
were solved for values of dp and dj and plotted for various values of the slope 
(tana). The graph is given in figure 7. Having given the depth of either the 
supercritical or subcritical flow, the conjugate depth of the hydraulic jump 
can be determined. Then from figure 6 the value of Lj can be ascertained. 
By trial and error, the position of the jump can be Selected which will satisfy 
the initial condition of momentum at the start of the jump, and the tailwater 
condition below the jump. 


Case of jump formed both on the slope and the horizontal portion of the apron. 
(Case B, Proc. Paper 1405) 


As to the case in which the jump starts on the sloped part of the basin 
(figure 3), the treatment has been for the intermediate case between hori- 
zontal and sloping apron. 

From Table 9 in Proc. Paper 1405, pees have been taken to form Table 


Ill of this discussion. In this table the parameter Xo= a has been com- 
puted from the experimental data and appears in column & The partial engl 
of the jump formed on the slope is noted as 1 as in the original paper. The 
length L of the jump have been obtained from Figure 6. The ratio 1/L is 
shown for each case in column 9, Table II. ‘ 
The conjugate depth for do, namely dj, is supposed to consist of two : 
terms: one dj, as if the whole jump were to be formed on a sloping apron, 
and a second term dj as if the whole jump were to form itself on the hori- 
zontal apron. - 


d=d)x [- xcosta+d} (1 ) @) 
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r divided by d, 


OE eit Tis +: x costa+X\ (I- + ) (10) 


30, has been introduced for sloping aprons equal or above 


The factor cos 
1.200. : 

In columns 10 and 11 are the values of X, and X] taken from figure 7 for 
he different slopes, and in column 15 the resulting conjugate depths. In 
-olumn 16 the actual conjugate depths (tailwater) for each experience. The 
lifference in % is the last column of the table. It is shown that the average 
lifference in the 103 experiments, for the computed tailwater according with 


his theory and the actual values measured in the experiments is only 4.1%. 
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FIGURE 5 
TESTS WITH HORIZONTAL 
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THE HYDRAULIC DESIGN OF STILLING BASINS: SMALL BASINS FOR 
PIPE OR OPEN CHANNEL OUTLETS—NO TAIL WATER REQUIRED 
(BASIN VI)# 


Discussion by John R. Argue 


JOHN R. ARGUE.!—The Bradley-Peterka impact-type stilling basin was 
designed for use in Irrigation works where flow from the structure is car- 
ried away by a trapezoidal channel with invert at sill level and sides sloping 
at 1:1. Where the inflow into the structure is carried in a pipe, head in ex- 
cess of 2 diameters (above invert) would be extremely rare. The pipe would 
usually be running partly full (i.e. open channel flow in the pipe). The 
Bradley-Peterka design shows a rigid relationship between maximum dis- 
charge and basin dimensions and no variation in these dimensions in given 
with different erodible materials downstream from the structure. 

A similar problem to that considered by Bradley and Peterka, has con- 
fronted Local Government engineers in New South Wales, Australia. It con- 
cerns outlet structures for road pipe culverts. 

In many parts of the State, arid lands are occasionally subject to heavy 
rainfall, which causes gullying and severe erosion. It is common for a pipe 
culvert during flood conditions, to be subject to a head of 3 diameters above 
the upstream invert. Where discharges of this magnitude are passed onto 
erodible soils, severe damage is caused to both the roadway and the terrain 
downstream from the structure. 

The existing standard structure consists of a headwall and two wing walls 
set at 45° to the centre line. When this structure operates at high discharge 
under open channel flow conditions in the pipe, a fast moving jet passes di- 
rectly from the pipe onto the erodible bed. The structure itself provides no 
energy dissipation—the full force of the jet is taken by the moveable material 
downstream from the structure. 

The severe scour which occurs around this structure is entirely due to the 
strong eddies which form on either side of the fast moving jet. These two 
eddies rotate in opposite directions—the left anticlockwise, the right clock- 
wise —both being “powered” by the centre line jet passing from the pipe. 

A successful design must, therefore, attack this problem at the root cause 
xy ensuring that flow onto the moveable bed is not concentrated into a fast 
moving jet, but evenly distributed at a tolerable velocity across as wide a 
ront as possible. This is undoubtedly the criterion for the majority of hy- 
jraulic structures which discharge onto moveable beds. 

In pipe culvert outlet structures it is most important that the extent of 
3cour downstream from the structure be reduced as much as possible. It is 


1. Proc. Paper 1406, October, 1957, by J. N. Bradley and A. J. Peterka. 
|. Engr., Harbors and Rivers Branch, N. S. W. Dept. of Public Works, 
Sydney, Australia. 
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desirable to have the scour hole confined within the boundaries of the road- 
way reserve. Providing that the scour hole does not encroach upon the toe 
of the roadway embankment, the width of scour is immaterial. A squarish 
hole rather than an elongated hole would, therefore, best suit the condition 
required. 

For pipe culvert structures, it is preferable to have the pipe running full 
(i.e. pipe flow). The advantages which pipe flow holds over open channel 
flow in the pipe are threefold: - 


1. Greater discharge for a given head; 
2. Smaller average velocity; 
3. More even velocity distribution across the pipe. 


To investigate the problem, a research programme using a 4" pipe culvert 
model with downstream pipe invert at natural surface level is being carried 
out at the Hydraulic Laboratory of the Public Works Department, N. S. W. 
The object of this investigation is to determine a satisfactory outlet struc- 

‘ture for all pipe culverts with invert at natural surface level, in the proto- 
type range 18"' to 72". 


#» 


When details of the Bradley-Peterka impact-type stilling basin were pub- — 


lished an interesting comparison was made between this American structure © 


and that which has evolved, to date, at the P. W. D. Laboratory. The com- 
parison was made on the basis of a 54'' diameter pipe, 28 feet long, sloping 
at 1% grade. The Bradley-Peterka structure model corresponding to this 
size was built in accordance with the details given in A.S.C.E. Paper 1406. 
The alternative end sill design was used in order to reduce scour toa 
minimum. 

Comparative scour tests were carried out in a moveable bed of Sydney 
sand—a fine grained, yellow, non-cohesive beach sand. 

The following discussion gives details of the comparison:- 


Testing Procedure: 


To assess the effect of tailwater on the functioning of the structures, two 
sets of tests were carried out as follows:- 


a) Tailwater at 0.5 feet (prototype) above sill level; 
b) Tailwater at 1.0 feet (prototype) above sill level. 


This level was maintained constant, during testing, by a tailgate and a 
point gauge located in the still water clear of the structure and near the toe 
of the embankment. 


Each test was run for a model time duration of 30 minutes, the build-up of 


discharge taking the first five minutes of this period. The 30 minute test 
duration represents approximately 85% stability of scour. 
Both structures were tested with discharges ranging from 115 C.F.S. to 


¢ 


f. 


4 


7 


330 C.F.S. The lower discharge corresponds to the head at which pipe flow 


would commence, i.e. 115 C.F.S. is the maximum open channel flow discharge. 


The latter discharge (330 C.F.S.) corresponds, approximately, to a head of 3 
pipe diameters. 

The maximum design discharge nominated for the Bradley-Peterka 54" 
diameter outlet, is 190 C.F.S. Thus, it was necessary for the test to cover a 
range which is nearly 75% greater than the design maximum for the Ameri- 
can structure. 


——————L ee 
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tESULTS 


COMPARISON OF STRUCTURE FUNCTION 


tating. The most obvious difference between the two outlet structures was 
he variation in discharge. For a given headwater level, the Bradley- 
eterka structure given approximately 10% less discharge —at high flows— 
han does the P.W.D. design (see fig. 1) this is probably due to the back 
ressure from the water throttled between the hanging baffle and the head- 
vall in the Bradley-Peterka outlet, reducing the head difference through the 
jipe. 


Jpen Channel flow (in pipe) 


3radley-Peterka structure: 0< Q < 240 CFS. 


The hanging baffle had no effect until the flow reached 220 CFS. Up to 
his discharge the flow became progressively more concentrated towards the 
entre of the end sill. At this discharge (220 CFS), the hanging baffle began 
o take effect. At 240 CFS, the baffle caused self-priming—it being impos- 
sible to run the model at open channel flow for higher discharges. 


>.W.D. Structure: 0< Q< 330 CFS. 


The P.W.D. structure will run at open channel flow over the whole range 
f heads up to 3 diameters (270 CFS). Up to 120 CFS the flow across the end 
jill is evenly distributed. Beyond this discharge, the flow, at high velocities, 
comes progressively more aerated at the centre. 


ipe flow 
radley-Peterka structure: 115< Q < 220 CFS. 


The hanging baffle again appeared to have very little function within this 
‘ange. Flow across the end sill was fairly evenly distributed. 


220< Q< 330 CFS. 


At 220 CFS the hanging baffle began to take effect and separation at the 45° 
nd walls commenced. Owing to this separation, high flows across the end 
ill tended to become concentrated into a band equal in width to the parallel 
yalled section of the outlet structure (see fig. 7). 


.W.D. Structure: 115 < Q < 330 CFS. 


This structure diverges the flow and a steady hydraulic jump and “surface 
Oller” occur at all discharges. Flow across the end sill is fairly evenly 
istributed throughout the full range of discharges. Aeration occurs for all 
lows above 270 CFS. 


1616-80 HY 2 April, 1958 
COMPARISON OF SCOUR PATTERNS : 


1. Maximum Scour Depth (see figs. 2A, 4A) 


It is clear from these graphs that the separation effect (at the 45° end — 
walls), which commences in the 220 - 230 CFS range has fully developed 
by 260 CFS, and that for higher discharges the flow across the end sill is - 
concentration into a width approximately equal to the width between the 
parallel sides of the structure. 

The P.W.D. design is less effective at small discharges but more ef- 
fective at high discharges. The end sill of this structure is 16.0 feet wide, 
compared with the 17.5 feet width of the Bradley-Peterka structure. If 
the P.W.D. structure were lengthened sufficient to give an end sill width 
of 17.5 feet, then it is probable that at low flows (i.e. not greater than 220 
CFS) the two curves would coincide. 

The Australian structure also has the advantage that the maximum 
depth of scour occurs at a moderate discharge (230 - 260 CFS) and that 
excessive overloading will not produce progressively deeper scour. 


2. Downstream Extremity of Scour (see figs. 2B, 4B) 


This aspect is of particular importance to roadway culvert outlets, as 
it is highly desirable to have the entire area effected by scour confined 
within the boundaries of the roadway reserve. 

The curves for the Bradley-Peterka structure again show the influence 
of separation, at the 45° end walls, under high discharge conditions. The 
transition range from zero separation to full separation being again 220 - 
260 CFS. 

The P.W.D. design shows a steady increase over the whole range, 
levelling off at a little over 40 feet from the end of the pipe. 

The two sets of curves corresponding to the two tailwater conditions 
are almost identical. 


3. Maximum Scour at Structure (see figs. 3A, 5A) 


The safety of the structure and of the embankment depends largely on 
this property. 

An important model alteration was made here in order to give a true 
comparison. 

The American structure has an end sill which extends down to a depth 
of 3 feet. With the erodible material used in these tests, it was found that 
after some flows of short duration, the floor of the structure under the sill 
was undercut. This produced a “model effect” scour pattern which made 
true comparison difficult. To prevent this undercutting, the walls of the 
structure were produced to a depth of 6 feet (prototype) below sill level— — 
the internal and plan dimensions remaining the same. In this way, scour 
depths to RL—6.0 could be tolerated. 2 

The Australian structure was treated similarly. ' 

The Bradley-Peterka design showed distinct advantages in this aspect, © 
though tailwater was observed to effect the relationship between the curves 
of the two structures. 


Comparative photographs of typical scour around the two structures 
are shown in figs. 9A, 9B. 
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4. Scour near toe of Embankment. (see figs. 3B, 5B) 


a 


b 
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This aspect is important from the point of view of safety of the em- 
bankment. 

The serious effects of the extension of scour in the upstream direction 
along the side walls of the structure and towards the roadway, depend on 
the following:- 


a) Maximum scour at structure. 
b) The location of deepest scour. 
c) Angle of repose of the erodible material. 


The curves for the two designs lie fairly well together —the deeper 
structure scour in the P.W.D. design being compensated by its greater 
distance away from the end of the pipe. 

The points corresponding to discharges in excess of 290 cusecs are 
doubtful as the boundary walls of the moveable bed (set 4' -6"' apart), 
caused an eddy at the toe of the embankment on each side of the structure. 
This was clearly a model effect. 


Discussion on Comparison Tests. 
The comparison tests were carried out to determine:- 


Whether the Irrigation outlet of Bradley and Peterka could be successfully 
applied to Australian road culverts. 


Whether the structure which had evolved from the research programme at 
the P.W.D. Laboratory was superior or inferior to the American struc- 
ture. 


From the viewpoint of first cost, the P.W.D. structure involves ap- 
proximately three quarters of the quantity of concrete in the corresponding 
American outlet. Also, the formwork in the Australian design is greatly 
simplified. 

As far as performance is concerned, it is very difficult to declare that 
either design is superior. 

Comparison of hydraulic models in moveable beds is usually made on 
the basis of “maximum depth of scour.” This is unreliable here as, in 
the Bradley-Peterka structure, the scour patterns within the design range 
(0 to 190 CFS) are quite different from those beyond this range. When 
tested beyond 0-190 CFS., the scour hole tends to become elongated and 
deep along the centre line. The scour holes of the two structures are 
compared for three flows in diag. 6. 

A deep scour hole downstream from a structure is not, in itself, bad 
unless it promotes erosion of material under or around the structure; the 
P.W.D. design, although causing generally a shallower scour hole, always 
has slightly more scour at the structure. 

The downstream extremity of scour, is again related to the shape of the 
scour pattern. The elongated hole which occurs when the Bradley-Peterka 
outlet is subject to high flows, is dangerous from the point of view of 
gully erosion. Small flows, when passed through a previously scoured 
elongated hole, tend to cut a channel along the centre line. This effect is 
much reduced when small flows pass through a “square ended” scour hole, 
as produced by the P.W.D. structure. 
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Scour at the structure and proximity of scour to the toe of the embank- 
ment are both important factors to be considered. The Bradley-Peterka 
structure is superior here i.e. the tendency to scour around the structure © 
is slightly less in the American than in the Australian. ‘ 

However, in both cases, rock protection would have to be carried out 
thoroughly around the structure. 

Scour along the side walls was due entirely to “slipping” of the non- 
cohesive bed material in order to lie at its appropriate angle of repose. 

The degree to which this failure occurs is due to the type of material in 
which the structure is built. If scour at the two ends of the end sill can 
be prevented or controlled by suitable rock placement, “angle of repose” 
scour will be eliminated. 

As mentioned in the preamble, the Bradley-Peterka design shows a 
rigid relationship between maximum discharge and basin dimensions—no 
alteration being made with different erodible materials downstream from 
the structure. The P.W.D. design permits flexibility in this respect. If f 
the scour velocity of the downstream erodible material is known, then the 
width of the end sill can be adjusted to ensure that this velocity is not ex- 
ceeded. This adjustment is achieved by suitably fixing the length of the 
structure—the side walls remaining at the same angle with the centre line — 
for all designs. ¢ 


th lle 
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Further Investigation. 


The present investigation, which will also include study on multiple pipe — 
culverts, flood gates, causeways and inlet structures, is far from com- 
plete. However, progress to date has been satisfactory. An attempt will 
be made to increase the angle which the side walls make with the centre 
line in order to further reduce the cost of the structure, and also to study 
the effects of different end sill shapes. 
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Bradley-Peterka basin P.W.D. Outlet Structure 
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Separation boundary 
at very high flows. 
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Diagram showing flow through Bradley- Peterka 
structure at high flows ie. Q> 260 CFS. for 54” pipe. 
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SKETCH OF P.W.D. CULVERT FOR PIPES 18" TO 72” 
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9A, 9B: Scour Around Structures often Test with High Discharge 
(is Q7260CFS) A: Bradley-Peterka Structure. B: PWD Structure. 


10A, 10B: Operation of PWD Structure with Discharge equal to 240CFS. 
A: Open Channel Flow in Pipe. B: Pipe Flowing Full. 


11A, 11B: Operation of Both Structures with Discharge equal to 330CFS. 
Pipe Flowing Full) A: Bradley-Peterka Structure. B: PWD Structure. 
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DISCHARGE CHARACTERISTICS OF RECTANGULAR THIN-PLATE WEIRS2 
Discussion by Turgut Sarpkaya 


TURGUT SARPKAYA, 1 J.M. ASCE.—The study of the characteristics of 
flow over an innocent looking weir has become a basic and challenging prob- 
lem for research workers in fluid mechanics, quite apart from the interest it 
holds for its purely practical engineering applications. Therefore, any engi- 
neer who is concerned with the measurement of discharge by the use of thin- 
plate rectangular weir along with a true understanding of the limitations of 
various approximate formulas should welcome this most noteworthy addition 
of experimental data. The precision of the authors’ measurements and the 
wide variation of the parameters involved, excluding the fluid properties, and 
the qualitative physical explanations of the observed characteristic curves of 
discharge coefficient for the raison d’étre of the proposed formula are the 
outstanding features of the paper. 

The authors associate the influence of viscosity and surface tension with a 
fictitious increase in head h, and in width b, in order to obtain a discharge 
coefficient independent of the size of the weir and of the magnitude of the 
head and embody the combined effects of the several phenomena attributed to 
viscosity and surface tension in the quantities k, and k, without fully ack- 
nowledging the limitations which these quantities imposed upon the proposed 
formula. This writer does not wish to imply that a formula developed by the 
adjustments of head and width or by any other method are not valuable as 
long as they satisfy the experimental results within the range they are de- 
veloped for. He would like, however, to see that the development of the em- 
pirical formulas be closely related with the natural variations of the parame- 
ters involved and be interpreted in that manner wherever possible. And at 
all times it should be kept in mind that an empirical equation derived by fit- 
ting a curve to a certain range of data can scarcely be expected to be general, 
or physically sound. This writer would like to remain faithful to the original 
definition of the discharge coefficient and would submit herewith a criticism 
of the authors’ presentation together with an analysis which he believes will 
clarify the problem further. 

The standard form of flow formula is based on the Bernoulli proposition of 
equality of the sum of pressure and kinetic energies in the two Sections con- 
cerned. Changes of energy due to fluid friction are allowed for only in their 
effect on the discharge coefficient, which is an experimentally determined 
value. In addition to the fluid friction, the discharge coefficient is influenced 
by, and allows for, other departures from ideal conditions, e.g., the nonuni- 
formity of velocity distribution, contraction of jet, and the mean curvature of 


a. Proc. Paper 1453, December, 1957, by Carl E. Kindsvater and Rolland W. 
Carter. 
1. Asst. Prof., Dept. of Eng. Mech., Univ. of Nebraska, Lincoln, Nebr. 
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the stream lines. The velocity distribution is determined by initial accelera~_ 
tions, fluid and boundary friction, and the geometrical characteristics of the 
weir. The increasing influence of viscosity shown by the lower Reynolds’ 
numbers influences the coefficient of discharge in three ways: (a) Internal 
fluid and boundary friction is increased, leading to greater resistance to 
shear in the acceleration through the notch, hence lowering of the discharge 
coefficient, (b) The velocity distribution of flow approaching the weir-notch ~ 
is changed in a direction which increases the central velocity in relationto 
the mean. This results in a greater total kinetic energy in the approach flow 
and a lower requirement of pressure energy to accelerate through the notch, 
thus raising the discharge coefficient, (c) Increased wall friction on the up- — 
stream face of the weir plate, causing a retardation of the layers close to the 
surface which form, by their momentum, the lower and the side boundaries 
of the nappe. This results in an enlargement of the latter quantities in com- 
parison with the ones in potential flow, hence an increased coefficient. The 
above reasoning is in close conformity with the one given by the authors. 
However, the authors have failed to set forward the limitations of their rea- 
soning. Is it true that the discharge coefficient will increase with viscosity 
and surface tension always? If it is so, then highly viscous fluids such as 
glucose, honey, and molasses would flow faster than water for a given head. 
This cannot be true. What then is the answer to the apparent reversal of the 
foregoing reasoning? The only explanation for this is that for values of 
Reynolds’ number smaller than a certain value the discharge coefficient de- 
creases rapidly. This is already a well-known phenomenon for flow in ori- 
fices, nozzles and venturi tubes, (Fig. 1). The humped shape of orifice co- 
efficient curves is quite comparable with rectangular notchweir discharge 
coefficient curves. Obviously, the predominance of the influences (b) and (c) 
over a limited flow phase for Reynolds’ numbers greater than Rer, creates 
the hump on the coefficient curve. For values of R smaller than Rey, the in- 
fluence (a) is predominant. In the case of a low ratio of area with practically 
zero approach velocity the hump is due to influence (c) alone. Since the flow 
over the weir is free surface flow, the surface tension has additional effect 
on the influences sighted above and makes the analysis of the problem even 
more complicated. However, since the viscosity seems to play a more im- 
portant role, for lower values of Reynolds’ number—or in other words as the 
viscosity increases—the coefficient of discharge will decrease instead of in- 
crease which would be expected. Therefore, the authors’ explanations over 
the influences of viscosity and surface tension should definitely be confined 
to a flow range where the Reynolds’ number is greater than a certain critical © 
value. Moreover, for certain values of Reynolds number some irregularities © 
can be expected in discharge coefficient due to the transitional instability of 
the flow in the approach channel at the transition between turbulent and vis- 
cous flows. The proof for the existence of a critical Weber number and its 
experimental verification is rather difficult. Firstly because, for the liquids 
commonly used, the one with a higher viscosity has a lower surface tension . 
at a given temperature, and secondly because for a given liquid the increase 
of surface tension with a decrease in temperature, is much slower than the 
increase of viscosity. Therefore, the viscosity becomes always predominant 
and makes it very difficult to obtain a wide range of Weber numbers when the 
viscosity and the flow characteristics are kept constant. | 
. 
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The equation proposed by the authors as a comprehensive equation of dis- 
charge for all rectangular thin-plate weirs could be written in the following 
form: 
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2 3/2 k 
Q=—/2c bh eee eye (1) 
3 ‘ b h 


in which C, is assumed to be a function of b/B and h/P alone. If the last 
parentheses is opened and higher orders of small quantities neglected, one 
obtains, 


Q ig, of, 
PSE ASIEN Fe - (1+ +b yCe (2) 
2 fe bh > 1g 
3 


Experiments of various investigators seem to establish the fact that Co 
could be written in the following functional form: 


C= £(b/B) + g(b/s) ~ (3) 


in which f and g are used as a symbol of “the function of.” This writer pro- 
poses that f and g be expressed in the following form: 


£(b/2) = OL+ B-(b/2) pe 


and 


e(b/s) = ¥ +8.) (4-b) 


in whicha,8, ”,46,m, and n are assumed to be constants independent both 
of the characteristics of the weir and of the properties of fluid. Hence, the 
coefficient of discharge becomes, 


4 * se) -{ep on +]¥+8 o/s) | +} | Tat sie er = +s} (5) 


b 


The mean velocity at the section of weir notch is given by, 
yes c.¥2gh (6) 
3 


The derivation of a general equation that should be applicable to any 
liquid flowing over a rectangular notch weir is aided by a dimensional analy- 
sis of the problem. Although, the Eq. (5) is empirical one, it should remove 
some of the objections to its use for liquids other than those actually tested 
if it is expressed in terms of Reynolds and Weber numbers and verified later 
by experiments. The coefficient of discharge Cg which is defined here as in- 
dependent of the acceleration of gravity is, 


ita =. se 


1616-96 HY 2 April, 1958 


Q 
C see Tel Del eh a, aE eo (7) 
a 2/ bh /2gh ie 


ee ee 


The influence of the ratio b/h will be neglected in view of the experimental 
results mentioned by the authors. The definitions of the Reynolds and Weber 
numbers for this particular problem is somewhat crucial. Since both the 
Reynolds and Weber number must be related to both the head and the width of 
the weir, and since there is no way of predicting which length is physically 
more significant for each number. In each’number, velocity can be replaced — 
by a dimensional equivalent /gh. However, it should be kept in mind that if — 
the particular function of Eq. (6) were known, this substitution would not be 
entirely correct. But since the particular function is not known, the replace- 
ment is possible. In this paper Reynolds and Weber numbers will be defined 
as follows: ¢ 


R«—— and We ‘ae (8) ; 
v (° 7 


: 
If now b/B is replaced by 7,,h/P by 7, 1/R by 773, and 1/W by 74, one has, 


c. ae (Tt, ,T1,,11,,11, ) (9) f 


Dimensional analysis states that the general solution of the problem is in the 
form, 
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in| which) P5768; t5.q, ss steske are constants to be determined. Letting a be 
the value of Cq as 71, 79, 73, and 74 approach zero, one has | 
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If Ey. (11) is compared with Eq. (5), one obtains, 


c, -{fe+s (b/B) ]+[¥- § (o/s) | +4 {. a a i =| (12) 
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low the constants in Eq. (4-a) and (4-b) are obtained by plotting the values of 
b/B) and g(b/B) versus b/B, obtained by the authors’ Fig. (10). A logarith- 
lic plot of these values gives a perfect straight line. However, here only 

1e plots of these values on an ordinary graph paper are shown, Fig. (2). 

he values obtained are given| below: 


L- 0.5g $= 0.0773 
Bs 09.015 Me 2.43 
Y= -0.0023 N= 2.43 
ubstituting, 
~6 
Ve nxo rt*/sec, O = 0.005 it/rt, 
kK = 0.003 ft PH 1.94 slug/tt? 


quation (12) takes the following final form: 


243 260 
i {fs + 0.015 (b/B) ]+[ 0.0023 + 0.773 (b/B) a +} x 


vac tp ams be 
/ 2 


ae 
Ro re W 


(14) 


his equation holds very good for all values of b/B and h/P except for b/B = 
.9. For this latter value experimental results gives about 1% greater Cg 
alues. In the above equation the value of & = 0.587 deserves more attention. 
he theoretical value of o for slots, orifices, sluice gates and rectangular 
eirs is 0.611. 

With the definitions accepted for Reynolds and Weber numbers, the varia- 
on of kp and kp with the fluid properties can be obtained from, 


oe Ge <3 2B (C 7 
ne ” Sow G6) Siena ie ae 


) which the subscripts x and w indicate respectively “for any liquid” and for 
ater.” It would be indeed very interesting to determine experimentally the 
ariation of k, and k, with viscosity and surface tension. If the Eqs. (15-a,b) 
old for various liquids tested, then the Eq. (14) could be used for all liquids. 
nd only then the figures (11) and (12) presented by the authors can have a 
-ue and broader meaning. Moreover, the usage of highly viscous fluids will 
smove the uncertainties existing in kp and particularly in kp. Since accord- 
1g to Eq. (15-a) the value of kpx for a liquid 8 times more viscous than wa- 
1r would be 0.012 ft., which is a value sufficiently great to determine more 
ecurately by means of experiments. It should be recognized that, after the 
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demonstration of the use of R and W as suitable parameters, the treatment 
is purely empirical. It is conceivable that several different functions mig’ 
exist which would satisfy the experimental data equally well. Extrapolation 
of these functions outside the experimental range is unwarranted and shou 
be made with caution or not at all. 

In conclusion, this writer hopes that further experiments will be per- 
formed with different liquids to make the proposed equation applicable to 
liquids other than water as well as water, and held over a wide range of b, 


h/P, R, and W. 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are identified by number below, Technical-division 
ansorship is indicated by an abbreviation at the end of each Paper Number, the symbols referring to: Air 
ansport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
tulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (SA), Soil 
chanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and Harbors - 
W), divisions. Papers sponsored by the Board of Direction are identified by the symbols (BD). For 
ies and order coupons, refer to the appropriate issue of “Civil Engineering,” Beginning with Volume 82 
muary 1956) papers were published in Journals of the various Technical Divisions. To locate papers in 
Journals, the symbols after the paper numbers are followed by a numeral designating the issue of a par- 
ular Journal in which the paper appeared. For example, Paper 1449 is identified as 1449 (HY 6) which 
jieaies that the paper is contained in the sixth issue of the Journal of the Hydraulics Division during 1957. 


: ° VOLUME 83 (1957) 


RIL: 1196(EM2), 1197(HY2), 1198(HY2), 1199(HY2), 1200(HY2), 1201(HY2), 1202(HY2), 1203(SA2), 1204 
), 1205(SM2), 1206(SM2), 1207(SM2), 1208(WW1), 1209(WW1), 1210(WW1), 1211(WW1), 1212(M2), 
13(EM2), 1214(EM2), 1215(PO2), 1216(PO2), 1217(PO2), 1218(SA2), 1219(SA2), 1220(SA2), 1221(SA2), 
> 1223(SA2), 1224(SA2), 1225(PO)c, 1226(WW1)¢, 1227(SA2)°, 1228(SM2)°, 1229(EM2)°, 1230 
(BY2)°. 


t¥:; 1231(ST3), 1232(ST3), 1233(ST3), 1234(ST3), 1235(IR1), 1236(IR1), 1237(ww2), 1238(wWw2), 1239(Ww2), 
124c(ww2), 1241(ww2), 1242(Www2), 1243(WW2), 1244(HW2), 1245(HW2), 1246(HW2), 1247(HW2), 1248 
(Ww 2), 1249(HW2), 1250(HW2), 1251(Ww2), 1252(WW2), 1253(IR1), 1254(ST3), 1255(ST3), 1256(HW2), 1257 
(@R1)°, 1258(HWa2)c, 1259(ST3)°, 


NE: 1260(HY3), 1261(HY3), 1262(HY3), 1263(HY3), 1264(HY3), 1265(HY3), 1266(HY3), 1267(PO3), 1268 

(PO3), 1269(SA3), 1270(SA3), 1271(SA3), 1272(SA3), 1273(SA3), 1274(SA3), 1275(SA3), 1276(SA3), 1277 

(HY3), 1278(HY3), 1279(PL2), 1280(PL2), 1281(PL2), 1282(SA3), 1283(HY3)¢, 1284(PO3), 1285(PO3), 1286 
), 1287(P03)¢, 1288(SA3)°. 


LY: 1289(SM3), 1290(EM3), 1291(EM3), 1292(EM3), 1293(EM3), 1294(HW3), 1295(HW3), 1296(HW3), 1297 
3), 1298(HW3), 1299(SM3), 1300(SM3), 1301(SM3), 1302(ST4), 1303(ST4), 1304(ST4), 1305(SU1), 1306 
1), 1807(SU1), 1308(ST4), 1309(SM3), 1310(SU1)°, 1311(EM3)°, 1312(ST4), 1313(ST4), 1314(ST4), 1315 
4), 1316(ST4), 1317(ST4), 1318(ST4), 1319(SM3)°, 1320(ST4), 1321(ST4), 1322(EM3), 1323(AT1), 1324 
1), 1325(AT1), 1326(AT1), 1327(AT1), 1328(AT1)°, 1329(ST4)¢. 


| : 1330(HY¥4), 1331(HY4), 1332(HY4), 1333(SA4), 1334(SA4), 1335(SA4), 1336(SA4), 1337(SA4), 1338 
, 1389(CO1), 1340(COL), 1341(CO1), 1342(CO1), 1343(CO1), 1344(PO4), 1345(HY4), 1346(PO4)¢, 1347 
BD1), 1348(HY4)¢, 1349(SA4)° 1350(P04), 1351(PO4), 


PTEMBER: 1352(IR2), 1353(ST5), 1354(ST5), 1355(STS), 1356(ST5), 1357(ST5), 1358(ST5), 1359(TR2), 1360 
(@R2), 1361(ST5), 1362(0R2), 1363(1R2), 1364(IR2), 1365(WW3), 1366(WW3), 1367(WW3), 1368(WW3), 1369 
(Ww3), 1370(ww3), 1371(HW’4), 1372(HW4), 1373(HW4), 1374(HW4), 1375(PL3), 1376(PL3), 1377(IR2)¢,1378 
carey, 1379(IR2), 1380(HW4), 1381(WWw3)°, 1382(ST5)°, 1383(PL3)°, 1384(IR2), 1385(HW4), 1386(HW4), 


R: 1387(CP2), 1388(CP2), 1389(EM4), 1390(EM4),1391(HY5), 1392(HY5), 1393(HY5), 1394(HY5), 1395 
), 1396(PO5), 1397(PO5), 1398(PO5), 1399(EM4), 1400(SA5), 1401(HY5), 1402(HY5), 1403(HY®), 1404 
v5), 14050H¥5), 1406(HY5), 1407(SA5), 1408(SA5), 1409(SA5), 1410(SA5), 1411(SA5), 1412(EM4), 1413 

4), 1414(PO5), 1415(BM4)©, 1416(POS)°, 1417(HY¥5)°, 1418(EM4), 1419(PO5), 1420(P05), 1421(PO5), 
(SA5)°, 1423(SA5), 1424(EM4), 1425(CP2), 


R: 1426(SM4), 1427(SM4), 1428(SM4), 1429(SM4), 1430(SM4)°, 1431(ST6), 1432(ST6), 1433(STS), 
(ST6), 1435(ST6), 1436(ST6), 1437(ST6), 1438(SM4), 1439(SM4), 1440(ST6), 1441(ST6), 1442(ST6)¢, 
(SU2), 1444(SU2), 1445(SU2), 1446(SU2), 1447(SU2), 1448(SU2)°. 


(BER: 1449(HY6), 1450(HY6), 1451(HY6), 1452(HY6), 1453(HY6), 1454(HY6), 1455(HY6), 1456(HY6)¢, 
457(PO6), Paseo, 1459(PO6), 1460(PO6)°, 1461(SA6), 1462(SAG), 1463(SAG), 1464(SA6), 1465(SA6), 
S(SAG)©, 1467(AT2), 1468(AT2), 1469(AT2), 1470(AT2), 1471(AT2), 1472(AT2), 1473(AT2), 1474(AT2), 
2), 1476(AT2), 1477(AT2), 1478(AT2), 1479(AT2), 1480(AT2), 1481(AT2), 1482(AT2), 1483(AT2), 
72), 1485(AT2)°, 1486(BD2), 1487(BD2), 1488(PO6), 1489(P06), 1490(BD2), 1491(BD2), 1492(HY6), 


1493(BD2). 
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AF + 1494(BM1), 1495(EM1), 1496(EM1), 1497(IR1), 1498(IR1), 1499(1R1), 1500(IRL), 1501(IR1), 1502 
(IR1) J Mivcsey at on 1505(IR1), 1506(IR1), 1507(IR1), 1508(ST1), 1509(ST1), 1510(ST1), 1511(ST1), 
12(ST1), 1513(WW1), 1514(ww1), 1515(WW4), 1516(WW1), 1517(WW1), 1518(WW1), 1519(ST1), 1520 
M1)°, 1521(1R1)°, 1522(ST1)°, 1523(WW1)°, 1524(HW1), 1525(HW1), 1526(HW1)°, 1527(HW1). 

au + 1528(HY1), 1529(PO1), 1530(HY1), 1531(HY1), 1532(HY1), 1533(SA1), 1534(SA1), 1535(SM1), 
[53 eu) Gest teuis: avai 1539(SA1), 1540(SA1), 1541(SA1), 1542(SA1), 1543(SA1), 1544(SM1), 
5(SM1), 1546(SM1), 1547(SM1), 1548(SM1), 1549(SM1), 1550(SM1), 1551(SM1), 1552(SM1), 1553(P01), 
554(PO1), 1555(PO1), 1556(PO1), 1557(SA1)°, 1558(HY1)°, 1559(SM1)°, 

RCH: 1560(ST2), 1561(ST2), 1562(ST2), 1563(ST2), 1964(ST2), 1565(ST2), 1566(ST2), 1567(ST2), 1568 
(wwe), re sped Np 1571(WW2), 1572(WW2), 1573(WWw2), 1574(PL1), 1575(PL1), 1576(ST2)°, 
(577(PL1), 1578(PL1)°, 1579(wwe2)°. 

HL: 1580(BM2), 1581(EM2), 1582(H¥2), 1583(HY2), 1584(HY2), 1585(HY2), 1586(HY2), 1587(HY2), 1588 
1Y2), atin ee sears): 1592(SA2), 1593(SU1), 1594(SU1), 1595(SU1), 1596(EM2), 1597(PO2), 
598(PO2), 1599(PO2), 1600(PO2), 1601(PO2), 1602(PO2), 1603(HY2), 1604(EM2), 1605 (SU1)°,1606(SA2), 
60 (SA2), 1608(SA2), 1609(SA2), 1610(SA2), 1611(SA2), 1612(SA2), 1613(SA2), 1614(SA2)°, 1615(1R2)°, 1616 
¥2)¢, 1617(SU1), 1618(PO2)°, 1619(EM2)°, 1620(CP1). 


of several papers, grouped by divisions. 
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